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1.0  SUMMARY  OF  SIGNIFICANT  ACCOMPLISHMENTS 

The  significant  accomplishments  from  this  research 
project  fall  into  five  areas. 

1.1  Real  time  observations  in  the  scanning  electron 
microscope  have  delineated  the  details  of  the  fracture 
processes  that  result  in  mode  I  and  mode  II  delamination  of 
composite  materials.  These  observation  give  clear 
explanation  of  why  Guc//Gic  rat*°  f°r  brittle  materials  is 
3.0  or  larger  while  for  ductile  materials  it  is  closer  to  1.1 
These  in-situ  observations  of  fracture  also  indicate  that 
distinctively  different  damage  zones  develop  ahead  of  growing 
cracks  for  mode  I  and  mode  II  delamination . 

1.2  Two  techniques  have  been  developed  to  measure  the 
strain  field  around  a  crack  tip.  Stereo-imaging  and  direct 
measurement  of  distortion  of  a  fine  array  of  dots  placed  on 
the  surface  have  both  been  found  to  be  effective  in  measuring 
the  strain  field  around  a  crack  tip.  One  of  the  surprising 
results  from  these  measurements  is  that  the  local  strain  to 
failure  at  the  crack  tip  is  much  greater  than  the  elongation 
measured  in  a  tensile  test  (48%  versus  8%  for  Hexcel  F185). 

1.3  A  linear,  orthotropic  finite  element  code  has  been 
used  to  calculate  the  stress  fields  around  the  crack  tip  for 
mode  I  and  mode  II  loading.  The  relative  size  and  shape  of 
the  stress  fields  for  these  two  different  loading  conditions 
explains  the  different  sizes  and  shapes  of  the 
damage/deformation  zones  see  in-situ  (as  noted  in  1.1).  A 
phenomenological  model  has  been  developed  to  try  to  explain 
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why  Glc  and  GIIc  for  delamination  of  a  composite  material 
made  from  a  tough  resin  are  so  much  lower  than  the  neat  resin 

GIc ' 

1.4  A  J-integral  approach  for  mode  I  has  been  used  to 

investigate  the  delamination  of  multi-directional  composites. 
A  comparison  of  the  end  notch  flexure  test  and  the  loaded 
flexure  test  for  measuring  Gllc  has  found  that  they  give 
consistent  results.  A  J-integral  approach  for  the 

measurement  of  for  composites  made  from  more  ductile 

resin  has  been  developed.  Initial  results  using  this 
approach  look  very  promising  for  characterizing  systems  which 
previously  could  not  be  characterized  with  conventional 
analyses . 

1.5  A  combination  of  in-situ  fracture  observations  with 
postmortem  fractography  has  clarified  some  long  standing 
questions  with  regard  to  f ractographic  interpretation  in 
failure  analysis. 

2.0  OBJECTIVES 

The  objectives  of  the  research  program  whose  results 
are  summarized  in  this  report  have  been: 

2.1  to  better  define  the  deformation  and  fracture  physics 
of  delamination  fracture  in  graphite/epoxy  composite 
materials  so  that  realistic  micromechanics  models  of  matrix 
dominated  fracture  can  be  developed;  and 

2.2  to  develop  reliable  experimental  and  analytical 
techniques  to  measure  mode  I,  mode  II,  and  mixed  mode 
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delamination  fracture  toughness  of  both  unidirectional  and 
mul tidi rectional  composite  laminates  to  provide  meaningful 
design  parameters  and  benchmarks  against  which  predictions  of 
the  various  micromechanics  models  may  be  tested. 

3.0  SUMMARY  WORK  STATEMENT 

The  following  areas  of  work  have  been  performed  to 
achieve  the  stated  objectives. 

3.1  In-situ  fracture  studies  in  a  scanning  electron 
microscope  (SEM)  have  been  undertaken  to  better  define  the 
deformation  and  fracture  processes  which  accompany  mode  I  and 
mode  II  delamination  fracture. 

3.2  Two  techniques  to  measure  the  displacement  fields 
around  crack  tips  in  neat  resins  and  composites  from  direct 
observations  made  in  the  scanning  electron  microscope  have 
been  developed  and  utilized. 

3.3  A  simple  phenomenological  model  consistent  with 
observations  made  in  3.1  and  3.2  which  is  capable  of 
predicting  the  results  obtained  in  3.4  is  being  developed. 

3.4  Testing  and  analytical  techniques  have  been  developed 
for  characterizing  mode  I,  mode  II,  and  mixed  mode 
delamination  of  split  laminates  which  give  very  nonlinear 
load-displacement  curves;  e.g.  ones  made  from  more  ductile 
resins  and/or  split  laminates  with  multi-directional  layups. 

3.5  The  results  from  the  in-situ  observations  of  fracture 
in  the  SEM  have  been  utilized  in  conjunction  with  post-mortem 
fractography  to  develop  a  failure  analysis  methodology. 
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4.0  SUMMARY  OF  ACCOMPLISHMENTS  BY  AREAS  STATED  IN  WORK 

STATEMENT 

4.1  In-Situ  Fracture  Observations 

A  summary  of  significant  results  will  be  given  in  this 
section  with  more  complete  results  and  discussion  available 
in  the  twelve  papers  that  have  been  published  based  on  this 
work.  These  papers  are  listed  in  section  5.0  with  copies  of 
each  in  Appendix  II. 

In  this  research  program,  in-situ  fracture 
observations  have  been  made  for  both  neat  resin  loaded  in 
mode  I  and  graphite/epoxy  composite  materials  loaded  to  give 
mode  I  and  mode  II  delamination.  The  primary  objective  of 
these  observations  have  been  to  determine  the  crack  tip 
processes  that  lead  to  fracture  and  to  relate  the  size  and 
shape  of  the  damage  zone  observed  around  the  tip  of  a  growing 
crack  to  the  macroscopically  measured  fracture  toughness. 

4.1.1  In-Situ  Observations  in  SEM  of  Fracture  of  Neat  Resins 

In-situ  observations  of  neat  resin  fracture  have  been 
made  on  standard  0 . 5T  compact  tension  specimens  (see 
ASTM-E399 ^ ) ,  except  that  the  thicknesses  were  typically  about 
3-4  mm  rather  than  the  specified  12.5  mm. 

Typical  in-situ  fracture  observations  are  presented  in 
Fig.  1  for  the  relatively  brittle  resin  Hercules  3502  (highly 
crosslinked  epoxy)  and  for  the  much  more  ductile  Hexcel  F185 
(moderately  crosslinked  epoxy  with  14%  rubber  added  to 
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enhance  toughness).  The  very  fine  microcracking  observed  in 
these  figures  is  due  either  to  deformation  and  fracture  of 
the  150A  thick  gold/palladium  film  that  is  sputter  coated 
onto  the  polished  surface  of  the  neat  resin  compact  tension 
specimens  to  avoid  charging  or  real  microcracking  in  the 
specimen.  It  has  been  determined  that  the  film  microcracking 
begins  at  about  3-5%  tensile  strain.  (This  was  done  by 
coating  a  rubber  band  and  stretching  it  in  the  SEM.)  Thus, 
the  microcracked  zone  gives  some  indication  of  the  volume  of 
material  surrounding  the  crack  tip  which  experiences  either 
significant  deformation  (at  least  3%)  or  real  resin 
microcracking . 

F185  Resin — The  "butterfly  wings"  pattern  of  the 
microcrack  distribution  as  well  as  the  very  fine  scale  of 
cracking  in  the  in-situ  observations  of  fracture  of  F185  neat 
resin  (Fig.  1)  indicate  this  to  be  only  film  cracking 
resulting  from  the  considerable  plastic  deformation  occurring 
in  this  very  ductile  resin  prior  to  crack  advance.  A 
comparison  of  the  fracture  surface  of  the  F185  observed  post 
mortem  to  the  fracture  surface  of  the  3502  which  fails 
brittlely  clearly  indicates  that  the  F185  experiences  a 
ductile  fracture  (see  Fig.  2).  This  further  substantiates 
our  contention  that  the  microcracking  observed  in-situ  in  the 
F185  is  coating  microcracking  due  to  resin  deformation. 

The  viscoelastic  nature  of  the  Hexcel  185  resin  gives 
a  very  discontinuous  crack  growth  behavior,  as  seen  in  Fig. 3. 
The  initial  loading  is  occurring  under  displacement  control 
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Figure  1 


v 


In-situ  fracture  of  (a)  3502  (4300X)  and  (b) 
F185(1000X).  The  white  region  along  the  crack 
is  due  to  charging  of  uncoated  resin,  exposed 
by  the  fracture  process.  The  point  lines 
especially  evident  in  the  3502  arc  polishing 
scratches  introduced  during  specimen 
preparation . 
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at  a  very  slow  rate.  Crack  tip  deformation  and  blunting  are 
noted,  with  the  deformation  zone  indicated  by  the 
microcracking  of  the  100A  thick  gold/palladium  film.  Once 
the  strain  in  the  crack  tip  reaches  a  critical  value,  a  sharp 
crack  emerges  from  the  blunted  notch.  The  material  is  unable 
to  support  this  slightly  longer  and  much  sharper  crack  at  its 
current  load  level.  Furthermore,  the  blunting  process 
previously  observed  requires  some  time  to  be  realized,  time 
which  was  formerly  provided  by  the  slow,  monotonic  loadup. 


With 

the  specimen 

now  already  at 

full  load,  the 

crack 

will 

tend 

to  propagate 

unstably  and 

quickly  with 

much 

less 

deformation  until  the  specimen  unloads  sufficiently  to  give 
crack  arrest.  For  slow  monotonic  loading  in  displacement 
control,  this  process  of  crack  tip  blunting,  rapid  advance, 
and  arrest  will  be  repeated  many  times. 

Post-mortem  fractography  of  the  F185  specimens 
indicates  a  considerable  difference  in  the  amount  of 
deformation  that  occurs  during  unstable  crack  propagation  as 
compared  to  reinitiation  of  crack  growth  during  monotonic 
loading  (in  displacement  control).  The  individual 
"platelets1'  seen  in  Fig.  2  represent  unstable  advance  of  the 
crack  tip.  (Compare  the  length  of  unstable  crack  advance  in 
Fig.  3  to  platelet  size  in  Fig.  2.)  The  steps  between 
"platelets"  are  the  result  of  the  crack  tip  blunting  that 
occurs  during  reloading  to  give  further  crack  advance. 
Ductile  fracture  is  indicated  for  both,  but  the  amount  of 
deformation  is  much  less  during  rapid  crack  advance.  The 
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Figure  3.  Microcracking  ahead  of  crack  tip  of  F 1 85  resin  showing  crack  growth 
sequence  under  constant  applied  displacement,  (a)  3000x  Detail  of 
crack  tip  blunt  crack,  (b)  9000x.  Crack  propagation  through 

?;aCT™k  Coaiescence  t0  form  a  very  Sharp  narrow  crack,  (c)  5000x 
Crack  extension  occurred.  FJew  crack  tip  is  very  sharp. 
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fine  holes  are  due  to  cavitation  of  the  rubber  particles 
added  to  the  F185  resin  to  enhance  toughness. 

3502  Resin — The  microcracking  in  the  brittle  3502  is 
much  more  localized  around  the  crack  tip,  and  the  observation 
of  these  microcracks  in  real  time  indicates  that  they 
coalesce  to  give  macrocrack  advance,  as  seen  in  Fig.  2.  With 
a  monotonic  increase  in  mouth  opening  displacement, 
microcracking  around  the  crack  tip  develops  which  partially 
shields  the  crack  tip,  reducing  locally  the  crack  tip 
stresses.  However,  once  these  microcracks  coalesce  giving 
local  crack  tip  advance,  unstable  crack  growth  follows  due  to 
the  presence  of  the  new  sharp  crack  tip  without  the  benefit 
of  microcrack  shielding.  Advance  proceeds  under  displacement 
control  until  sufficient  unloading  occurs  that  the  running 
crack  arrests.  Subsequent  quasi-static,  monotonic  increase 
in  the  mouth  opening  displacement  again  causes  the  formation 
of  new  microcracks  surrounding  the  crack  tip,  with  the 
density  increasing  with  mouth  opening  displacement  until 
coalescence  occurs,  again  giving  additional,  unstable  crack 
extension . 

The  river  pattern  seen  in  Fig.  2  is  due  to  coalescence 
with  growth  of  the  many  microcracks  that  develop  in  the 
vicinity  of  the  crack  tip  during  quasi-static,  monotonic 
loading.  Once  coalescence  has  occurred,  the  crack  advances 
rapidly  to  the  new  arrest  location,  with  little  opportunity 
for  new  microcrack  formation.  The  river  patterns  that  form 
at  the  arrest  lines  as  well  as  the  almost  smooth  fracture 


surface  in  between  arrest  lines  strongly  supports  the 
scenario  just  presented  for  crack  growth  in  brittle  resin 
systems . 

Comparison  of  3502  and  F185  Resins — The  relative 
deformation/damage  zone  sizes  seen  in  Fig.  1  for  F185  and 
3502  resins  may  be  better  understood  by  comparing  the 
constitutive  behavior  of  the  two  materials  as  determined  in 
tensile  tests,  the  results  of  which  are  presented  in  Fig.  4. 
It  is  clear  that  the  lower  yield  strength  of  the  F185  as  well 
as  the  much  more  nonlinear  constitutive  behavior  combine  to 
give  a  much  larger  deformation  zone  and  much  greater  load 
redistribution  away  from  the  crack  tip,  postponing  the 
fracture  event.  The  result  is  a  much  greater  resistance  to 
crack  propagation,  and  thus,  a  much  greater  mac roscopi cal ly 
measured  fracture  toughness. 

It  should  be  emphasized  that  the  strain  to  failure 
measured  in  tensile  tests  greatly  underestimates  the  local 
strain  to  failure  in  the  crack  tip  region.  For  example,  the 
F185  resin  gave  8-9%  elongation  in  a  2.5  cm  gage  length.  By 
contrast,  the  local  strain  to  failure  as  measured  in  a  10 
micron  gage  length  at  the  tip  of  the  crack  indicated  48% 
strain  (as  will  be  explained  in  more  detail  in  section  4.2  of 
this  report).  In  similar  fashion  the  tensile  elongation  in  a 
2.5  cm  gage  length  of  3502  resin  was  only  1.5%  whereas  the 
elongation  measured  in  a  10  micron  gage  length  at  the  tip  of 
a  mode  I  delaraination  crack  was  15%.  However,  to  the  degree 
that  crack  tip  strain  to  failure  is  some  constant  times  the 
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Figure  4.  Longitudinal  Stress-stra in  curve  of  F 1 85  and  3501-6  epoxy 


tensile  elongation,  or  more  likely  times  the  nonlinear 
component  of  the  tensile  elongation,  tensile  data  may  still 
be  useful  to  rank  resins,  though  the  tensile  elongation  would 
obviously  not  be  a  useful  fracture  parameter  in  a 
micromechanics  model. 

4.1.2  IN— SITU  OBSERVATIONS  IN  SEM  OF  DELAMINATION  FRACTURE 

OF  COMPOSITE  MATERIALS 

Direct  observation  of  the  delamination  fracture 
process  in  composite  materials  allows  one  to  determine  the 
micromechanical  processes  responsible.  For  example, 
delamination  crack  growth  may  occur  by  interfacial  debonding 
or  by  resin  cracking,  as  seen  in  Fig.  5.  Furthermore,  the 
size  and  shape  of  the  deformation/damage  zone  is  seen  to  be 
quite  variable,  depending  on  both  the  resin  toughness  and  the 
state  of  stress  (i.e.,  mode  I,  mode  II,  or  mixed  mode). 

Figure  6  presents  in-situ  fracture  observations  of  the 
deformation/damage  zone  sizes  for  mode  I  and  mode  II 
delamination  fracture  of  brittle  AS4/3501-6.  The  extent  of 
the  damage  zone  ahead  of  the  crack  tip  is  much  greater  for 
the  mode  II  loading  than  for  the  mode  I  loading.  A  much 
slower  decay  in  the  stress  field  ahead  of  the  crack  tip  (to 
be  reported  in  section  4.3  of  this  document)  for  mode  II 
loading  than  for  mode  I  loading  is  responsible  for  the 
difference  in  the  respective  damage  zone  sizes. 

The  deformation/damage  zone  for  mode  I  and  mode  II 
delamination  fracture  of  the  much  more  ductile  T6T145/F185  is 
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Figure  5 
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Crack  extension  by  resin  fracture  and 
interfacial  debond  are  seen  in  three 
composites  made  with  ductile  resins. 
AS41Q6,  Dow  Chemical  (1100X)1;  (b) 
T6T145/F185,  Hexcel  (1000X),  and  (c) 
T6145/F155,  Hexcel  (1200X). 
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seen  in  Fig.  7.  It  is  obvious  that  the  deformation/damage 
zone  for  both  mode  I  and  mode  II  delamination  of  the 
T6T145/F185  is  much  greater  than  that  observed  in  the  more 
brittle  AS4/3501-6.  A  schematic  representation  of  the 
relative  damage  zone  sizes  for  these  two  composite  systems, 
loaded  in  mode  I  and  mode  II,  is  seen  in  Fig.  8. 

The  deformation/damage  zone  size  has  been  measured  for 
several  systems  and  an  effort  has  been  made  to  see  how  the 
size  of  the  deformation/damage  zone  correlates  with  the 
delamination  fracture  toughness.  The  results  are  presented 
in  Table  1.  In  general,  larger  deformation/damage  zone  sizes 
correlate  with  higher  delamination  fracture  toughness; 
however,  the  correlation  is  not  monotonic. 

Hicromechanisms  of  Mode  I  Delamination--The  details  of 
the  fracture  process  have  also  been  noted  during  in-situ 
observations  of  fracture.  In  most  systems  studied,  mode  I 
crack  advance  occurs  after  the  development  of  a 
deformation/damage  zone  ahead  of  the  crack  tip.  In  brittle 
systems,  this  zone  consists  of  microcracks  which  coalesce 
with  each  other  and  the  macrocrack  to  give  crack  advance. 
Sometimes,  crack  advance  occurs  by  debonding  before 
microcrack  coalescence  occurs.  When  this  takes  place,  fiber 
bridging  and  fiber  breakage  will  be  observed.  In  ductile 
systems,  local  deformation  around  the  crack  tip  preceeds 
crack  advance,  but  with  no  true  specimen  microcracking,  only 
the  superficial  coating  cracks  which  are  indicative  of  large 
deformation.  Crack  advance  generally  occurs  by  ductile 


15 


Figure  6.  Damage  zones  ahead  of  the  crack  tip  of  AS4/ 3501-6  under  (a) 
mode  I  loading,  and  (b)  mode  II  loading. 
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igure  7.  Damage  zone  anead  of  crack  tip  of  T6T145/FTS5  unaer  (a)  750x ,  mode 
I  loading,  (b)  luOOx,  mode  II  loading,  (c)  IOOOx,  mode  II  loading. 
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Schematic  of  relative  damage  zone  sizes  for  mode 
3502  neat  resin  and  ASM  ,'3502  composite,  (b)  F T 85 
T6 T 1 „  p  1 85  composite. 
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tearing,  but  crack  advance  by  interfacial  debonding  was  also 
commonly  observed  in  composites  made  with  more  ductile 
resins  . 


i 


Micromechanisms  of  Mode  II  Delamination--Mode  II 
delamination  of  brittle  systems  occurs  in  a  very  distinctive 
way.  Microcracks  form  for  a  considerable  distance  ahead  of 
the  crack  tip  at  a  45°  angle  to  the  plane  of  the  ply,  which 
is  the  plane  pe rpendi cul a r  to  the  principal  normal  stress. 
These  cracks  grow  until  they  reach  the  fibers  which  bound  the 


resin 

rich 

region 

between  plies 

It  sometimes  appears  that 

they 

stop 

in 

the 

resin  short 

of  the  nearest  visible  fiber 

(  see 

Fig . 

6, 

for 

example ) . 

However,  there  are  certainly 

fibers  just  beneath  the  surface  which  are  responsible  for 

arresting  these  growing  microcracks.  Coalescence  of  these 

microcracks  is  required  for  macrocrack  advance.  This 

coalescence  generally  occurs  at  the  fiber/matrix  interface, 

giving  a  corrugated  roof  appearance  along  with  the  formation 

of  "hackles"  in  the  resin  between  fibers,  as  seen  in  Fig.  9. 

It  is  apparent  that  the  final  orientation  of  the  microcracks 

and  the  resulting  hackles  is  greater  than  45°.  This  is  a 

result  of  the  fact  that  hackles  aie  often  rotated  into  a  more 

nearly  upright  position  by  the  mode  II  loading  prior  to  final 

microcrack  coalescence,  as  shown  schematically  in  Fig.  10. 

It  has  been  noted  by  several  investigators  (e.g.,  Bradley  and 
2 

Cohen  )  that  the  GTT_/G_  ratio  for  delamination  of 

1 1  c  I  c 

composites  made  with  brittle  resins  is  typically  3  or 
greater.  This  is  understandable  when  one  considers  the 
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following  factors:  (1)  a  much  more  extensive  damage  zone 
ahead  of  the  crack  tip;  (2)  the  greater  difficulty  in  getting 
microcrack  coalescence;  (3)  the  more  torturous  path  taken  by 
the  crack;  and  (4)  the  much  greater  evidence  of  crack 
interaction  with  the  fibers  via  the  microcracks  which  form 
and  arrest  temporarily  at  the  fibers. 

Mode  II  delamination  of  the  more  ductile  T6T145/F185 
gave  no  evidence  of  resin  microcracking.  Rather,  only  a  very 
fine  coating  cracking  was  noted  (see  Fig.  7),  indicating 
resin  deformation.  Thus,  the  fracture  process  for  mode  II 
delamination  (like  that  for  mode  I  delamination)  appears  to 
be  ductile  rupture,  with  occasional  fiber  debonding.  Thus, 
it  is  not  surprising  that  the  gIIc/gjc  ratio  for  delamination 
of  composites  made  with  ductile  resins  is  approximately  1. 

Mixed-Mode  Delamination--A  series  of  mixed  mode 
delaminaticn  experiments  were  run  in  which  the  percentage  of 
mode  II  energy  release  rate,  Gllc/  was  systematically  varied 
from  0*  to  100%.  These  tests  were  performed  on  several  resin 
systems.  where  the  brittle  AS4/3502  was  used,  the  density  of 
the  hackles  as  well  as  their  orientation  relative  to  the 
delamination  planes  were  found  to  increase  with  increasing 
percentage  of  mode  II  loading,  as  seen  in  Fig.  11.  The  total 
critical  energy  release  rate  was  also  found  to  mono ton i ca 1 ly 
increase  with  increasing  percentage  of  mode  II  energy  release 
rate,  as  seen  in  Fig.  12  for  composites  made  with  brittle 
resins.  No  such  systematic  variation  in  either  fracture 
surface  appearance  or  critical  energy  release  rate  was  noted 


4-  in  situ  and  post  mortem  fraetography  of  AS4/ 3502  under  mode  II 

del  ami nation .  (a)  3000x.  Microcrack  formation  at  45  degrees  followed 
by  opening  of  microcracKs.  (b)  3000x .  Microcrack  coalescence. 

(c«  30oijx .  Microorack  coalescence  with  rotation  at  root  of  hackles. 


9.  (Continued)  (d)  2000x.  (e)  200 Ox.  Final  coalescence  of  microcracks 
results  in  macroscopic  crack  formation  showing  the  typical  hackles 
characteristic  of  mode  II  delamination  of  brittle  composites. 


igure  9.  (Continued).  Post-mortem  fractography  showing  the  highly  hackled 
fracture  surface,  (f)  lOOOx.  (g)  lOOOx.  (h)  2000x. 
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Resin  rich 
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(a)  Microcrackinq  ahead  of  macroscopic  crack 


Resin  rich 
reqion  between 
plies 


(b)  Microcrack  coalescence  (microcracks  beqin 
to  open  up  first) 


Resin  rich 
reqion  between 
pi  i  es 


(c)  Final  microcrack  coalescence  is  accompanied 

with  hackle  rotation  at  the  root  of  the  hackles. 


Fiqure  10.  Schematic  of  mode  II  delamination  crack  formation 
(hackle  formation  process). 


for  the  ductile  T6T145/F185  as  the  percentage  of  mode  II 
loading  was  increased. 

While  most  of  the  laminates  studied  in  this  work  were 
unidirectional  layups,  a  few  multi-directional  layups  were 
also  investigated.  Typical  results  for  mode  I  delamination 
of  a  specimen  containing  +45°  plies  across  the  delamination 
plane  are  seen  in  Fig.  13. 

4.2  Strain  Field  Mapping 

In  the  previous  section,  the  physical  details  of  the 
mode  I  and  mode  II  delamination  fracture  process  as 
determined  by  direct,  real-time  observations  in  the  SEM  were 
summarized.  It  was  decided  early  in  this  program  that 
quantifying  the  strain  field  around  the  crack  tip  would  be 
essential  to  assist  in  the  development  and  verification  of 
micromechanics  models.  Two  different  approaches  were  taken 
to  try  to  achieve  this  objective:  (1)  the  stereo-imaging 
approach  developed  for  metals  by  Davidson  and  Lankford  at 
Southwest  Research  Institute  and  (2)  the  direct  measurement 
of  the  distortion  of  a  dot  pattern  placed  on  the  surface  with 
the  electron  beam  of  the  SEM.  These  two  approaches  will  be 
explained  and  results  from  each  presented  next. 

4.2.1  Strain  Field  Mapping  by  Stereo-Imaging 

Recent  works  by  Davidson  and  Lankford3-'’  have 
described  a  method  of  accurately  determining  the  in-plane 
displacements  and  strain  fields  developed  during  the 
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Figure  12.  Total  G  vs.  %  mode  II  showing  increase  in  brittle  systems 

( AS1!/ 3501 -6 )  and  little  change  in  ductiie  systems  (T6T145/F185 


deformation  of  a  specimen.  The  technique  involves 
photographing  the  crack  tip  region  of  the  specimen  both 
before  loading  and  while  loaded,  and  then  using 
stereo-imaging  to  determine  the  displacement  field.  The 
in-plane  strains  are  then  calculated  by  differentiating  this 
displacment  field.  The  advantage  of  this  technique  is  that 
high-resolution  strain  field  maps  can  be  obtained  for  any 
level  of  loading,  including  loading  that  produces  nonlinear 
deformation  around  the  crack  tip. 

Previously  this  experimental  approach  has  only  been 
utilized  to  determine  the  strain  field  around  the  crack  tip 
of  metals  subjected  to  fatigue  stressing3  ^ .  In  this  study, 
similar  stereo-imaging  techniques  have  been  developed  to 
examine  the  strain  fields  associated  with  delamination  in 
graphite/epoxy  composites.  Stereo-photograph  pairs  (one  of 
the  loaded  specimen,  a  second  of  the  unloaded  specimen)  of  a 
delamination  crack  tip  for  both  mode  I  and  mode  II  loading 
conditions  are  made  in  a  scanning  electron  microscope  (SEM). 
From  the  displacement  and  strain  field  maps  generated,  the 
effects  of  macroscopic  loading  as  well  as  the  role  of  fiber 
constraint  on  the  deformation  and  resin  yielding  in  the 
region  of  the  crack  tip  were  examined.  All  data  reduction 
routines  used  in  the  strain  field  analysis  were  performed  on 
an  Apple  Computer. 

4. 2. 1.1  Stereo-Imaging 

We  perceive  three  dimensions,  or  depth  of  field, 
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because  our  eyes  view  our  surroundings  from  slightly 
different  angles.  The  apparent  position  of  foreground 
objectives  relative  to  the  distant  background  seen  from  one 
eye  is  slightly  displaced  from  that  seen  through  the  other 
eye.  As  the  brain  melds  these  two  images  into  one,  it  takes 
note  of  the  relative  displacements  of  the  foreground  objects 
and  interprets  it  as  depth. 

In  determining  the  in-plane  delamination 
displacements,  an  area  well  ahead  of  the  crack  tip  in  an 
unloaded  specimen  is  photographed,  ensuring  an  unloaded 
reference  condition.  The  specimen  is  then  loaded  and  the 
crack  allowed  to  grow  into  the  field  of  view,  which  is 
subsequently  rephotographed.  In  the  second  case,  due  to  the 
loading,  localized  deformation  in  the  crack  tip  region  will 
appear  as  displacements  relative  to  the  surrounding 
unstrained  regions,  as  captured  in  the  photograph  of  the 
specimen  prior  to  loading.  When  the  two  photographs  are 
observed  in  a  stereoscope,  these  relative  surface 
displacements  are  perceived  are  as  relief  in  the  out  of  plane 
or  2-direction.  It  is  noted  that  a  state  of  triaxial 
displacement  may  exist  at  the  loaded  crack  tip,  but  only  the 
in-plane  displacement  can  be  observed  using  this  technique. 

4. 2. 1.2  Experimental  and  Analytical  Procedures 

The  graphite-epoxy  systems  examined  in  this  study  were 
Hercules  AS4/3501-6  and  AS4/D0W  P-7.  These  resin  systems 
were  selected  because  of  their  differences  in  ductility  and 


A, 
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interfacial  adhesion  strengths.  The  Hercules  AS4/3501-6 

2 

system  consisted  of  a  brittle  resin  (Gj  =  70  J/m  )  along 

with  unsized  fibers,  resulting  in  a  composite  delamination 

2  2 

toughness  of  Gj  -  140  J/m  and  Gj  =  634  J/m  .  The  AS4/Dow 
P7  system  used  sized  fibers  with  a  tougher  Novolac  resin, 
modified  to  increase  both  toughness  and  resin/fiber  adhesion, 

resulting  in  a  composite  delamination  toughnesses  of  G  - 

2  2 
340  J/m  and  Gjj  -  862  J/m  . 

Small  (25  mm  x  11  cm),  undi rectional  double  cantilever 

beam  specimens,  14  plies  in  thickness,  were  cut  from  cured 

laminate  panels  of  the  two  systems.  A  thin  0.03  mm  teflon 

sheet  inserted  at  one  edge  during  the  layup  of  the  laminate 

panels  before  curing  provided  a  starter  crack  between  the 

mid-plies.  Delamination  loading  of  the  specimens  within  the 

SEM  was  accomplished  by  using  a  tensile  stage.  Mode  I 

loading  was  achieved  by  driving  the  precracked  portion  of  the 

specimen  onto  a  fixed  wedge  (see  Fig.  14).  The  wedge  was 

sufficiently  blunt  to  ensure  that  it  remained  well  behind  the 

crack  tip.  A  three  point  bend  fixture  was  used  to  provide 

pure  shear  or  mode  II  loading  at  the  crack  tip  (see  Fig.  14). 

For  neat  resin  testing,  compact  tension  specimens  of 

size  IT  were  used  as  per  ASTM  E399^,  except  that  the 

thickness  was  3-4  mm  instead  of  the  12.5  mm  prescribed  by  the 

standard.  The  same  loading  fixtures  were  originally  used  for 

the  in-situ  fracture  observations  reported  in  section  4.1  of 

this  report. 

As  described  before,  an  area  well  ahead  of  the  crack 
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tip  was  photographed,  and  then  rephotographed  after  the  crack 
was  allowed  to  grow  into  the  field  of  view,  representing  the 
unloaded  and  loaded  conditions  respectively.  These 
photographs  were  than  viewed  with  a  TOPCON  stereoscope.  A 
1/4-inch  reference  grid  is  marked  on  one  photograph  to  enable 
the  location  on  the  specimen  to  be  identified.  A  parallax 
bar  is  then  used  to  measure  the  apparent  height  ( 2x )  of  the 
surface  at  each  of  the  points  of  the  reference  grid.  A 
reference  point  far  ahead  of  the  crack  tip  outside  of  the 
deformation  zone  is  selected  to  provide  a  reference  point 
( RP )  of  zero  displacement  (see  Fig.  2).  Displacement  ( Dx )  is 
calculated  at  each  point  on  the  reference  grid  using  a 
standard  stereo  imaging  analysis  by  the  following 
relationship^ : 

Dx  -  ( 1/M)  (dP)/( P  +  dP)  ( 1  ) 

where  M  is  the  photographic  magnification,  dP  is  the 
difference  in  the  parallax  bar  measurement  at  a  given  point 
on  the  reference  grid  minus  the  parallax  bar  measurement  at 
the  reference  point  (Zx-RP),  and  P  is  the  average  separation 
between  the  conjugate  center  and  conjugate  principal  points 
for  the  two  pictures  used  in  the  stereo  pair.  A  more 
complete  discussion  by  Avery^  of  this  technique  is  available. 

Since  only  displacements  along  the  eye  axis  can  be 
observed  in  the  stereoscope,  the  photographs  must  be  rotated 
90  degrees  and  the  procedure  repeated  to  obtained  the 
Y-direction  displacements  ( Dy ) .  Additional  details  may  be 
found  in  a  soon  to  be  published  dissertation  by  Hibbs. 
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(a)  Schematic  of  mode  I  delamination  of  scanning  electron  of 
microscopy  specimens,  (b)  Schematic  of  mode  II  delamination 
of  scanning  electron  microscopy  specimens,  (c)  Schematic  of  mode 
I  fracture  of  scanning  electron  microscopy  CT  specimens. 
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gore  Ad. 


Relative  displacement  of  points  near  the  cracK  tip 
resulting  from  the  fading  of  the  specimen.  When  view  with  a 
stereoscope  these  relative  displacements  are  seen  as 
vertical  relief.  The  point  xQ  is  selected  far  ahead  of  the 
deformation  zone  and  is  used  as  a  reference  point  of  zero 
d i sp .^cement . 
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4. 2. 1.3  Experimental  Results 


Results  for  mode  I  loading  of  the  AS4/'Dow  P7  system 
and  mode  II  loading  of  the  AS4/3501-6  system  will  be 
presented  in  this  report  with  more  extensive  results  to  be 
presented  in  Michael  Hibbs'  dissertation. 

Mode  I  Displacement  and  Strain  Fields--The 
unloaded/loaded  stereo  pair  for  the  AS4/Dow  P7  composite 
system  loaded  in  mode  I  conditions  is  given  in  Fig.  15.  To 
show  the  displacement  field  more  clearly,  the  displacement 
vectors  as  determined  from  the  photographs  in  Fig.  15  using 
the  TOPCON  stereoscope  are  represented  graphically  in 
Fig.  16.  The  displacement  field  is  not  symmetric  about  the 
crack  tip.  This  probably  results  from  asymmetric  constraints 
due  to  the  asymmetric  distribution  of  fibers  on  either  side 
of  the  plane  of  the  crack.  The  maximum  displacements  are 
seen  to  be  approximately  0.25  microns.  At  the  crack  tip  the 
loading  is  shown  to  be  essentially  pure  mode  I. 

The  resulting  strain  field  seen  in  Fig.  17  shows  that 
the  maximum  shear  strain  is  for  the  most  part  confined  to  the 
resin  rich  area  between  the  fibers,  with  a  maximum  of  about 
20%  strain  at  the  crack  tip.  A  more  ductile  resin  would  be 
expected  to  have  deformation  outside  the  resin  rich  region 
between  plies.  with  increasing  distance  ahead  of  the  crack 
tip,  the  magnitude  of  the  strain  decreases,  reaching  a  value 
of  7%  strain  at  a  distance  of  33  microns.  It  is  interesting 
to  note  that  the  tensile  elongation  measured  for  this 
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Figure  <5.  The  unloaded  loaded  stereo  photographs  of  the  composite 
AS^  Dow  P7  loaded  in  mode  I  conditions.  (1Q00X) 
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relatively  brittle  resin  is  less  than  2%.  We  should  hasten 
to  add  that  the  indicated  crack  tip  shear  strain  of  20%  is 
not  true  resin  strain  because  of  the  fine  microcracking  that 
develops  in  the  crack  tip  region.  The  true  resin  strain 
would  be  less.  Nevertheless,  displacements  due  to  resin 
deformation  and  microcracking  combine  to  give  deformation 
shear  equivalent  to  20%  strain. 

The  displacement  field  that  results  for  mode  II 
loading  of  AS4/3501-6  is  given  in  Fig.  18,  where  a  symmetric 
displacement  field  with  essentially  pure  mode  II  loading  at 
the  crack  tip  is  seen.  In  this  field,  the  magnitude  of  the 
displacements  are  on  the  order  of  0.75  microns. 

The  resulting  strain  field  (Fig.  19),  as  in  the  case 
for  the  mode  I  loading  of  the  Dow  P7  system,  is  seen  to  be 
confined  in  the  resin  rich  region  between  the  fibers,  again  a 
result  that  would  not  be  true  for  a  composite  made  with  a 
more  ductile  resin.  A  maximum  shear  strain  of  approximately 
100%  occurs  at  the  crack  tip,  with  a  very  slow  reduction  in 
magnitude  with  increasing  distance  ahead  of  the  crack  tip. 
At  the  edge  of  the  observed  field,  33.0  microns  ahead  of  the 
crack  tip,  the  magnitude  of  the  shear  strain  was  still  85%. 

It  is  difficult  to  see  how  a  composite  with  such  a 
brittle  resin  as  the  AS4/3501-6  could  have  such  a  large  local 
strain  to  failure  until  one  recognizes  that  most  of  the 
measured  displacements  are  (as  was  the  case  for  mode  I)  due 
to  extensive  microcracking,  hackle  formation,  and  hackle 
rotation.  These  can  be  seen  in  Figs.  6,9  and  10  in  section 
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4.1. 

■ 

These  results  for  the  strain  field  for  mode  I  and  mode 

II  delaminaticn  are  consistent  with  the  observations  of 

deformation/damage  zone  size  previously  reported  in  4.1  and 

will  be  seen  later  in  this  report  (section  4.3)  to  be 

|  qualitatively  consistent  with  finite  element  calculations  of 

f  the  strain  field  decay  ahead  of  a  crack  for  these  two  loading 

> 

>  conditions. 

i 

4.2.2  Direct  Measurement  of  Strain  Field  Using  a  Dot  Pattern 

The  purpose  of  this  study  has  been  to  develop  a 
technique  which  would  permit  the  direct  measurement  of  the 
strain  field  around  a  crack  tip  with  sufficient  resolution  to 
make  a  comparison  of  the  fields  for  cracks  in  neat  resin 
,  versus  delamination  cracks  in  composites  made  from  the  same 

resin.  The  approach  used  has  been  to  burn  a  very  fine  array 
of  dots  onto  the  polished  surface  (perpendicular  to  the  plane 
of  the  crack)  of  a  specimen  using  the  electron  beam  in  a 
scanning  electron  microscope.  The  exact  coordinates  for  each 
of  these  points  are  determined  before  and  after  loading  the 
specimen,  allowing  the  determination  of  the  displacement  of 
each  point.  This  displacement  field  may  then  be 

differentiated  to  give  the  strain  field  around  the  crack  tip. 

4. 2. 2.1  Experimental  Procedure 

With  considerable  trial  and  error,  a  technique  has 
been  developed  to  burn  a  very  systematic  dot  pattern  onto  the 
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polished  surface  of  the  specimens.  The  beam  of  the  electron 
microscope  was  adjusted  to  give  dots  that  were  sufficiently 
large  to  remain  easily  recognizable  after  considerable 
deformation  but  small  enough  to  give  adequate  resolution  and 
minimal  damage  to  the  specimen.  The  very  regular  dot  pattern 
seen  in  Fig.  20  was  burned  onto  the  surface  using  a  JEOL-35 
scanning  electron  microscope  with  the  assistance  of  an  image 
processing  system  manufactured  by  Kontron.  The  dot  spacing 
seen  in  Fig.  20  is  approximately  10  microns.  A  finer  dot 
spacing  is  possible,  but  not  very  useful  unless  the  size  of 
the  dots  can  also  be  reduced,  while  retaining  their 
visibility  after  large  deformation. 

The  resin  system  selected  for  use  in  this  study  was 
Hexcel  FI 8 5 .  This  resin  was  selected  because  it  is  known  to 
be  relatively  ductile  (8%  elongation  in  a  tensile  test),  and 
therefore  would  be  expected  to  give  significant  strain  in  the 
vicinity  of  the  crack  tip  prior  to  failure. 

Hexcel  F185  resin  was  cast  into  rectangular  plates. 
Compact  tension  specimens  were  machined  from  these  plates  and 
subsequently  fatigue  precracked.  Then  the  CT  specimens  were 
polished  on  a  microprocessor-controlled  automatic  polisher  to 
a  0.03  micron  finish.  Finally,  a  dot  pattern  was  burned  on 
the  surface  using  the  electron  beam  in  the  SEM,  as  shown  in 
Fig.  20.  The  crack  tip  is  seen  at  the  bottom  of  Fig.  20, 
just  below  the  dot  pattern.  After  loading  to  grow  the  crack 
a  small  amount  into  the  dot  array,  the  deformation  field  was 
determined  using  the  distortion  of  the  dot  pattern.  The 
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displacement  field  and  the  associated  strain  field  were 
determined  using  the  image  processing  system  in  conjunction 
with  a  computer  program  written  to  numerically  differentiate 
the  displacement  field  data. 

4. 2. 2. 2  Results  and  Discussion 

The  results  are  presented  in  Figs.  21-23.  The  normal 
strains  as  well  as  the  principal  shear  strain  are  seen  as  a 
function  of  x-y  coordinates.  The  crack  tip  has  coordinates 
of  y-57.48,  x»59.14.  The  compact  tension  specimen  was  loaded 
in  the  x-direction  with  crack  growth  in  the  x-direction. 

It  should  be  emphasized  that  the  three  dimensional 
plotting  routine  used  to  present  the  results  graphically  does 
local  smoothing  prior  to  plotting.  The  polynominal  curve  fit 


for 

the 

smoothing 

operation  interpolates 

as  well,  giving  a 

much 

finer 

grid 

in  the  graphical  results 

than  was  actually 

used 

on 

the 

specimen  surface.  (Contrast 

actual  grid  size, 

Fig. 

20, 

to 

grid 

size  seen  in  the  three 

dimensional  plots, 

Figs  . 

21- 

23.  ) 

The  strain  at  yielding  in  F185  as  measured  in  a 


tensile  test  is  about  1.2%.  Thus,  the  region  of  nonlinear 
deformation  around  the  crack  tip  is  seen  to  be  quite 
extensive.  If  it  is  remembered  that  carbon  fibers  typically 
have  diameters  of  about  7  microns,  then  the  nonlinear 
deformation  zone  would  extend  to  at  least  five  fiber 
diameters  above  and  below  the  plane  of  delamination  if  the 
composite  has  a  similar  strain  field  around  the  tip  of  a 
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growing  delamination  crack  to  that  observed  in  Figs.  21-23 

for  the  F 1 8 5  neat  resin.  A  qualitative  indication  of  a 

deformation  zone  of  this  size  in  a  delaminating  T6T145/F185 

composite,  is  seen  in  Fig.  7  (section  4.1)  and  has  previously 

2 

been  noted  by  Bradley  and  Cohen  . 

A  second  important  observation  concerns  the  magnitude 
of  the  local  strain  at  the  crack  tip.  The  elongation 
measured  in  a  tensile  test  of  F185  is  8%,  with  no  apparent 
necking.  Thus,  one  might  have  assumed  for  micromechanics 
modelling  that  the  local  strain  to  failure  would  be  8%,  or 
possibly  less  due  to  the  effects  of  constaint.  However,  it 
is  clear  that  the  tensile  elongation  greatly  underestimates 
the  true  local  strain  to  failure  at  the  tip  of  a  growing 
crack  in  F185.  The  measured  value  of  48%  should  be 
considered  a  lower  bound  since  it  is  measured  over  a  10 
micron  gage  length  at  the  crack  tip  in  a  very  steep  strain 
gradient  and  would  increase  to  an  even  higher  value  as  crack 
advance  caused  ductile  fracture.  The  critical  value  for 
local  strain  to  failure  for  delamination  crack  growth  is 
probably  somewhat  lower  than  the  value  for  crack  growth  in 
the  neat  resin  measured  here  due  to  debonding  and  constraint 
produced  by  the  fibers.  Comparable  results  for  crack  growth 
in  a  delaminating  T6145/F185  composite  were  not  available  at 
the  time  of  the  preparation  of  this  report  but  should  be 
available  soon  in  the  new  grant  period. 

One  concern  that  should  be  mentioned  is  the  question 
of  whether  the  surface  strain  measurements  are  in  any  way 
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representative  of  the  bulk  strain  distribution.  In  materials 
where  cavitation  or  mi c r oc r a ck i ng  does  not  occur,  one  would 
expect  a  significant  difference  between  the  strain 
distribution  for  the  plane-stress  state  of  stress  that  exists 
at  the  surface  and  the  strain  distribution  that  would  exist 
for  the  more  nearly  plane-strain  state  of  stress  in  the  bulk. 
However,  F185  with  rubber  particle  toughening  does  cavitate, 
relaxing  the  through-thickness  stress.  Thus,  for  this  system 
we  believe  the  surface  strain  distribution  is  representative 
of  the  bulk  strain  distribution. 

Direct  comparison  of  results  from  stereo-imaging 
with  results  based  on  dot  pattern  distortion  will  also  be 
available  soon  under  the  new  contract. 

4.3  Analysis  and  Modeling 

Analysis  and  modeling  have  been  used  in  two  ways  to 
better  understand  the  experimental  results  of  sections  4.1, 
4.2,  and  4.4.  First,  the  very  different  size  and  shape  of 
the  damage  zone  for  mode  I  loading  compared  to  mode  II 
loading  has  been  studied  using  a  finite  element  analysis  in 
which  the  constitutive  relationship  for  a  linear,  orthotropic 
continuum  is  assumed.  While  ignoring  the  details  that  would 
be  treated  explicitly  in  a  micromechanics  model,  this 
approach  has  been  quite  successful  in  accounting  for  the 
major  distinctives  between  mode  I  and  mode  II  delami na t i on . 
Second,  a  simple  model  has  been  proposed  to  account  for  the 
observed  differences  in  neat  resin  toughness  compared  to  mode 
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I  delamination  toughness  for  composites  made  with  ductile 
resins.  The  results  of  these  two  areas  of  investigation  will 
be  presented  in  this  section. 

4.3.1  Linear,  Orthotropic  Finite  Element  Analysis  of  Mode  I 

and  Mode  II  Delamination 

In  this  investigation,  the  stress  fields  ahead  of  the 
crack  tip  of  a  split  laminate  specimen  loaded  to  give  mode  I 
delamination  and  mode  II  delamination  have  been  determined  by 
means  of  a  finite  element  analysis.  The  results  are  compared 
to  the  observed  damage  zones  developed  ahead  of  the  crack  tip 
for  the  same  loading  conditions,  as  reported  in  section  4.1. 
A  comparison  with  the  direct  strain  field  measurements  of 
section  4.2  is  made. 

4. 3. 1.1  Materials  and  Methods 

The  determination  of  the  stress  field  ahead  of  the 
crack  tip  of  the  split  laminate  specimen  loaded  under  mode  I 
and  mode  II  was  accomplished  by  first  generating  a  two 
dimensional  mesh  152.4  mm  long  by  2.54  mm  thick  consisting  of 
791  nodes  and  294  elements  as  shown  in  Fig.  24a.  Triangular 
six-noded  elements  were  used  around  the  crack  tip  with 
mid-side  nodes  displaced  to  the  quarter  point7,  and  a 
substantial  refinement  of  the  mesh  around  the  crack  tip  was 
made  to  overcome  the  difficulty  imposed  by  the  stress 
singularity  present  at  the  crack  tip  .  This  linear 
orthotropic  analysis  has  been  made  with  a  finite  element 
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(b)  Mode  I  boundary  conditions. 
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(c)  Mode  II  boundary  conditions. 


Figure  24.  Finite  element  mesh  and  boundary  conditions 


g 

algorithm  developed  by  Henriksen  .  The  algorithm  is  based  on 

a  nonlinear  code  updated  with  Lagrangian  formulation  using 

six  and/or  eight  noded  isoparametric  elements  with  two 

degrees  of  freedom  per  node.  Elastic  constants  for  a 

unidirectional  laminate  of  AS4/3502  graphite/epoxy  composite 

9 

to  be  used  in  the  analysis  were  obtained  from  Hercules  . 

Mode  I  loading  was  simulated  by  applying  a  symmetric 
load  at  the  cracked  end  of  the  mesh  as  shown  in  Fig.  24b. 
Mode  II  loading  was  introduced  by  asymmetrically  loading  the 
cracked  end  (see  Fig.  24c).  The  load  level  applied 

corresponds  to  approximately  the  load  at  onset  of  crack 
growth  for  these  two  modes  of  failure,  as  determined  from 
experimental  measurements  in  combination  with  beam  linear 
beam  theory. ^  Stresses  were  obtained  and  were  plotted  as  a 
function  of  distance  ahead  of  the  crack  tip. 

4. 3. 1.2  Results  and  Discussion 

Figures  25  and  26  show  the  finite  element  results  for 
both  mode  I  and  mode  II  loading.  Figure  25a  is  a  Syy 
(normal)  stress  contour  plot  from  the  vicinity  of  the  crack 
tip  for  mode  I  loading.  The  normal  stress  Syy  which  for  this 
loading  and  geometry  is  the  principal  normal  stress  S^, 
rapidly  decreases  ahead  of  the  crack  tip.  Figure  25b 
corresponds  to  the  Sxy  (shear)  stress  contour  plot  for  a  mode 
II  condition.  Note  how  the  shear  stress  drops  off  much  more 
slowly  with  distance  ahead  of  the  crack  tip  than  does  the 
normal  stress  for  mode  I  loading  (Fig.  25a).  Furthermore, 
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(a)  Syy  *  S2,  mode  I. 


(  b)  Sxy  s  S , ,  mode  1 1 . 


Figure  25.  Stress  contour  plots  for  mode  I  and  mode  II 
loadi ng 
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the  shape  of  the  stress  field  is  more  narrow  and  elongated. 
For  this  loading  condition,  Sxy  =  maximum  shear  stress  =  S^. 
Figure  26  shows  the  stress  distribution  as  a  function  of 
distance  ahead  of  the  crack  tip  for  mode  I  and  mode  II.  In 
the  case  of  the  normal  stress  for  mode  I,  the  stress  drops 
off  rapidly  until  it  is  compressive  at  a  distance  of  0.76  mm 
ahead  of  the  crack  tip.  Finally,  it  gradually  approaches  a 
zero  stress  level  at  approximately  8  mm  from  the  crack  tip 
which  is  maintained  all  the  way  to  the  end  of  the  beam.  It 
should  be  noted  that  the  compressive  stresses  observed  for 
mode  I  loading  are  not  expected  to  significantly  influence 
the  fracture  mode  at  the  crack  tip  because  they  occur  far 
enough  away  from  the  crack  tip.  The  mode  I  damage  zones 
sizes  presented  in  section  4.1  (see  Table  1)  were  between 
.020  mm  and  .200  mm,  well  below  the  .76  mm  indicated  in  the 
analysis  as  the  distance  ahead  of  the  crack  beyond  which  the 
stresses  are  compressive. 

In  the  case  of  the  shear  stress  distribution  ahead  of 
the  crack  tip  for  mode  II,  it  monotonically  decreases  to  a 
constant  value.  As  can  be  seen,  the  shear  stress  ahead  of 
the  crack  tip  for  mode  II  loading  decays  much  more  slowly 
than  the  normal  stress  for  mode  I  loading.  These  results 
indicate  that  the  stress  concentration  at,  the  crack  tip  is 
distributed  over  a  larger  distance  for  mode  II  loading  than 
for  mode  I  loading.  For  mode  II,  the  extent  of  the  stress 
field  above  and  below  the  center  plane  is  much  smaller  than 
for  mode  I,  as  seen  in  the  stress  contour  plots  in  Fig.  25. 
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Figure  27.  Von-Mises  stress  distribution  ahead  of  the  crack 
tip  for  mode  I  and  mode  II  loading 
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A  careful  comparison  of  the  stress  fields  of  Figs.  25 
and  26  with  the  damage  zone  shown  in  Figs.  6-8  for  mode  I  and 
mode  II  delamination  reveals  a  direct  correspondence  between 
them.  The  much  longer  and  more  narrow  stress  field  predicted 
by  the  linear,  orthotropic  analysis  is  clearly  reflected  in 
the  observed  damage  zones  for  mode  I  and  mode  II 
delamination.  There  is  also  quantitative  agreement  with 
measured  strain  fields,  which  also  indicate  a  much  more 
gradual  decay  of  the  strain  field  ahead  of  the  crack  tip  for 
mode  II  loading  than  for  mode  I.  It  should  be  emphasized 
that  quantitative  agreement  between  the  directly  measured 
strain  field  and  the  calculated  stress  field  would  not  be 
expected.  The  stress  field  calculations  ignored  both 
nonlinear  behavior  and  microstructural  details  (fibers  and 
resin)  since  they  assumed  linear,  orthotropic  behavior. 
Furthermore,  there  is  no  easy  way  to  get  from  the  measured 
stcain  to  the  associated  stress,  since  the  nonlinear 
constitutive  relationship  is  unknown. 

When  dealing  with  damage  zone  development,  the  effect 
of  the  state  of  stress,  namely  plane  stress  versus  plane 
strain,  may  be  important.  This  is  particularly  true  where 
the  principle  shear  stress  plays  a  role  in  .the  deformation  of 
fracture  process.  The  damage  zones  shown  in  Figs.  6  and  7 
were  obviously  made  at  the  edge  of  the  specimen  where  a  plane 
stress  condition  occurs.  Therefore,  in  order  to  have  a  first 
order  approximation  as  to  the  effect  the  difference  in  the 
state  of  stress  at  the  surface  of  the  specimen  where  our 
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direct  observations  are  made,  and  at  the  center  were  plane 
strain  conditions  prevail,  the  Von-Mises  stress  for  mode  I 
and  mode  II  were  plotted  as  functions  of  distance  ahead  of 
the  crack  tip.  The  results  in  Fig.  27a  indicate  that  for 
mode  I,  the  difference  in  the  Von-Mises  stress  (or  root  mean 
square  shear  stress)  in  plane  stress  versus  plane  strain  is 
approximately  10%.  In  the  case  of  mode  II,  the  difference  is 
negligible,  as  is  evident  in  Fig.  27b.  Therefore,  it  appears 
that  for  these  orthotropic  materials,  the  difference  between 
plane  stress  and  plane  strain  is  not  too  significant  and  our 
surface  observations  of  damage  zone  size  and  detail  should 
not  be  very  different  than  subsurface  behavior.  Previous 
studies  in  which  the  postmortem  appearance  of  the  fracture 
surface  near  a  free  edge  to  that  at  the  center  of  a  specimen 
are  consistent  with  this  interpretation^.  A  final 
observation  from  Fig.  27  is  that  the  stress  field  ahead  of  a 
crack  tip  decays  much  more  slowly  for  an  orthotropic 
materials  with  the  crack  parallel  to  the  fibers  than  for  an 
isotropic  material. 

4.3.2  Modeling  the  Effect  of  Fibers  on  Mode  I  Delamination 

Toughness 

A  first  order  estimate  of  the  delamination  fracture 
toughness  may  be  obtained  by  assuming  that  the  delamination 
fracture  process  is  essentially  the  same  as  the  fracture 
process  in  the  neat  resin,  except  that  the  fibers  act  as  a 
rigid  filler,  reducing  the  volume  of  resin  available  to 
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deform  in  the  crack  tip  region.  It  has  previously  been 

2 

suggested  by  Bradley  and  Cohen  that  the  energy  dissipation 
per  unit  area  of  crack  extension  may  be  calculated  by 
picturing  a  small  tensile  specimen  ahead  of  the  crack  tip, 
which  is  slowly  stretched  as  the  crack  tip  approaches  and 
finally  breaks  as  the  crack  tip  passes  (Fig.  28).  The  energy 
absorbed  per  unit  area  of  crack  extension  would  be  calculated 
for  such  a  model  by  summing  the  energy  absorbed  per  unit 
volume  of  material  over  the  volume  of  the  hypothetical 
tensile  specimen.  Such  a  summation  may  be  written 
mathematically  as  follows: 

Gic  -  ffndendy  (1) 

where  GIC  is  the  fracture  toughness  of  the  material,  2h  is 
the  height  of  the  hypothetical  tensile  bar,  and  is  the 
component  of  the  stress  tensor  normal  to  the  plane  of  the 
crack.  Until  recently,  neither  the  strain  field  nor  the 
constitutive  relationship  required  to  calculated  stress  from 
strain  was  known.  Thus,  a  direct  evaluation  of  Eq .  1  was  not 
possible.  A  phenomenological  approach  has  been  taken  to  try 
to  evaluate  this  energy  release  rate  for  delamination . 

The  addition  of  fibers  to  the  neat  resin  can  perturb 
the  energy  calculation  in  Eq.  1  in  several  ways.  First,  the 
extent  of  the  deformation  zone  might  be  changed  due  to  the 
additional  constraint  imposed  by  the  fibers,  changing  the 
value  of  h  (i.e.,  the  height  of  our  hypothetical  tensile 
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specimen)  in  Eq .  1.  Second,  rigid  filler  would  further 

reduce  the  effective  gage  length  of  material  capable  of 

deforming  in  our  hypothetical  tensile  specimen.  Third,  fiber 

debonding  and/or  constraint  could  change  the  local  strain  to 

fracture.  These  three  factors  could  potentially  account  for 

the  observed  decrease  in  fracture  toughness  from  8100  J/m*"  in 

2 

the  F185  resin  to  1900  J/m  in  the  composite  delamination 
fracture  toughness.  The  first  two  of  these  three  factors  can 
be  quantified  based  on  actual  observations.  If  one  assumes 
that  coating  microcracking  begins  for  strains  above  the 
threshold  strain;  (say  3s),  then  the  extent  of  the 

microcracked  zone  can  be  used  to  approximately  quantify  the 
magnitude  of  the  h  in  Eq .  2. 

A  cross-section  of  the  composite  prepared 

me ta 1 1 og raph i ca 1 1 y  to  reveal  the  microstructure  (Fig.  29)  may 
be  used  to  determine  the  volume  fraction  of  fibers  in  the 
hypothetical  tensile  specimen  ahead  of  the  crack  tip.  Since 
the  volume  fraction  is  quite  nonuniform,  the  microstructure 
can  be  divided  into  three  regions:  the  resin  rich  region 
between  plies  with  a  fiber  volume  fraction  of  19%;  the  ply 
region  with  a  volume  fraction  of  fibers  of  76%,  and  a 
transition  zone  with  a  volume  fraction  of  approximately  33%. 
The  relative  heights  of  these  three  regions  are  shown  in  the 
schematic  in  Fig.  30  along  with  the  height  of  the 

microcracked  zone  in  the  F185  resin.  Using  a  simple  rule  of 
mixtures  approach,  and  taking  account  of  the  smaller 
microcracked  (and  deformed)  zone,  2h,  one  may  estimate  the 
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delamination  fracture  toughness  to  be  4000  J/m^  as  shown 
below: 


GIC  =  (ha(l-VAf)  +  hfa  (l-VBf)  +  hc  (  l-VCf  )  ]  (2) 

(composite  hQ 
delamination ) 

where  GIC  refers  to  composite  delmination  and  GIC  refers  to 
the  neat  resin.  This  calculation  implicitly  assumes  that  the 
local  strain  to  fracture  in  the  F185  resin  and  the  composite 
are  the  same  and  that  the  stress  distribution  is  also 
similar,  at  least  on  average.  It  is  recognized  that  these 
are  crude  assumptions  at  best. 

The  calculated  value  of  Mode  I  delamination  fracture 

2 

toughness  of  4000  J/m  is  clearly  much  larger  than  the 

2 

observed  value  of  1900  J/m  .  This  discrepancy  is  due  to  the 
assumptions  in  the  calculation  that  the  strain  field 
distribution  as  well  as  the  local  strain  to  failure  are 
similar.  In  reality,  fracture  at  the  fiber/resin  and/or 
constraint  due  to  the  presence  of  the  fibers  should  cause  the 
strain  to  failure  in  the  composite  to  be  somewhat  less  than 
that  in  the  neat  resin.  Furthermore,  the  decay  of  the  strain 
field  ahead  of  the  crack  tip  is  different  for  an  orthotropic 
material  than  for  an  isotropic  one  (See  Fig.  27). 

Thus,  the  incorporation  of  the  recent  measurements  of 
the  strain  field  described  in  section  4.2  into  this  model 
should  bring  the  predictions  more  in  line  with  the  observed 
delamination  toughness.  This  is  currently  being  done  under 
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the  new  contract.  A  constitutive  relationship  will  have  to 
be  either  determined  or  assumed  before  this  approach  can  be 
fully  implemented.  With  a  constitutive  relationship  and  the 
measured  strain  fields,  it  would  be  possible  to  evaluate 
directly  Eq.  1  rather  than  using  the  more  indirect, 
phenomenological  approach  expressed  in  Eq.  2. 

4.4  FRACTURE  MECHANICS  CHARACTERIZATION  OF  MATERIALS 

In  this  section  our  work  to  develop  better  test 
methods  and  analysis  to  determine  the  mode  I  and  mode  II 
delamination  fracture  toughness  of  composite  materials  will 
be  described.  The  newer  resin  systems  are  significantly 
tougher  than  the  first  generation  resins  for  which  the 
current  test  methodology  was  developed.  For  these  early 
systems  which  were  quite  brittle,  linear  load-displacement 
curves  for  split  laminate  testing  is  the  norm  and  linear 
elastic  analysis  is  adequate  to  extract  meaningful  fracture 
toughness  parameters  from  the  experimental  data.  As  it  has 
become  necessary  to  characterize  some  of  the  newer,  tough 
resins  such  as  PEEK,  both  nonlinear  and  inelastic  behavior 
are  frequently  observed  in  the  load-displacement  records,  and 
thus,  the  usual  linear  beam  theory  is  not  adequate  to  obtain 
meaningful  fracture  mechanics  parameters  from  these  test 
resul ts . 

It  has  also  been  found  that  mode  II  testing  with  the 
generally  much  higher  critical  energy  release  rates  observed 
for  this  mode  of  fracture  is  also  much  more  likely  than  mode 
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I  resting  to  exhibit  nonlinear  behavior.  The  ueiamination 
fracture  toughness  characterization  of  multi-axial  laminates 
with  the  associated  reduction  in  stiffness  along  the  axis  of 
the  split  laminate  being  tested  has  also  been  found  to  often 
give  nonlinear  behavior  unless  very  thick  specimens  are 
utilized. 

Our  efforts  in  response  to  these  challenges  have  been 
to  develop  a  J-integral  approach  to  delamination  fracture 
characterization  for  both  mode  I  and  mode  II  loading.  These 
two  efforts  will  be  described  in  the  two  subsections  that 
follow.  The  presentation  in  these  subsections  will  be 
somewhat  more  detailed  than  those  in  4.1,  4.2  and  4.3  as  most 
of  the  results  have  not  yet  been  published. 

4.4.1  Mode  I  Delamination 

The  subject  of  the  work  reported  in  this  section  is 
the  means  of  determining  the  mode  I  delamination  fracture 
toughness  of  multidirectional  composites.  Material 
characterization  is  a  prerequisite  for  systematic  design  and 
analysis  of  composite  structures.  The  intent  is  to  develop 
and  verify  a  theoretically-based  method  for  characterizing 
delamination  fracture  despite  the  complicating  factors  which 
will  be  discussed. 

4. 4. 1.1  Current  Status  of  the  Problem 

Delamination  of  Unidirectional  Composites — The 
delamination  of  unidirectional  composites,  though  not 
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directly  the  subject  of  the  present  study,  neveiLlieless  fornix 

the  background  for  it.  Therefore,  this  discussion  of  the 

current  status  of  the  problem  begins  with  its  review. 

The  great  majority  of  studies  of  delamination 

toughness  have  involved  unidirectional  composites.  A  major 

motivation  for  this  activity  has  been  to  understand  the 

relationship  between  neat  resin  fracture  toughness  and 

1 2 

composite  interlaminar  toughness  .  Several  different 
specimen  designs  have  been  used  by  various  investigators 
Of  special  interest  because  of  its  wide  use  and 
apparent  utility  is  the  double  cantilevered  beam  ( DCB ) 
specimen,  as  shown  in  Fig.  31.  Numerous  investigators  have 

1  £1  O 

studied  it  as  applied  to  unidirectional  composites  .  The 

DCB  specimen  has  also  been  applied  to  a  planar  isotropic 

1 9 

short-fiber  reinforced  composite  .  In  these  studies  the 
necessary  analyses  have  been  developed  and  the  procedure's 
viability  has  been  demonstrated  by  experiment. 

Several  different  ways  of  analyzing  DCB  data  have  been 
proposed.  Various  assumptions  about  material  and  specimen 
behavior  are  made  in  each  method.  Analyses  based  on  linear 
beam  theory,  large  displacement  theory,  empirical  compliance 
calibration,  the  area  under  the  load-displacement  curve,  and 
an  energy  rate  interpretation  of  the  J  integral  have  been 
used^ '  ^ .  The  effectiveness  of  a  particular  method 
depends  on  the  situation  at  hand.  For  the  case  where  there 
is  permanent  deformation  due  to  material  damage,  Keary  et 
al.^  determined  that  the  J  integral  is  superior  to 
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approaches  based  on  linear  beam  theory  and  empirical 
compliance  calibration.  This  observation  is  salient  to  this 
study  since  such  damage  has  been  observed  in  multidirectional 
composi tes . 

A  similar  conclusion  can  be  drawn  from  Fig.  32  which 
shows  data  generated  in  this  program  using  DCB  specimens. 
The  delamination  fracture  toughness  determined  using  the  area 
method  (G)  and  using  a  damage-dependent  J  integral  (J)  are 
plotted  versus  crack  length.  One  G  and  one  J  curve  were 
determined  from  each  of  three  specimens.  Not  only  does  J 
appear  to  give  results  which  are  more  consistent  over  the 
range  of  crack  lengths,  but  it  gives  lower  toughness  values. 
The  overestimation  of  fracture  energy  using  the  area  method 
reflects  that  methods'  unsuitability  when  damage  is  present 
away  from  the  crack  tip.  The  stacking  sequence  used  results 
in  damage  in  the  legs  of  the  specimens;  the  material  was 
Hexcel's  T2C145/F155  in  a  24-ply  plate  with  the  stacking 
sequence  [-+45/(+-45)2/(-+45)2/(+-45)2/(-+45) 2/+-45 ] . 

4. 4. 1.2  Delamination  of  Multidirectional  Composites 

In  actual  structures,  laminate  layups  are  not 

generally  unidirectional.  Therefore  data  from  unidirectional 

tests  can  only  be  used  directly  in  structural  design  and 

analysis  if  it  can  be  shown  that  the  unidirectional  data 

reflects  the  response  of  the  more  complex  material.  Recent 

studies  show  that  this  is  not  always  the  case.  For  example, 

2 1 

Nicholls  and  Gallagher  have  observed  in  their  DCB  tests 
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Figure  32.  Fracture  energy  versus  crack  length  using  the  area  method  (G) 
and  a  damage-dependent  J  integral  (J)  to  analyze  DCB  data. 

The  material  is  T2C145/F155  with  a  lay-up  of 
[±45/(T45)2/(±45)2/(T45)2/(±45)2/(T45)2/±45!. 
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almost  a  twofold  increase  in  fracture  energy  over  the 
unidirectional  value  when  a  crack  travels  at  a  +30°  or  +45° 
interface  in  ASl/3502  graphite/epoxy.  The  fracture  surfaces 
evidenced  a  change  in  mechanism  which  was  responsible  for  the 
greater  toughness. 

In  this  study  just  mentioned,  several  notable  results 
were  obtained  which  were  reflected  in  the  recommendations  for 
further  work.  First,  it  was  observed  that  when  the  crack  was 
propagating  at  the  interface  of  two  plies,  neither  of  which 
was  parallel  to  the  propagation  direction,  the  crack  tended 
to  jump  between  plies  by  forming  another  crack  parallel  to 
the  original  crack,  but  on  a  different  plane.  They  suggested 
that  a  more  detailed  knowledge  of  the  crack  front  shape  in 
this  region  would  be  helpful  in  understanding  this 
phenomenon.  A  second  subject  for  additional  work  was  the 
transition  between  a  region  where  the  crack  broke  through  a 
ply  and  a  region  where  growth  occurred  between  two  plies,  one 
of  which  had  fiber  parallel  to  the  propagation  direction.  As 
shown  in  Fig.  33,  the  transition  zone  was  V-shaped,  with  the 
first  region  tapering  off  into  the  second.  This  shape 
suggests  an  edge  effect,  and  thus  a  width  effect.  Knowledge 
of  the  stress  state  across  the  width  of  the  specimen  would  be 
needed  to  explain  this  behavior.  Thirdly,  a  wide  range  of 
fracture  toughness  values  was  observed.  There  is  therefore 
reason  to  consider  how  to  interpret  the  results  in  a  way  that 
is  relevant  to  analysis  of  crack  growth  in  structures. 

Given  the  complex  delamination  behavior  of 
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mul  t  i  di  r  ec  t  i  nns  1  composites,  how  can  one  compare  the  results 
from  the  delamination  testing  between  various  layups?  Also, 
how  do  the  toughness  measurements  correlate  with  the  values 
for  the  neat  resin?  By  understanding  the  relationship 
between  the  behavior  of  various  layups,  one  hopes  to  be  able 
to  predict  the  delamination  performance  of  a  proposed  layup 
based  on  the  results  of  a  few  baseline  tests,  rather  than 
needing  to  make  and  test  each  potential  layup.  Also, 
understanding  the  critical  parameters  which  determine 
toughness  would  lead  to  guidelines  for  the  choice  and 

development  of  materials. 

2  2 

Chai  also  used  the  DCB  geometry  to  study 
delamination  at  various  interfaces  in  multidirectional 
layups.  Like  Nicholls  and  Gallagher,  he  found  that 
delamination  energy  correlated  strongly  with  fracture 
morphology.  He  observed  that  various  mechanisms  operated 
during  the  course  of  a  single  test,  and  that  transitions 
between  fracture  mechanisms  occurred  unpredictably .  In 
general,  the  crack  did  not  continue  to  grow  in  its  plane,  but 
shifted  interfaces  as  it  propagated.  Large  differences  in 
the  angles  of  the  plies  at  the  crack  plane  were  noted  to 
encourage  ply  jumping.  Recent  tests  by  Goetz  have 
demonstrated  that  in  the  absence  of  0°  fibers  to  contain  the 
crack,  it  can  leave  the  midplane  and  jump  from  ply  to  ply 
until  finally  reaching  the  surface;  during  the  process  the 
specimen  behavior  is  radically  altered.  Such  behavior  is  not 
seen  in  the  testing  of  unidirectional  composites. 
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Figure  33.  Schematic  of  transition  between  two  regions  on  a  fracture  surface. 
Region  I  corresponds  to  the  crack  front  breaKing  through  a  ply. 
Region  II  appears  as  the  crack  propagates  between  two  dissimilarly 
oriented  plies  one  of  which  is  oriented  parallel  to  the  direction 
of  crack  propagation.  (After  [12].) 
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Chai  observed  local  toughness  changes  on  the  order  of 
100%  due  to  transitions  in  fracture  mechanism.  High  fracture 
energies  were  not  necessarily  due  to  fiber  breakage;  for 
example,  even  when  the  crack  traveled  between  two  0°  plies, 
the  multiplanar  character  of  the  fracture  lead  to  variations 
in  roughness  prompting  his  statement  that,  "The  range  in 
fracture  energies  found  in  the  region  of  the  plies  as 
opposed  to  a  single  value  may  be  cause  for  concern  in  the 
fracture  testing  of  unidirectional  laminates." 

A  major  conclusion  in  Chai's  paper  was  that  the  mode  I 
interlaminar  fracture  energy  was  independent  of  test  specimen 
geometry  (crack  length  and  width),  independent  of  ply 
orientations  at  the  delamination  surface,  and  dependent  only 
on  the  matrix  constituent.  It  should  be  noted,  however,  that 
in  order  to  make  this  statement  Chai  had  to  exclude  from 
consideration  all  regions  of  the  fracture  surface  where  fiber 
breakage,  pull-out,  and  ply-jumping  were  observed. 
Essentially,  by  excluding  fracture  processes  where  the  fibers 
made  a  difference  (on  the  basis  of  fracture  surface 
morphology),  Chai  insured  that  the  matrix  consistence  was 
determinative  in  the  toughnesses  measured.  Indeed,  in  the 
absence  of  such  a  screening  of  data  on  the  basis  of  surface 
morphology,  a  constant  value  of  delamination  toughness  is  not 
obtained.  When  he  compared  the  results  of  other 
investigators  with  the  hypothesis,  it  was  largely 
unsupported.  In  general,  fracture  energy  appears  to  be 
dependent  on  interface  angle  and  potentially  on  fiber 
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DroDerti.es  as  well. 


The  evidence  he  presents  suggests  that  a  less  strong, 
more  precise  conclusion  is  that  minimum  toughness  is  obtained 
for  matrix-dominated  mechanisms,  and  that  similar  toughness 
is  sometimes  obtained  despite  varying  fracture  morphology 
details.  Because  a  lower  bound  gives  the  most  conservative 
result  in  analysis,  this  is  very  significant.  However, 
identifying  a  constant  lower  bound  does  not  fully  explain  the 
data.  Understanding  the  higher  energy  processes  which  lead 
to  higher  delamination  resistance  is  still  necessary. 

Another  aspect  of  the  paper's  conclusion  bears 

scrutiny.  The  delamination  toughness  is  stated  to  be 

independent  of  DCB  width.  The  results  for  specimens  with  a 

range  of  widths  agreed.  For  one  material,  the  width  was 

varied  from  4.2  to  8  mm  and  for  another  from  7.2  to  10.4  mm. 

However,  the  thickness  of  the  48  ply  layups  were  6.1  and  6.7 

mm,  respectively.  Since  edge  effects  would  be  expected  to  be 

2  3 

seen  over  a  distance  on  the  order  of  the  thickness  ,  these 
specimens  were  not  nearly  wide  enough  to  produce  a  change  in 
behavior,  even  if  such  a  transition  does  occur  in  wider 
specimens.  Furthermore,  evidence  of  edge  effects  is  seen  in 
some  of  the  fracture  morphologies  shown  in  the  paper.  Aside 
from  the  effect  of  width,  Chai  cites  the  effect  of  specimen 
geometry  on  the  crack  tip  damage  zone  as  a  complicating 
factor  making  interpretation  of  the  toughness  he  reported 

more  difficult. 

2  4 

Jordan  conducted  delamination  tests  on 


77 


mnlt-irH  rpctional  laminates .  Tests  were  performed  using 

layups  with  0°  fibers  except  at  the  fracture  plane,  at  which 

there  was  a  +10  or  +45  degree  interface.  For  these  cases, 

the  stiffnesses  of  the  laminates  were  similar  to  that  of  the 

unidirectional  layup,  and  the  toughness  values  were  also 

similar.  However,  he  noted  that  as  the  percentage  of 

off-axis  plies  was  increased  and  the  axial  modulus  decreased, 

the  size  of  the  zone  of  microcracking  at  the  crack  tip 

increased,  as  did  the  amount  of  damage  in  the  specimen  away 

from  the  crack  tip.  Analytical  methods  which  did  not  account 

for  the  effect  of  the  far  field  damage  gave  results  showing  a 

large  increase  in  apparent  fracture  toughness.  However,  when 

data  reduction  was  performed  using  a  damage-dependent  J 

integral,  the  intrinsic  fracture  energy  was  found  to  be 

comparable  to  that  for  the  unidirectional  composite  in  three 

of  the  four  systems  studied.  This  suggests,  in  accord  with 

Chai's  conclusions,  that  there  is  an  intrinsic  fracture 

energy  which  is  independent  of  interface  angle.  In  another 
2  5 

recent  study  ,  a  similar  result  was  found  by  applying  the  J 

integral  to  delamination .  Such  layup  insensitivity  of 

toughness  suggests  that  the  local  normal  interface  stress, 

rather  than  local  layup-induced  shear  stress,  is  the  primary 

2  4 

factor  in  delamination  .  It  should  be  noted  that  this 
conclusion  is  based  on  limited  data.  Whether  local  shear 
stress  was  insignificant  or  whether  any  amount  of  local  shear 
stress  would  be  unimportant  is  not  clear.  Jordan's  results 
from  mixed-mode  tests  on  one  composite  suggest  that  when 
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externally-produced  is  over  about  40%  of  the  total  energy 
release  rate  ( GT=»GT +Gj  j  )  ,  GT  will  become  sensitive  to  the 
percentage  of  Gjj  and  will  increase.  Thus,  if  the  results 
from  load-induced  mixed-mode  tests  are  relevant  to  the 
analysis  of  delamination  based  on  the  local  stress  state,  one 
would  expect  local  mode  II  conditions  to  be  evidenced  by  an 
increase  in  toughness. 

The  determination  of  mode  I  delamination  fracture 
toughness  of  multidirectional  composites  poses  numerous 
questions--many  of  which  do  not  exist  for  unidirectional 
composites.  When  performing  delamination  tests,  one  is 
seeking  a  geome t ry- i ndependen t  material  parameter.  The 
question  of  geometry  dependence  has  not  been  adequately 
addressed  for  multidirectional  DCB  specimens.  The  crack 
front  often  takes  the  shape  seen  in  Fig.  34a.  The  change  in 


curvature 

near 

the  edge  seems 

to 

be  an 

edge  effect.  The 

di f  f e  rence 

in 

c rack  length 

a  t 

the 

edge 

and  the 

center  if 

typically 

.15 

to  .2  inches 

for 

a 

one  inch  wide 

specimen . 

Compare  this  to  the  relatively  mild  curvature  seen  for  a 
unidirectional  specimen  (Fig.  34b).  The  difference  in 
curvature  implies  that  the  crack  tip  stress  field  is 
significantly  different  for  the  two  cases.  To  my  knowledoe, 
the  causes  and  significance  of  this  have  n-t  been  addressed 
in  the  open  literature.  However,  related  work  is  currently 
r-emg  done  at  Texas  A&M  University. 

Figure  35  shows  a  record  of  the  load  versus  applied 
:  :  lament  curve  for  a  typical  DCB  test.  The  smooth  curve 
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Load  (N) 


Displacement  (cm) 


Figure  35.  Load  vs.  applied  displacement  record  for  a  typical  DCB 
test  illustrating  stable  crack  growth 
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Load  (lbs) 


T6T145/F 155  Graphite/Epoxy 
i -45/0g/ (p  4 5 ) 2  'Og / ±4  5 ] 


Figure  36  .  Load  vs.  applied  displacement  for  a  DCB  test  illustrating 
unstable  crack  growth 
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is  the  result  of  stable  crack  growth.  In  contrast,  the  curve 
in  Fig.  36  shows  evidence  of  intermittent  unstable  growth 
during  the  displacement-controlled  test.  Such  data  are  more 
difficult  to  interpret  than  those  in  Fig.  35.  Tests 
conducted  on  the  current  grant  have  shown  that  regions  of 
fast  crack  growth  are  generally  regions  of  lower  fracture 
toughness.  What  is  the  significance  of  the  arrest  energy? 
What  is  its  geometry  dependence?  A  second  question  is  raised 
by  the  parts  of  the  curve  which  show  increasing  load  with 
increasing  displacement  coming  from  what  is  evidently  damage 
zone  growth  or  stable  crack  growth.  Is  there  some 
significance  of  this  behavior  analogous  to  that  seen  in  the 
resistance  curves  commonly  determined  for  metals?  Such 
questions  have  yet  to  be  answered. 

As  mentioned  previously,  studies  have  shown  that  the 

fracture  energy  of  a  specimen  can  vary  widely  with  crack 

length  if  fracture  mode  transitions  (e.g.,  ply  jumping)  take 

place.  This  creates  a  problem  in  data  interpretation  and 

reporting.  Chai  addressed  this  by  screening  his  results  to 

isolate  the  lowest  energy  mechanisms,  and  then  reporting  the 

lowest  value  as  the  interlaminar  fracture  toughness.  Often 

data  reduction  techniques  are  used  which  smooth  out  the  data, 

2  5 

such  as  the  area  method  .  There  is  a  need  for  data  analysis 
which  is  able  to  capture  the  details  of  the  various  fracture 
processes . 

It  has  been  reported  that  in  some  instances 
multidirectional  fracture  toughness  is  similar  to  the 
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unidirectional  value  and  in  some  it  is  not.  The  limits  of 

the  generality  of  the  former  statement  have  yet  to  be 

explored.  The  potential  significance  of  such  a  correlation 

is  very  great,  since  it  would  allow  one  to  perform  a  simple 

test  to  obtain  information  for  many  complex  layups. 

Previously,  the  generation  of  far-field  damage  (damage 

away  from  the  crack  tip,  e.g.,  in  the  legs  of  a  DCB  specimen) 

was  cited  as  a  complicating  factor  in  delamination  testing. 

Care  must  be  taken  to  differentiate  between  energy  going  into 

the  fracture  process  and  energy  going  into  generation  of  the 

far-field  damage.  Only  a  handful  of  investigators  have 

attempted  to  develop  analyses  and  tests  which  compensate  for 
2  0  2  5 

it  '  .  The  degree  of  far-field  damage  is  known  to  be 

dependent  on  stacking  sequence.  More  matrix-dominated  layups 
and  layups  with  large  differences  in  interface  angles 
generate  more  damage.  In  a  DCB  specimen,  damage  can  be  due 
not  only  to  bending  stresses,  but  also  to  tensile  stresses  in 
the  legs  near  the  loading  points  when  the  angle  of  rotation 
of  the  legs  is  large. 

Besides  affecting  the  introduction  of  far-field 

damage,  the  stacking  sequence  can  potentially  affect  the 

delamination  process  by  producing  local  shear  and  twisting  at 

the  crack  tip,  which  may  lead  to  partial  crack  closure  or 

2  6 

cause  local  mode  II  or  mode  loading 

4. 4. 1.3  Objectives  of  Current  Study 

The  primary  objective  of  the  present  study  has  been  to 
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develop  a  theoretically-based  procedure  for  determining  the 
mode  I  delamination  fracture  toughness  of  multidirectional 
composite  materials. 

The  nature  of  the  questions  to  be  answered  was 
expressed  in  the  previous  section.  Here  are  some  of  the 
fundamental  questions  underlying  the  objective  of  the 
proposed  research: 

o  Under  what  conditions  can  a  geometry-independent 
measure  of  fracture  energy  be  obtained? 

o  What  is  the  cause  and  significance  of  crack  front 
cu  rvatu  re  ? 

o  what  data  analysis  is  necessary  to  capture  the 
details  of  the  delamination  process? 

o  What  is  the  significance  and  geometry  dependence  of 
the  intermittent  growth  reflected  in  data  such  as 
seen  Fig.  34?  What  is  the  significance  of  the 
increasing  crack  growth  resistance? 

o  Under  what  conditions  do  various  fracture 
morphologies  occur? 

o  How  do  observed  toughness  values  relate  to  the 
corresponding  fracture  morphologies? 

o  For  a  given  crack  tip  damage  state  (constant 

intrinsic  fracture  energy),  what  is  the  effect  of 
varying  stacking  sequence,  i.e.,  changing  the 
amount  of  far-field  damage? 

o  What  is  the  effect  of  changing  the  crack  tip  stress 
state  through  altering  the  stacking  sequence? 


o  How  does  damage  due  to  tensile  stress  (as  opposed 
to  bending  stress)  impact  the  apparent  toughness? 
o  How  does  the  fracture  energy  depend  on  crack  speed? 

4. 4. 1.4  Experimental  Investigation 

To  date,  the  focus  of  the  work  has  been  on  application 
of  the  J  integral  to  double  cantilever  beam  ( DCB )  tests  in 
order  to  investigate  the  importance  of  geometry,  layup,  rate 
effects,  and  fracture  morphology.  The  goal  is  to  synthesize 
principles  which  can  guide  the  prediction  of  delamination 
performance  using  limited  material  data.  The  following 
discussion  highlights  the  method  used  and  results  obtained  to 
date . 

4. 4. 1.4.1  Materials  and  Layups — The  material  used  in  this 
study  is  T2C145/F155,  manufactured  by  Hexcel .  It  uses  a 
rubbe r- toughened  epoxy  matrix  with  about  6%  rubber  by  weight. 
The  nominal  fiber  volume  fraction  for  all  the  material  tested 
was  57%.  This  material  system  was  chosen  in  order  to 
accentuate  possible  non-linear  and  viscoelastic  responses  for 
a  commercially  available,  toughened  composite.  It  is  also 
currently  being  used  for  other  work  on  this  contract  and 
other  contracts  at  Texas  A&M,  allowing  sharing  of  material 
data.  Three  layups  were  used.  Hereafter,  they  will  be 
referred  to  using  the  following  designations: 

Designation _ Layup _ No.  Plies 

unidirectional  24 
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fiber-dominated  ( +-4 5/0 g/-+4 5 ] ant i sym  24 

angle-ply  [ +-45/( -+45 ) 2/( +-45 ) 2/-+45 ] antisym.  24 
The  unidirectional  layup  is  the  one  which  has  been  most 
commonly  used  in  the  literature  to  characterize  delamination. 
In  this  study  it  was  used  not  only  to  examine  delamination 
between  zero  degree  plies  where  ply  interpenetration  could 
occur,  but  to  give  a  baseline  for  evaluating  data  for 
multidirectional  layups.  The  fiber-dominated  layup  was  used 
to  study  delamination  at  a  +-45  degree  interface.  Each  leg 
of  the  DCB  specimen  was  made  balanced  and  symmetric  in  order 
to  eliminate  stretching-shearing  and  stretching-bending 
coupling.  Antisymmetry  about  the  midplane  of  the  laminate 
was  designated  in  order  to  put  the  delamination  plane  at  a 
+-45  degree  interface.  The  stacking  sequence  for  the 
angle-ply  layup  was  chosen  to  minimize  bending-twi sting 
coupl ing  as  well . 

4. 4. 1.4. 2  Experimental  Results--All  data  shown  in  the 

following  figures  were  generated  using  DCB  specimens.  The 

fracture  toughness  was  calculated  using  the  J  integral  method 
2  7 

of  Schapery  .  In  this  formulation  J  is  calculated  using  the 
moment  at  the  crack  tip  during  crack  extension  and  the 
moment-curvature  relationship  for  one  leg  of  the  specimen. 

The  J  integral  is  twice  the  area  to  the  left  of  the 
moment-curvature  curve,  as  shown  in  Fig.  37.  In  the 
following  discussion,  J  will  correspond  to  this  path  integral 
definition.  The  energy  release  rate  was  also  calculated  for 
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comparison.  The  symbol  G  will  signify  the  energy  release 
rate  calculated  by  the  area  method^.  The  area  method 
derives  from  the  derivative  definition  of  energy  release 
rate,  and  provides  an  average  value  for  an  increment  of  crack 
growth.  The  symbol  will  refer  to  the  energy  release  rate 
calculated  using  the  equation 

3P„ 
c  c 

G,  =  -  (3) 

2ba 

where  P  is  the  load  at  crack  extension,  S  is  the  load  line 
c  c 

displacement  at  crack  extension,  b  is  the  specimen  width,  and 
a  is  the  crack  length.  This  equation  is  based  on  the 
assumptions  of  linear  beam  theory  and  material  linearity. 

Comparison  of  Jc  and  Gc  for  Delamination  of  Angle  Ply 
Composites — Motivation  for  using  a  J  integral  comes  from  the 
need  for  an  analysis  which  would  account  for  continuum  damage 
such  as  might  be  generated  throughout  the  legs  of  a  DCB 
specimen  during  testing.  The  area  method  of  calculating  the 
energy  release  rate  allows  for  geometric  and  material 
nonlinearity,  as  does  J,  but  does  not  differentiate  between 
energy  which  goes  into  driving  the  delamination  and  that 
which  causes  general  damage  away  from  the  crack  tip,  leading 
to  an  overestimation  of  the  delamination  fracture  toughness. 
The  difference  is  illustrated  in  Fig.  32,  where  G  and  J  are 
plotted  versus  crack  length  for  three  angle-ply  DCB 
specimens.  Note  that  for  each  specimen,  G  is  significantly 
higher  than  J.  The  presence  of  damage  is  illustrated  by  the 
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p  moment 


Figure  37.  Double  cantilever  beam  specimen,  integration  path,  and  moment 
curvature  diagram  for  calculation  of  J  [after  1] 


moment-curvature  relationship  for  the  angle-ply  layup,  as 
shown  in  Fig.  38.  Note  that  the  unloading  portion  of  the 
moment-curvature  plot  does  not  retrace  the  loading  portion. 

For  layups  containing  a  high  percentage  of  0  deqr-e 
fibers,  one  would  expect  that  G  and  J  would  give  better 
agreement,  since  continuum  damage  would  be  minimized.  The 
moment-curvature  relationships  for  the  unidirectional  and 
fiber-dominated  layups  were  linear,  and  the  unloading 
portions  of  the  curves  retraced  the  loading  portions.  When 
this  is  the  case  and  linear  beam  theory  applies,  it  can  be 
shown  that  J  is  equal  to  calculated  using  Eq.  3.  This 
formulation  is  helpful  when  comparing  energy  release  rate  and 
J  because  G^  is  calculated  for  a  particular  point  (as  J  is 
for  the  linear  and  nonlinear  cases),  avoiding  the  averaging 
of  the  area  method.  Figure  39  is  the  load-displacement 
record  from  a  fiber-dominated  specimen  which  displayed  the 
development  and  breakdown  of  a  tie  zone,  causing  the 
delamination  resistance  to  increase  and  then  suddenly  drop 
with  crack  growth.  The  corresponding  G^  and  J  are  plotted  in 
Fig.  40.  Both  initiation  and  arrest  values  are  shown.  The 
two  methods  are  in  reasonably  close  agreement,  supporting  the 
interpretation  of  the  difference  between  energy  release  rate 
and  J  for  the  angle-ply  layup  as  being  caused  by  damage. 

Even  when  damage  is  not  present,  the  simple  expression  in  Eq . 
3  does  not  always  hold.  When  geometric  nonlinearity  is 
present  due  to  large  rotations  of  the  legs,  the  energy 
release  rate  must  be  calculated  using  either  the  area  method 
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(possibly  losing  information  due  to  averaging)  or  using 
nonlinear  beam  theory^. 

Geometry  Dependence  of  J — The  usefulness  of  J  as  a 
characterizing  parameter  depends  on  its  independence  of 
geometry,  such  as  specimen  width.  To  investigate  the  effect 
of  specimen  width,  angle-ply  specimens  1/2,  1,  and  2  inches 
wide  were  tested.  Representative  results  are  given  in  Fig. 
41.  The  high  values  for  the  1  inch  wide  specimen  for  shorter 
crack  lengths  were  found  to  be  associated  with  the  complex 
fracture  morphology  which  developed  near  the  starter  crack. 
For  more  nearly  self-similar  crack  advance,  it  appears  that 
width  does  not  have  much  effect  on  the  fracture  toughness 

measured  over  the  range  of  widths  and  crack  lengths  studied. 

2  8 

Recent  work  at  Texas  A&M  indicates  that  the  state  of  stress 
in  the  legs  of  the  DCB  specimen  is  in  the  transition  range 
between  plane  strain  and  plane  stress  over  the  practical 
range  of  crack  length  to  width  aspect  rations.  However,  the 
effects  of  the  transition  on  J  were  not  seen  for  the  aspect 
ratios  used  in  our  tests.  However,  there  is  another  aspect 
of  specimen  behavior  which  does  show  a  marked  width  effect. 
Unlike  unidirectional  composites,  multidirectional  composites 
exhibit  significantly  curved  crack  fronts  in  DCB  tests.  The 
curvature  depends  on  the  width  of  the  specimen.  Figure  42 
compares  the  normalized  crack  front  profile  for  1/2  inch,  1 
inch,  and  2  inch  wide  specimens  for  several  crack  lengths. 
Apparently  the  crack  front  curvature  has  two  sources:  the 
antielastic  curvature  of  the  DCB  leg  and  free  edge  effects. 
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PARAMETER 


Figure  40.  Gi  and  J  vs.  crack  length  for  the  fiber-dominated  specimen  of 
figure  4.4.  Arrest  values  are  marked  by  squares  around  data 
points.  The  layup  Is  I +-45/0g/(-+45)2/0g/-+45| . 
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Figure  41,  g  and  J  vs.  crack  length  for  three  angle-ply  specimens 


The  general  trend  appears  to  be  decreasing  curvature  with 
increasing  crack  length,  but  there  is  considerable  scatter  in 
the  behavior.  The  differences  between  the  shapes  of  the 
three  plots  are  plain.  The  curvature  associated  with  edge 
effects  is  dominant  for  the  1/2  inch  specimen,  but  becomes 
less  significant  for  increasing  widths.  These  profiles  were 
taken  from  fracture  surfaces  where  the  crack  arrested  after 
jumping  forward,  marking  the  surface. 

Effect  of  Stacking  Sequence  on  J — The  effect  of  layup 

on  J  was  studied.  It  was  found  that  when  complex  mechanisms 

of  fracture  involving  fiber  breakage,  ply  jumping,  extensive 

fiber  bridging,  etc.,  were  absent,  the  results  for  various 

layups  were  often  comparable.  However,  multidirectional 

layups  often  displayed  such  mechanisms.  Other 
21  22 

investigators  '  have  reported  similar  behavior.  Figure  43 

illustrates  the  variety  of  behaviors  observed.  The 

unidirectional  specimen  had  a  uniform  toughness  over  a  large 

range  of  crack  lengths.  The  fiber-dominated  layup  led  to 

large  variations  in  J  associated  with  tie  zone  development 

and  breakdown.  Typically  the  fiber-dominated  layup  produced 

data  with  the  most  scatter.  Data  interpretation  was  often 

difficult.  Such  large  variations  in  toughness  for  similar 

21 

layups  have  been  reported  elsewhere  .  Practically,  this  is 
unfortunate,  since  one  would  like  to  use  this  layup  to  study 
the  effect  of  various  fiber  angles  at  the  delamination  plane 
without  resorting  to  making  measurements  of  the 
moment-curvature  relationship  (since  the  relationship  is 
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1/2  INCH  WIDE  SPECIMEN 


1  INCH  WIDE  SPECIMEN 


Figure  42.  Normalized  crack  front  profiles  for  three  angle-ply  DCB 

specimens  of  different  widths.  The  profiles  for  various  crack 
lengths  have  been  superimposed. 
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2  INCH  WIDE  SPECIMEN 


Figure  42.  (continued) 
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linear  so  that  Eq .  3  can  be  used).  The  angle-ply  layup, 
being  a  "softer"  layup,  and  with  distributed  (energy 
absorbing)  damage,  exhibited  less  crack  jumping  and  scatter 
in  the  data.  However,  as  shown  in  Fig.  43,  complex  fracture 
mechanisms  were  sometimes  present  early  in  the  test  before 
the  delamination  plane  had  been  established. 

The  effect  of  fracture  morphology  is  further 
illustrated  in  Fig.  44  where  features  on  the  fracture 
surfaces  are  related  to  the  toughness  measured.  Note  that 
once  a  more  or  less  uniform  fracture  morphology  was 
established,  the  toughness  values  settled  into  better 
agreement.  The  two  angle-ply  specimens  shown  were  next  to 
each  other  in  the  plate  of  material  before  cutting.  Dark  and 
light  bands  were  formed  when  the  crack  tip  jumped  forward, 
then  arrested.  Dark  regions  were  formed  during  rapid  crack 
advance;  light  regions  correspond  to  slow  crack  growth. 

These  studies  of  the  effect  of  layup  illustrate  a  major 
difference  between  unidirectional  and  multidirectional 
laminate  delamination  behavior.  The  opportunity  for  more 
mechanisms  of  fracture  leads  to  a  greater  complexity  of 
behavior.  It  is  therefore  essential  that  attention  be  paid 
to  the  fracture  surface  morphology  when  interpreting  the 
data . 

Effect  of  Crack  Speed  on  J — Tests  were  done  to 
investigate  the  dependence  of  J  on  crack  speed.  The  results 
for  a  unidirectional  layup  and  a  fiber-dominated  layup  are 
given  in  Fig.  45.  In  both  cases  there  is  a  slight  decreasing 
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Figure  43.  j  vs.  crack  length  for  three  different  layups 
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Figure  44.  Correlation  of  fracture  surface  morphology  with  measured 

delamination  fracture  toughness  for  two  angle-ply  specimens 
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trend  in  J  with  increasing  crack  speed.  This  trend  agrees 
with  the  observation  that  when  sudden  crack  jumping  occurred, 
a  low  value  of  toughness  was  usually  measured.  It  is 
significant  to  note  that  no  strong  trend  can  be  established 
for  the  range  of  crack  speeds  experienced  by  one  specimen. 
Thus,  rate  effects  are  not  important  in  interpreting  data  for 
individual  specimens  in  this  material  tested  at  ambient 
tempe  rature . 

Fracture  Surface  Morphology  and  J — The  importance  of 

considering  the  fracture  morphology  when  interpreting  the 

delamination  toughness  data  has  already  been  mentioned. 

Macro-  and  mi c r o f r ac tog r aphy  were  done  to  document  the 

fracture  mechanisms  which  were  observed.  Two  items  are  of 

particular  interest.  First,  when  off-axis  plies  were  present 

at  the  fracture  plane,  toughness  values  were  sometimes 

elevated  even  in  the  absence  of  complex  fracture  mechanisms. 

Apparently  this  was  due  to  increased  surface  roughness  caused 

when  bundles  of  fibers  were  pulled  away  from  the  fracture 

surface.  This  feature  was  only  observed  for  layups  with 

off-axis  plies  at  the  delamination  plane.  For  a  more  brittle 

system  the  phenomenon  was  found  to  give  differences  of  a 

factor  of  two  in  toughness  between  unidirectional  and 

2  5 

off-axis  specimens  .  Figure  46  is  a  cross-section  of  a 
specimen  (normal  to  the  direction  of  crack  propagation)  at 
the  fracture  surface.  A  bundle  of  fibers  is  separated  from 
the  surface  at  this  plane,  but  is  attached  to  it  at  another 
point . 
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A  second  phenomenon  which  is  of  significance  is  that 
in  almost  all  specimens  tested,  delamination  was  found  to 
occur  within  a  ply  rather  than  between  two  plies.  Figure  47 
shows  a  typical  situation  on  a  polished  cross  section  at  a 
fracture  surface.  Note  that  the  fracture  plane  does  not  pass 
through  the  interply  region,  but  stays  a  few  fiber  diameters 
away  from  it.  Practically,  the  presence  of  resin-rich 
regions  and  ply  interfaces  has  only  an  indirect  impact  on  the 
mode  I  delamination  toughness  when  such  steady-state 
morphology  is  present. 

Future  work  under  the  follow-on  grant  is  expanding  on 
these  investigations  of  the  effect  of  geometry,  layup,  rate 
effects  and  fracture  surface  morphology.  Further  experiments 
are  being  performed  to  explore  the  effect  of  specimen  width, 
especially  as  it  relates  to  crack  front  curvature.  The 
effect  of  the  number  of  plies  and  the  bending  stiffness  of 
the  legs  of  the  DCB  specimen  is  being  addressed.  Also, 
tests  are  being  conducted  with  various  fiber  angles  at  the 
delamination  plane.  The  dependence  of  delamination  toughness 
on  crack  speed  for  matrix-dominated  layups  is  also  being 
investigated  over  a  wider  range  of  displacement  rates.  An 
attempt  will  be  made  to  systematize  the  correlation  between 
fracture  surface  morphology  and  delamination  toughness.  In 
addition  to  these  continuing  activities,  new  studies  are 
beginning  to  improve  our  understanding  of  the  effect  of 
continuum  damage  on  the  apparent  toughness.  This  will 
probably  include  delamination  tests  which  simulate  mechanical 
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Specimen 


Figure  46.  J  vs.  crack  speed  for  fiber-dominated  DCB  specimens.  The 
layup  Is  |+-45/08/(-+45)2/08/-+45 | 


Figure  47a  Micrograph  of  a  polished  section  of  a  failed  DCB  specimen 
showing  the  fracture  surface.  A  bundle  of  fibers  has  been 
pulled  away  from  the  surface.  (150X) 


Figure  47b.  Micrograph  of  a  polished  section  of  a  failed  DCB  specimen 

showing  the  fracture  surface.  The  large  circles  at  the  top  are 
the  material  used  to  mount  the  specimen  for  polishing.  (100X) 


states  in  real  structures,  such  as  tests  using  a  modified  DCB 
specimen  with  bending  and  stretching  inputs. 

4.4.2  Nixed  Mode  and  Mode  II  Delamination  of  Composite 

Materials 

Our  efforts  in  the  characterization  of  mixed  mode  and 
mode  II  delamination  of  composite  materials  have  fallen  into 
three  categories:  (1)  a  comparison  of  the  end  notch  flexure 
and  the  end  loaded  flexure  tests  to  see  if  they  give 
consistent  values  of  GIIc;  (2)  development  of  a  J-integral 
approach  to  evaluate  mode  II  delamination  fracture  toughness 
for  composites  made  of  resins  which  have  significant 
nonlinear  behavior  and/or  split  laminate  specimens  which  give 
significant  geometric  nonlinearity  during  testing;  and  (3) 
determination  of  the  G  total  as  a  function  of  the  percentage 
of  Gjj  as  well  as  the  state  of  stress  ahead  of  the  crack  tip 
for  mixed  mode  loading.  The  results  of  our  efforts  in  these 
areas  will  be  summarized  in  the  three  subsections  that 
follow, 

4. 4. 2.1  Comparison  of  the  End  Notch  Flexure  Test  and  the 

End-Loaded  Split  Laminate  Test  for  Delamination  of 
Composite  Materials 

Several  test  configurations  have  been  proposed  for 
measuring  the  resistance  to  delamination  under  mode  II 
loading,  including  the  end-notched  flexure  (ENF)  test  and  the 
end-loaded  split  laminate  (ELS)  test.  The  end-notched 
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)  ENF  test  configuration  and  associated  parameters. 


p°i 


(b)  ELS  test  configuration  and  associated  parameters. 


Figure  48.  Mode  II  delamination  tests. 
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flexure  test  configuration  can  be  seen  in  Fig.  48a.  This 
test  consists  of  a  laminate  beam  with  a  split  on  one  end, 
loaded  in  three  point  bending.  The  end  loaded  split  laminate 
test  consists  of  a  split  laminate  beam  specimen  which  is 
loaded  at  the  cracked  end  and  fixed  at  the  opposite  end  (see 
Fig.  48b).  The  experimental  results  from  the  ENF  and  ELS 
test  methods  may  be  evaluated  to  give  Gj  using  equations 
derived  using  linear  beam  theory;  i.e.,  assuming  geometric 
linearity  and  linear  elastic  material  behavi o r 3 ®  '  3 3  .  A 
review  of  the  current  literature  indicates  that  both  test 
configurations  are  currently  being  used  to  measure  GIIC3^_3^. 

During  a  recent  finite  element  analysis33,  the  ENF 
test  was  found  to  be  a  pure  mode  II  fracture  test  within  the 
constraints  of  small  deflection  theory.  The  study  also 
revealed  that  the  interlaminar  normal  stress  is  identically 
zero  along  the  beam  center  line  and  the  interlaminar  shear 
stress  exhibits  the  expected  singularity.  A  similar  analysis 
was  performed  for  the  ELS  test  configuration  in  this  project 
using  the  code  described  in  4.3.  The  results  which  were 
presented  in  section  4. 3.1.2  also  indicate  pure  mode  II 
stressing  along  the  center  line. 

4. 4. 2. 1.1  Analysis 

End-Notched  Flexure  (ENF)  Test — The  ENF  test  consists 
of  loading  a  split  laminate  beam  specimen  using  a  three  point 
bend  fixture.  Figure  48a  shows  the  ENF  test  configuration 
and  associated  test  parameters.  Calculation  of  the  strain 
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Syy  or  Sxy  (GPa) 


Figure  49.  Stress  distribution  ahead  of  the  crack  tip  for 
mode  I  and  mode  II  loading. 
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energy  release  rate  Gjj  assuming  linear  beam  theory  and 

linear  elastic  material  behavior  gives  the  following 
i  -34 

result  , 


2  2 

9PzazC 


II 


2b  ( 2L  3  +  3a  3 ) 


(  4  ) 


where  a  is  the  initial  crack  length,  b  is  the  width,  C  is  the 
measured  compliance,  L  is  half-span  length,  and  P  is  the 
applied  load  at  the  center  pin.  The  strain  energy  release 
race  Gj  can  also  be  obtained  by  means  of  an  experimental 
method  for  determining  the  relationship  between  compliance 
and  crack  length.  In  this  case,  GIX  is  given  by34 
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2b 
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(5) 


where  m  is  obtained  from 

C  -  ma3  +  const.  (6) 

This  equation  is  obtained  by  making  a  least  squares  fit  of  a 
straight  line  to  experimental  data  from  measured  compliance 
against  crack  length  cubed.  This  data  is  obtained  from 
load-deflection  curves  at  different  crack  lengths  of  one  of 
the  tested  coupons  of  width  bQ . 

End-Loaded  Split  Laminate  (ELS)  Test — The  energy 
release  rate  for  pure  mode  II  in  this  case  is  determined  by 
asymmetrically  loading  a  spl i t- lami nate  beam  specimen. 


Ill 


Figure  48b  shows  the  ELS  test  configuration  and  associated 


test  parameters.  Linear  beam  theory  and  linear  elastic 
material  behavior  are  assumed  by  the  analysis.  A  complete 
explanation  of  this  analysis  can  be  found  in  reference  31. 
From  this  analysis,  Gjj.  can  be  expressed  by 
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II 
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P  is  the  load  applied  at  the  cracked  end  of  the  beam  (see 
Fig.  48b).  The  modulus  E.^  measured  in  bending  can  be  shown 
to  be 
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where  6  is  the  total  beam  tip  deflection. 

Area  method — In  all  the  methods  of  data  reduction 
previously  described,  the  critical  energy  release  rate  was 
calculated  by  knowledge  of  the  critical  load  Pc  at  onset  of 


crack  growth 

,  and  the 

crack  length  'a'  prior  to 

crack 

extension . 

Therefore , 

Gtt  is  a  measure 

1 1  c 

of  the 

energy 

required  per 

unit  area 

of  crack  growth  to 

initiate 

crack 

propagation . 

However,  as 

the  crack  propagates, 

G T  -  does  not 
lie 

necessarily  remain  constant,  in  which  case  distinct 
propagation  and  arrest  values  of  Gjj  can  exist. 

Where  nonlinear  load-displacement  curves  are  observed 
in  mode  II  testing,  the  above  relationships  based  on  linear 


analysis  are  no  longer  accurate.  An  alternate  approach  for 
evaluating  the  fracture  toughness  for  nonlinear  behavior  is 


the  area  method.  In  this  case,  Gjj  is  still  the  energy 
required  to  create  a  new  crack  per  unit  area  of  crack  surface 
created.  can  then  be  given  by 
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U 
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U  is  the  area  between  the  load-deflection  curve  for  loading 
and  unloading  for  a  very  small  change  Aa  of  crack  growth  (see 
Fig.  50).  An  important  advantage  of  this  method  is  that  only 
elastic  material  behavior  is  required  to  predict  G^. 
Therefore,  for  geometrically  non-linear  and/or  non-linear 
elastic  material  responses,  this  method  gives  an  average 
release  rate  for  the  observed  crack  extension.  For  unstable 
crack  growth,  it  gives  an  average  value  for  GIIc  which 
typically  falls  between  the  initiation  and  arrest  values 
measured  for  linear  behavior.  The  load-displacement  record 
should  return  to  the  origin  to  guarantee  that  no  significant 
far  field  damage  is  included  in  U  in  Eq .  9.  The  inclusion  of 
energy  dissipation  in  far  field  damage  in  the  fracture  energy 
term  U  in  Eq .  9  would  give  an  erroneously  high  estimate  of 
G j j c .  For  materials  where  G^  is  independent  of  crack  growth 
rate  and  crack  growth  increment  (i.e.,  systems  with  minimal 
fiber  bridging  and/or  plastic  wake),  the  average  and 
initiation  Gjj  should  be  identical.  Furthermore,  where  slow 
stable  crack  growth  occurs,  the  average  Gjj  calculated  from 


<S 

Figure  50. Alternate  interpretation  of  Gllc  (area  method) 


114 


the  area  method  should  equal  the  Guc  for  initiation 
calculated  from  linear  beam  theory.  Another  advantage  of 
this  method  is  that  Gjj  is  obtained  without  knowledge  of  the 
materials's  elastic  properties  and  it  is  always  given  by  Eq. 
9,  regardless  of  the  test  configuration  used. 

4. 4. 2. 1.2  Experimental  Procedure 

Split  laminate  beam  specimens  approximately  152.4  mm 
long,  25.4  mm  wide,  and  4.83  mm  thick  for  the  end-notched 
(ENF)  test  and  292.1  mm  long,  25.4  mm  wide,  and  3.05  mm  thick 
for  the  end-loaded  split  laminate  (ELS)  test  were  cut  from  24 
ply  laminated  panels  of  Hercules  AS4/3502  and  Hexcel 
T6C145/F185  graphite/epoxy  composites.  The  starting  crack 
was  introduced  in  both  panels  by  putting  a  teflon  insert 
0.127  mm  thick  at  the  midplane  during  layup.  Loading  tabs 
were  bonded  to  the  cracked  end  of  the  ELS  test  specimens  with 
a  two-part  structural  adhesive.  A  Tenslac  UN-10A  gray 
brittle  coating  was  also  applied  to  the  edges  of  the  specimen 
to  facilitate  measurement  of  crack  length  by  making  the 
location  of  the  crack  tip  more  evident. 

End-Notched  Flexure  Testing — Five  specimens  for  each 
graphite/epoxy  system  were  tested  using  a  three  point  bend 
fixture  with  outer  loading  points  101.6  mm  apart  attached  to 
a  closed  loop  servo-hydraulic  test  machine.  Prior  to 
testing,  the  specimens  were  precracked  in  order  to  provide  a 
sharp  crack  tip  from  which  to  initiate  the  mode  II  fracture. 
Load  point  displacements  were  measured  by  the  ram 
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displacement  and  loads  were  monitored  by  a  2225  N  load  cell. 
Real-time  analog  plots  of  the  load-deflection  curve  were  made 
on  an  x-y  recorder.  The  tests  were  conducted  by  positioning 
the  crack  tip  approximately  midway  between  the  center  and  one 
of  the  outer  loadings  pins.  Under  displacement  control 
conditions  and  a  ram  rate  of  2.54  mm/min,  the  specimens  were 
loaded  until  delamination  crack  growth  occurred.  Then,  the 
specimens  were  unloaded. 

Since  machine  stroke  displacement  was  used  to  measure 
specimen  deflection,  a  machine  compliance  experiment  was 
performed  to  make  the  necessary  deflection  corrections  to  the 
data  obtained  during  the  mode  II  delamination  tests. 

Gtt  calculations  were  made  using  Eq .  4  and  measured 

compliances  after  appropriate  machine  compliance  corrections 

had  been  made.  It  should  be  noted  that  corrections 

associated  with  shear  compliance  were  not  taken  into 

consideration  because  no  reliable  values  of  the  shear  modulus 

G12  were  available.  However,  these  corrections  are  usually 

small  and  the  omission  was  not  thought  to  have  a  significant 

effect  on  the  measured  GTT  values.  The  measured  compliances 

lie 

were  obtained  by  curve  fitting  the  best  straight  line  to  the 
load-displacement  data  from  the  loading  part  of  the  test. 
Equations  5  and  6  were  also  used  to  calculate  GTT  .  To  use 
these  two  equations,  a  compliance  calibration  was  performed 
with  one  of  the  tested  specimens.  Finally,  Gllc  values 
calculated  by  the  area  method  were  determined  by  means  of  Eq . 
9.  The  energy  absorbed  in  the  creation  of  a  new  cracked 
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surface  was  approximated  numerically  by  integrating  the  area 
under  the  load/deflection  curves  obtained. 


End-Loaded  Split  Laminate  Testing — Two  specimens  for 
each  graphite/epoxy  system  chosen  for  this  study  were  tested 
in  a  closed  loop  servo-hydraulic  machine.  A  mode  I  precrack 
was  introduced  to  the  specimens  to  sharpen  the  initial  crack 
created  by  the  teflon  insert  and  to  provide  an  initial  crack 
length  to  beam  length  (a/L)  ratio  of  approximately  0.55, 
where  stable  crack  growth  is  expected^. 

Mode  II  delamination  testing  was  done  under 
displacement  control  at  a  rate  of  10.16  mm/min  and  6.35 
mm/min.  Loads  were  recorded  with  a  445  N  load  cell. 
Displacements  were  monitored  using  the  ram  displacement  of 
the  test  machine.  The  tests  were  conducted  by  loading  the 
specimens  until  the  crack  was  allowed  to  grow  approximately 
10  mm  to  38  mm.  At  that  point  the  test  was  stopped  and  the 
new  crack  location  marked  and  recorded.  This  procedure  was 
repeated  several  times  with  the  same  specimen  until  the  crack 
tip  was  approximately  12.7  to  25.4  mm  from  the  uncracked  end. 

The  critical  energy  release  rate  GIIc  was  calculated 
by  means  of  Eqs.  7  and  9.  The  modulus  was  calculated 
from  the  load-deflection  test  data  and  Eq  3.  The  energy 
absorbed  in  the  creation  of  a  new  cracked  surface  was 
approximated  numerically  by  integrating  the  area  under  each 
loop  in  the  load-deflection  curves. 

4. 4. 2. 1.3  Results  and  Discussion 
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End  Notched  Flexure  (ENF)  Test  Results — Table  2  shows 
the  ENF  test  results  for  both  graphite/epoxy  systems 
investigated.  Method  I  corresponds  to  GIIc  values  obtained 
using  beam  theory  equations  ( Eq .  4)  and  measured  compliances. 
Method  II  results  are  based  on  the  compliance  calibration 
(Eqs.  5  and  6).  Method  III  refers  to  data  reduction  based  on 
the  area  method  (Eq.  9)  and  the  energy  dissipated  under  the 
load-deflection  curve.  Delamination  crack  growth  occurred  at 
typical  load  values  of  556  to  612.5  N  for  AS4/3502 
graphite/epoxy  and  was  unstable  for  this  material.  Figure  51 
shows  a  typical  load-deflection  curve  of  AS4/3502  graphite 
epoxy  showing  the  sudden  drop  of  load  as  unstable 
delamination  crack  growth  occurs.  This  indicates  that  the 
GIIc  values  measured  using  the  area  method  will  be  an  average 
of  the  initiation  value  and  the  value  for  rapid,  unstable 
crack  growth.  The  very  small  nonlinearity  just  prior  to 
crack  extension  is  probably  associated  with  the  development 
of  a  crack  tip  damage  zone. 

The  average  critical  energy  release  rate  GIIc  for 

AS4/3502  using  methods  I,  II,  and  III  of  data  reduction  were 

2  2  2 
560  J/m  ,  648  j/m  ,  and  613  J/m  ,  respectively. 

Figure  52  shows  a  typical  load-deflection  curve 
obtained  for  mode  II  delamination  testing  of  T6C145/F185 
graphite/epoxy.  Crack  growth  was  stable  as  seen  in  the  upper 
part  of  the  curve.  However,  the  point  at  which  crack  growth 
occurred  is  not  well  defined.  Therefore,  the  maximum  load 
reached  during  each  test  was  used  to  estimate  the  onset  of 
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Table  2.-  End-notched  Flexure  (ENF)  Test  Results. 


GIIc  (J/m2) 


Material 

Specimen 

I 

II 

III 

AS4/3502 

1 

X 

525 

643 

613 

2 

508 

630 

595 

3 

560 

682 

613 

4 

5J3 

648 

630 

Av. 

543 

648 

613 

★  * 

T6C145/F1S5  1  -  -  2240 

2  -  -  2695 

3  -  -  2435 

4  _  1995 

Av.  -  -  2398 


*  No  ca'cj'azion  of  3jT-  is  pcsio’e  due  to  significant  nonlinearity  in 
load/def lect ion  curves.  lr  one  uses  linear  portion  of  curves  to  measure 
compliance  and  ?n)ax  for  crack  extension  (ignoring  nonlinear  behavior  near 
?max)’  G 1 1  c  bound  va’ues  of  1225  may  be  calculated. 
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crack  growth.  These  values  were  between  1390  and  1632  N . 

The  curve  also  indicates  inelastic  behavior  of  the  material 

(i.e.,  the  unload  curve  did  not  return  to  the  origin). 

Therefore,  the  assumptions  of  the  test  have  been  violated. 

It  should  be  noted  that  the  load-deflection  curves  at 

different  crack  lengths  used  for  the  compliance  calibration 

also  showed  non-linear  behavior.  However,  the  energy  release 

rate  obtained  by  the  area  method  gives  a  better  approximation 

of  the  critical  energy  release  rate  Gjj  because  at  least  the 

geometrically  non-linear  behavior  of  the  specimen  is  taken 

into  consideration  by  this  method.  The  calculated  average 

2 

GIIc  of  2400  J/m  should  be  considered  an  upper  bound 
estimate  of  the  mode  II  critical  energy  release  rate  of 
T6C1 45/F185  because  the  inelastic  behavior  is  not  accounted 
for  in  this  analysis.  A  J-integral  approach  to  be  discussed 
in  the  next  section  has  been  developed  to  obtain  a  more 
meaningful  measure  of  the  mode  II  delamination  fracture 
toughness  of  this  system. 

End-Loaded  Split  Laminate  (ELS)  Test  Results — The  ELS 
test  results  are  seen  in  Table  3.  Method  I  corresponds  to 
Gj j c  calculated  using  Eq.  7  and  with  calculated  from  Eq. 

8  and  test  data.  Method  II  corresponds  to  the  area  method 
(Eq.  9).  An  average  E ^  of  125.5  GPa  for  AS4/3502  and  109.9 
GPa  for  T6C145/F185  were  calculated.  Figure  53  shows  a 
typical  load-deflection  curve  obtained  by  testing  one  of  the 
AS4/3502  specimens.  As  seen  in  this  figure,  unstable  crack 
growth  is  observed  for  the  first  crack  extension.  This  is 
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LOAD  (N) 


DEFLECTION  (ran) 


Figure  53.  Typical  mode  II  delaminat ion  load/deflection  record  of 
ASL  '3502  using  the  end-loaded  split  laminate  test. 


Figure  5A.  Typical  mode  II  deiamination  load  deflection  record  of 
T6CIA5  F 1 85  using  the  end-loaded  split  laminate  test. 


Table  3.-  End-loaded  Split  Laminate  (ELS)  Results. 


Material 

Specimen 

GIIc  (J/^2) 

I  II 

AS4/3502 

1 

630 

595 

613 

560 

613 

543 

648 

543 

2 

595 

525 

595 

525 

508 

508 

Av. 

595 

545 

★ 

T6C145/F1S5 

1 

1645 

2153 

2030 

2030 

2275 

2415 

2275 

2363 

2188 

1873 

2275 

2118 

2 

1960 

2643 

2223 

2240 

2415 

2083 

2590 

2170 

2713 

2765 

2800 

2328 

Av. 

2275 

2260 

*  Nonlinearity  is  too  great  to  satisfy  linear  beam  theory  assumptions  implicit 
in  the  approach 
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surprising,  since  an  a/L  ratio  of  at  least  0.55  was 
introduced  to  all  ELS  specimens  during  the  mode  I  precracking 
done  prior  to  mode  II  testing;  according  to  stability 
analysis,  crack  growth  should  have  been  stable"^.  This  is 
probably  the  result  of  the  additional  energy  storage  in  the 
fixtures.  Also,  the  linear  beam  theory  used  in  the  stability 
analysis  did  not  take  into  account  the  nonlinearity  observed, 

especially  at  displacements  greater  than  50.8  mm.  An  average 

2  2 

GIIc  of  595  J/m  (method  I)  and  543  J/m  (method  II)  were 
found  for  AS4/3502.  These  results  are  within  10  percent  of 
each  other.  However,  considerable  scatter  was  observed  in 
the  GIIC  data  from  each  specimen  (Table  3).  This  scatter  is 
believed  to  be  the  result  of  the  significant  amount  of  error 
involved  in  measuring  crack  locations  which  is  done  by 
monitoring  the  crack  as  it  propagates.  In  contrast,  this 
scatter  was  very  small  in  the  ENF  results.  A  typical 
load-deflection  curve  for  T6C145/F185  is  seen  in  Fig.  54. 
Stable  crack  growth  is  observed  for  all  load/unload  loops. 
This  is  evidenced  by  the  smooth  decrease  in  loading  rate 
followed  by  a  gradual  decrease  in  load.  It  should  be  noted 
that  it  is  not  known  when  crack  growth  started  during  this 
process,  which  was  common  to  both  types  of  mode  II  tests  for 
this  material.  At  least  some  of  the  observed  nonlinear 
behavior  is  due  to  inelasticity  since  the  curves  fail  to 
return  to  the  origin  as  the  specimen  is  unloaded.  Nonlinear 
behavior  due  to  large  displacements  and  rotations  is  also 
observed  for  deflections  greater  than  50.8  mm. 
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The  average  mode  II  energy  release  rate  values  for 

T6C145/F185  were  found  to  be  2275  J/m^  (method  I)  and  2260 

j/m^  (method  II).  However,  a  systematic  or  monpton i c 

increase  in  Gllc  for  method  I  can  be  observed  throughout  the 

test.  This  is  the  result  of  the  nonlinearity  in  the 

load/deflection  curve  seen  in  Fig.  54.  Although  this 

nonlinearity  is  less  severe  than  the  nonlinearity  observed 

for  this  material  tested  using  the  ENF  test,  it  is  sufficient 

to  invalidate  the  data  analysis  done  using  linear  beam  theory 

(method  I).  The  scatter  of  the  results  for  method  II  (area 

method)  is  more  random  and  is  associated  with  the 

experimental  difficulty  in  obtaining  a  precise  measure  of  the 

change  in  crack  length.  However,  these  values  are  an 

overestimation  of  GTT  because  the  curves  did  not  return  to 

lie 

the  origin  as  the  specimens  were  unloaded.  Thus,  some  far 
field  energy  dissipation  is  unavoidably  included  in  Gjj  . 

A  comparison  of  the  results  predicted  using  both  test 
configurations  indicate  that  for  AS4/3502  (a  brittle 
composite)  similar  and  valid  results  were  obtained  for  the 
ENF  and  ELS  tests.  In  contrast,  for  T6C145/F185  neither  the 
ENF  nor  the  ELS  tests  give  valid  results,  since  the 
assumptions  associated  with  each  test  method  were  violated. 
Therefore,  in  order  to  properly  characterize  the  Cracture 
toughness  of  this  material  and  tougher  graphite/epoxy 
composites  in  general,  an  analysis  which  allows  for  material 
and  geometric  nonlinearity  such  as  a  J-integral  analysis  is 
needed.  It  is  worth  noting  that  similar  results  were 


obtained  for  the  ENF  and  ELS  when  data  was  analysed  using  the 
area  method.  However,  both  results  should  be  considered 
upper  bounds  of  the  mode  II  critical  energy  release  rate  of 
the  material.  In  the  next  subsection,  an  analysis  for  J 
will  be  presented. 

4. 4.2.2  J-Integral  Approach  for  Mode  II  Delamination  Fracture 
Toughness  Evaluation 

Most  of  the  current  analyses  to  measure  delamination 
resistance  in  composites  are  elastic  analyses  based  on  linear 
fracture  mechanics  (i.e.  energy  release  rate).  However, 
these  analyses  may  become  inappropriate  when  the  inelastic 
material  behavior  observed  in  the  delamination  process  is 
sufficient  to  give  nonlinear  behavior  in  the  measured 
load-displacement  curves  and/or  permanent  deformation  on 
unloading.  This  inelastic  behavior  may  be  due  to 

microcracking  or  extensive  resin  deformation.  In  fact, 
Corleto  and  Bradley  have  noted  that  under  mode  II  loading, 
composite  systems  with  tough  matrices  show  considerable 
non-linear  behavior  due  to  extensive  microcracking  and 
yielding  of  the  resin  prior  to  macroscopic  delamination  crack 
extension^  . 

Successful  fracture  characterization  of  metals  with 

inelastic  behavior  (i.e.  small  scale  yielding)  has  been  done 

3  8 

by  means  of  Rice's  path- independent  J  integral  .  The  use  of 

this  parameter  as  a  fracture  criterion  has  been  based  on  a 

39_40 

consideration  of  the  Hutchinson-Ri ce-Rosengren  ( HRR ) 


crack  tip  model  where  the  stress  and  strain  fields  ahead  of  a 

4  1 

crack  may  be  characterized  in  intensity  by  J.  McClintock 
has  further  demonstrated  that  for  deformation  plasticity 
theory  a  singularity  in  r  (the  radial  distance  from  the  crack 
tip)  exists  which  is  strongly  related  to  the  J  integral. 
However,  the  path  independence  of  this  integral  may  not  hold 
true  when  variable  damage  develops  as  a  result  of  unloading. 
Therefore,  the  J  integral  is  limited  to  linear,  nonlinear 
elastic,  and  elastic-plastic  behavior  under  monotonic 
loading.  Alternatively,  the  critical  value  of  J  may  be 
thought  of  as  giving  the  energy  release  rate  required  for 
crack  growth  for  nonlinear  elastic  materials. 

In  the  case  of  fiber  reinforced  composite  materials, 
the  use  of  the  J  integral  as  a  failure  criterion  has  been 
very  limited  and  confined  to  experimental  approaches  (i.e., 

indirectly  measuring  J  from  load-displacement  data  via  its 

42-45  42 

energy  interpretation  ).  However,  Agarwal  et.al  ,  have 

shown  that  for  short  fiber  composites,  the  stress-strain 

behavior  of  the  composite  may  be  closely  approximated  by  the 

Rambe r g-Osgood  relation  (i.e.,  =  a^(e^)n)x,  justifying 

the  use  of  the  J  integral  as  a  fracture  criterion  for  short 

fiber  reinforced  composite  materials.  Furthermore,  they 

argued  that  since  the  local  unloading  of  the  matrix  resulting 

from  microcracking  and  debonding  at  the  interface  does  not 

necessarily  also  occur  in  the  fibers,  this  observed  crack  tip 

damage  does  not  lead  to  material  unloading  to  the  same  extent 

as  a  crack  extension  does  in  metals.  Therefore,  although  the 
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requirement  of  no  unloading  is  not  completely  met,  its 

influence  will  be  quite  small  for  some  appl i ca t i ons ,  and 

thus,  may  be  ignored  for  some  applications. 

In  a  recent  investigation,  an  additional  approach  for 

cha racte r i z ing  fracture  of  inelastic  composite  materials  has 

46 

been  made  by  Schapery  .  He  uses  strain-energy  like 

potentials  to  characterize  material  response  (i.e.,  a ^  = 

3W/3e^j).  Theoretical  support  for  the  existence  of  these 

potentials  is  given  by  the  thermodynamics  of  reversible 

processes.  The  potentials  are  free  energy  and  internal 

energy  for  isothermal  and  adiabatic  processes,  respectively. 

Schapery  argues  that  if  a  potential  of  this  nature  exists,  it 

3  8 

is  often  possible  to  use  Rice's  J  integral  theory  to 
simplify  fracture  analysis. 

In  the  present  analysis,  an  approach  similar  to 
Schapery's  is  proposed  to  evaluate  the  mode  II  fracture 
toughness  of  composite  laminates  with  inelastic  behavior.  J 
is  defined  for  the  end-loaded  split  laminate  test  (mode  II) 
specimen  and  the  implications  of  the  analysis  reviewed. 
Also,  some  limited  assessment  of  the  proposed  method  has  been 
done.  A  more  complete  evaluation  of  the  utility  of  Jjj  to 
evaluate  mode  II  delamination  fracture  is  being  undertaken  in 
the  new  contract. 

4. 4. 2. 2.1  Theoretical  Model 

Figure  5  shows  the  unsymmetric  end-loaded  split 
laminate  test  used  to  determine  the  critical  energy  release 
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rate  G„  for  the  mode  II  delamination  of  fiber-reinforced 
lie 

composite  materials.  The  test  consists  of  loading  an  split 

laminate  at  the  cracked  end,  with  the  opposite  end  fixed  to 

achieve  a  mode  II  stress  state  at  the  crack  tip.  As  has  been 

previously  noted,  inelastic  behavior  has  been  observed  for 

29 

composite  systems  with  tough  matrices  .  Further,  inelastic 

behavior  is  expected  if  a  significant  volume  fraction  of  off 

angle  fibers  are  present  in  the  layup.  This  inelastic 

behavior  then  invalidates  the  standard  elastic  methods  used 

to  measure  fracture  toughness.  In  order  to  develop  an 

analysis  where  inelastic  behavior  is  accounted  for,  consider 

Rice's  J  integral,  as  extended  to  crack  growth  in  inelastic 

4  7 

media  with  large  deformation  by  Schapery 

J  =  C1  [wQdx2  -  T^  ( 9ui/3x1  )  dL  (10) 

where  wq  is  the  work  potential  density,  T^  the  tractions 
along  ,  and  u^  the  displacements.  A  rigorous  proof  and 
support  for  the  use  of  work  potentials  to  characterize  the 
behavior  of  materials  and  structures  with  damage  is  given 
elsewhere.  Figure  56  shows  the  counterclockwise  integration 
path  C^.  The  contribution  to  J  from  this  contour  comes  from 
the  left  and  right  vertical  segments  only,  since  for  the 
outer  horizontal  segments,  dx2  ■  T^  =  T2  =■  0  which  also  holds 
true  for  the  inner  vertical  segments  (crack  surfaces)  if 
friction  is  assumed  to  be  negligible  between  the  crack 
surfaces.  Further,  assuming  small  strains  and  that  the  left 
vertical  segment  experiences  small  rotations  and  right 
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vertical  segment  is  evaluated  far  enough  away  from  the  crack 
tip  and  fixed  end  as  to  avoid  nonuniformity  in  the  stress 
field,  using  small  rotation  theory 


2  M  du2  1  M.  du2 

J  -  [  fk  dM  -  P  -  ]  -  -  [  J  k  dM  -  2P -  ]  (11) 

11  B  A3  dx .  LB  ~0  dx1  R 

where  M  and  du^/dx^  are  the  moment  and  slope,  respectively. 

The  subscript  L  and  R  indicate  evaluation  of  the  integrals 

corresponding  to  the  left  and  right  vertical  segments, 

respectively.  B  is  the  width  of  the  specimen.  If  k  « 
2 

d2U2/dx^  is  used  to  obtain  the  slope,  Equation  11  reduces  to 


J 


1 1 


1 

B 


[2 


Ma 

J  k  dM]  - 
J0  L 


(12) 


where  Mg  corresponds  to  the  moment  at  the  crack  tip.  (See 
Appendix  I  for  a  detailed  derivation  of  this  equation). 
Equation  12  indicates  that  the  result  is  independent  of  the 
location  of  the  segments  and  both  are  equal  to  that  evaluated 
at  the  crack  tip.  This  path  independence  will  be  met  as  long 
as  the  crack  propagates  at  a  constant  maximum  moment  M  .  For 
a  long  laminate  that  is  initially  homogenous  in  the  x^ 
direction,  constant  M  indicates  that  the  stress  and  strain 

Cl 

state  in  the  neighborhood  of  the  crack  tip  is  independent  of 
crack  length  (self  similar  growth  condi* ions).  Under  this 
condition,  the  energy  release  rate  can  be  shown  to  be 
equal  to  J11  calculated  using  Eq.  12.  Further,  since  small 


131 


rotations  are  only  required  from  the  fixed  end  up  to  the 
neighborhood  of  the  crack  tip,  Eq .  12  will  also  be  valid  for 

large  rotations  at  the  loaded  end.  See  Ref.  46  for  further 
explanation  of  the  ability  of  this  type  of  analysis  to 
characterize  fracture. 

Equation  12  can  be  shown  to  reduce  to  G  as  given  by 
linear  beam  theory  by  substituting  k=M/EI,  and  noting  that 
M-M  -Pa  for  the  left  segment,  M  =  M  =  2Pa  for  the  right  segment, 

3  3 

48 

and  I  =1/8  .  This  gives, 

L  K 


j 


1 1 


-  n2  2 

6  Pa 


B  E  I, 


(13) 


where  I  is  the  moment  of  inertia  of  the  uncracked  segment  of 
K 

the  beam. 


4. 4. 2. 2. 2  Experimental  Methods 

Limited  experimental  verification  of  as  given  by 
Eq .  12  has  been  done  with  a  1.0  x  .121  in.  x  8.0  in.  and  a 
1.0  in.  x  .188  in.  x  6.0  in.  split  laminate  of  unidirectional 
AS4/3502  and  T6T145/F185  composites,  respectively.  The 
former  one  has  a  brittle  matrix  and  the  latter  one  a  tough 
matrix.  Typical  mode  II  delamination  tests  have  been 
performed  on  both  composite  systems,  carefully  recording  the 
moment  arm  (distance  from  load  line  to  crack  tip)  and  load  at 
onset  of  crack  growth.  Moment-curvature  relationships  (e.g., 
Fig.  57)  for  these  composites  have  been  obtained  with  beams 
having  one-half  (cracked  portion)  and  full  lay-up  (uncracked 
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portion)  of  the  split  laminates  tested  in  mode  II 

delamination.  Finally,  Jjj  calculations  where  obtained  with 

1 

J  =  -  [  2  (shaded  area  L)  -  (shaded  area  R)]  (14) 

1 1  B 

where  shaded  area  L  and  shaded  area  R  are  shown  in  Fig.  57. 

4. 4. 2. 2. 3  RESULTS  AND  DISCUSSION 

Figure  58  shows  the  load/deflection  curves  obtained 

for  the  delamination  tests.  Stable  as  well  as  unstable  crack 

growth  occurred  on  both  systems.  Figures  59  and  60 

correspond  to  the  moment-curvature  relationships.  For  both 

systems  these  relationships  were  linear.  Figures  61  and  62 

show  the  J  _  _  and  GTT  values  measured,  plotted  as  a  function 
lie  lie 

of  crack  length  for  AS4/3502  and  T6T145/F185,  respectively. 
Also,  Tables  4  and  5  show  the  overall  results  obtained  during 
the  tests.  An  average  J  c  of  2.5  lb/in  and  a  Gllc  of  2.7 
lb/in  for  AS4/3502  and  Jjlc  of  12.0  lb/in  and  GIIc.  In  the 
case  of  AS4/3502  the  difference  between  J  and  G  amounts  to  no 
more  than  10  percent.  This  small  discrepancy  can  be 
explained  by  carefully  observing  the  load/ deflection  record 
in  Figure  58a.  Notice  how  the  curves  consistently  have  a 
steeper  slope  as  the  load  increases  indicating  a  stiffening 
effect.  Obviously,  ^  which  is  based  on  linear  beam  theory 
will  measure  higher  values  of  the  energy  release  rate.  This 
stiffening  effect  is  not  well  understood  yet  but  in  this 
case,  it  seems  that  it  is  not  due  to  large  rotations  and 
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(a)  Moment-curvature  for  loading  of  left  segment 
of  split  laminate  beam  (cracked  ) 


(b)  Moment-curvature  for  loading  of  right  segment 
of  split  laminate  beam  (uncracked) 


Figure  57.  Moment-curvature  test  needed  to  experimentally  evaluate  Jjj- 
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S4/3502 


record  of' 
sites . 


5/FI85 


Figure  58.  (Continued). 


AS4/3502  (Thickness  =  h/2) 


Figure  59a.  Moment-curvature  relationship  correspond ing  to  the  cracked 
portion  of  unidi  rsct  io3nl  AS4>/3502  composit0. 
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MOMENT  (IN-LB) 


250 


AS4/3502  (Thickness  =  h) 


Figure  59b.  Moment -curvature  relationship  corresponding  to  the  uncracked 
portion  of  unidirectional  AJ9/3502  composite. 
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MOMENT  (IN-LB) 


T6T1 45/FI  85  (Thickness  =  h/2) 


Figure  ooa .  Moment -curvature  relationsmp  corresponding  to  the  cracked 
portion  of  unidirectional  7611145  F  ’  3  5  composite. 
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MOMENT  (IN-LB) 


CURVATURE  (1/IN) 


Figure  cut.  Moment-curvature  relationship  corresponding  to  the  uncracked 
portion  of  unidirectional  T6T1il5.'F  185  composite. 


AS4/3502  (Unidirectional) 


Figure  61.  Mode  II  fracture  toughness  as  a  function  of  crack  length  of 
un  1 J  i  reet  ions  1  AS4/3502  composite. 
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deflections  since  the  a/L  ratios  for  this  test  are  within  the 
linear  regime.  Instead,  friction  between  the  cracked 
surfaces  can  be  a  factor  contributing  to  the  problem, 
especially  as  the  crack  grows.  If  this  is  the  case,  the 
proposed  j  will  not  take  into  cons i de r a t i on  these  friction 
effects,  and  thus,  will  give  an  artifically  high  indication 
of  the  energy  required  per  unit  area  of  crack  growth. 

For  T6T145/F185,  Gjlc  i s  at  least  20  to  25%  higher 
than  J  .  Observing  the  load /de f 1 ec t i on  record  (Fig.  58b) 
for  this  system  during  de 1  ami na t i on ,  it  can  be  seen  that  the 
curves  failed  to  return  to  the  origin  upon  unloading, 
indicating  inelastic  deformation  during  the  test.  It  is  very 
clear  that  Gllc  will  yield  a  overestimate  of  the  critical 
energy  release  rate  because  some  energy  was  dissipated  in  the 
far  field,  but  was  included  in  the  energy  associated  with 
creating  an  new  surface.  Corleto  and  Bradley  have  noted  that 
the  inelasticity  for  these  composite  systems  with  tough 
matrices  is  the  result  of  considerable  matrix  deformation  and 
yielding  ahead  of  the  crack  tip  (long  damage  process  zone)^. 
In  contrast,  jijc  appears  to  successfully  measure  the 
fracture  toughness  of  this  material  as  evidenced  by  the  lower 
JTT  values  measured.  Further,  the  small  scatter  and 
consistency  found  in  the  measured  values  indicate  a 
constant  crack  tip  moment  when  the  crack  was  growing. 
Therefore,  the  path  independence  and  self  similar  growth 
conditions  assumed  by  the  analysis  were  not  violated. 

Finally,  even  these  this  limited  results  are 


Table  4.  Mode  II  Fracture  Toughness  of  Unidirectional  AS4/3502  Composite. 


Crack  Length 
( in) 

Moment  Arm 
( in) 

pc 

Mc 

(in-lb) 

( in-^b ) 

^fb/in) 

JIIc 

6.53 

6.4 

9.9 

31.68 

63.36 

2.53 

2.38 

6.6 

6.45 

9.95 

32.09 

64.18 

2.61 

2.44 

6.72 

6.5 

9.9 

32.18 

64.35 

2.68 

2.46 

6.79 

6.6 

9.86 

32.54 

65.08 

2.7 

2.51 

6.88 

6.7 

9.69 

32.46 

64.92 

2.68 

2.50 

7.05 

6.8 

9.5 

32.1 

64.2 

2.71 

2.44 

7.17 

6.75 

9.57 

32.3 

64.6 

2.85 

2.47 

Avg. 

2.68 

2.46 

Pc  -  Critical  load. 

-  Crack  tip  moment  applied  to  uncracked  portion  of  split  laminate. 
Mc  -  Crack  tip  moment  applied  to  each  of  the  cracked  portions  of 
split  laminate 


Table  5.  Mode  II  Fracture  Toughness  of  Unidirectional  T6T145/F185  Composite. 


Crack  Length 
(  in) 

Moment  Arm 
(in) 

(1§) 

( in-Jb) 

( in-fb) 

c,V/i„>,,c 

3-3 

3.24 

81.68 

132.24 

264.48 

14.45 

11.66 

3.5 

3.44 

79.61 

136.91 

273.82 

15.45 

12.45 

5. 18 

5.05 

52.6 

132.81 

265.63 

14.78 

11.76 

5.33 

5.3 

50.6 

134.09 

268. 18 

14.31 

11.99 

5.55 

5.48 

49. 

134.26 

268.52 

14.35 

12.02 

Avg. 

14.7 

12. 
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T6T1 45/FI  85  (Unidirectional) 


Figure  62.  Mode  II  fracture  toughness  as  a  function  of  crack  length  of 
unidi reet  ional  T6T145/F185  composite. 
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encouraging,  a  complete  assessment  of  this  analysis  will 
require  testing  composites  with  different  extents  of  damage 
induced  by  layups  dominated  by  off  axis  fiber  orientations. 
Further,  the  effects  of  variable  ply  thickness  and  specimen 
widths  on  J  should  be  assessed. 


4.4.2. 3  MIXED  MODE  DELAMINATION 

Using  a  asymmetrically  loaded  split  laminate  (see  Fig. 
63),  one  can  achieve  a  variety  of  combinations  of  mode  I  and 
mode  II  loading.  The  total  energy  release  rate  is  given  by 
the  sum  of  the  mode  I  and  mode  II  values,  Gj  +  GIIc' 
calculated  using  the  relationships: 


GT  =  GI  +  GII 


(PL  ) ' 

SC 

BE  I 


3IPaLc^ 

4BEI 


(15) 


where  I  is  the  moment  of  inertia  of  one  leg  of  the  split 
portion  of  the  laminate,  P  =  (P  +P  )/2  and  P  -  ( P  -P  )/2. 
For  composites  made  from  brittle  resins,  the  GT  is  seen  to 
increase  with  increasing  percentage  of  mode  II  loading, 
especially  above  about  40%  mode  II.  By  contrast,  the  G^,  for 
a  composite  made  from  a  ductile  resin  system  is  found  to  have 
very  little  dependence  on  the  %GII  used  in  achieving 
delamiantion .  Typical  results  have  previously  been  presented 
in  Fig.  12  (section  4.1)  with  additional  results  provided  in 
Table  4.  The  reason  for  this  behavior  has  previously  been 
explained  in  detail  in  section  4.1. 
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Figure  63  Schematic  showing  how  asymmetric  loading  of 

split  laminate  can  introduce  a  mixed 
mode  I/mode  II  state  of  stress  at  the  crack 
tip. 
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4.5  FRACTORGRAPHY  AND  FAILURE  ANALYSIS 


There  has  been  considerable  confusion  over  whether  it 
is  possible  to  determine  crack  growth  direction  in  an 
unambigeous  way  from  the  direction  the  hackles  point.  The 
results  of  the  in-situ  observations  in  this  study  clearly 
indicate  that  hackles  can  point  either  in  the  direction  of 
crack  growth  or  opposite  to  the  direction  of  crack  growth,  as 
previously  seen  in  Fig.  9.  Furthermore,  the  development  of 
damage  ahead  of  the  crack  tip  gives  some  serious  question 
about  how  reliable  some  of  the  surface  features  may  be  in 
determining  crack  growth  direction.  Unless  all  of  the  damage 
ahead  of  the  crack  tip  develops  in  a  very  consistent  way, 
using  the  artifacts  that  result  from  this  damage  can  be 
misleading  in  f r ac tog r aphi c  interpretation.  A  great  deal 
more  work  needs  to  be  done  to  establish  that  such  a  regular 
pattern  is  present  in  the  development  of  the  damage  zone  that 
unambigous  information  concerning  the  crack  growth  direction 
can  be  inferred  form  fracture  surface  details.  More  details 
in  this  regard  are  contained  in  papers  5.4  and  5.6  section 
5.0  of  this  report. 
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and  Delamination  Fracture  Toughnesses  of 
Graphite/Epoxy  Composites,"  ASTM  Symposium  on 
Fractography  of  Modern  Engineering  Materials, 
Nashville,  November,  1985,  invi ted . 

"Micromechanisms  of  Fracture  in  Toughened 
Graphite/Epoxy  Laminates,"  ASTM  Symposium  on  Toughened 
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6.16 


6 . 17 

6 .18 

6 . 19 

6  .20 

6 . 21 


7.0 
7  .  1 


Composites,  Houston,  TX,  March,  1985. 

"Delamination  Fracture  Studies  of  Some  Toughened 
Graphite/Epoxy  Composites,  Including  Real-Time 
Fracture  Observations  in  Scanning  Electron 
Microscope,"  ASTM  Symposium  on  Toughened 
Composites,  Houston,  TX,  March,  1985. 

"Delamination  Fracture  Behavior  of 

Graphite/Thermoplastic  Composite  Systems,"  invited , 
NASA  Langley,  Hampton,  VA.,  February,  1985. 

"Mi c romechani cs  Modelling  of  Delamination  Fracture  in 
Graphite/Epoxy  Composites,"  Gordon  Conference  on 
Composite  Materials,  Santa  Barbara,  CA.,  Jan.,  1985. 
"In-Situ  Studies  in  Delamination  of  Graphite/Epoxy: 
Recent  Results",  Air  Force  Composite  Materials  Review 
Dayton,  Ohio,  Oct.,  1984,  invited . 

"In-Situ  Studies  of  T6T145/F185  Composite  Materials," 
presented  at  the  6th  International  Fracture  Symposium, 
New  Delhi,  India,  Dec.,  1984. 

"In-Situ  SEM  Studies  of  Fracture  of  Plastics  and  Fiber 
Reinforced  Plastics".  Dow  chemical,  Research  Center, 
Freeport,  TX,  July  1984,  invi ted . 

PROFESSIONAL  PERSONNEL  ASSOCIATED  WITH  RESEARCH  EFFORT 
William  M.  Jordan,  Ph.D.,  awarded,  December,  1985, 
dissertation  title:  "The  Effects  of  Resin  Toughness 
on  the  Delamination  Fracture  E-'havior  of 
Graphite/Epoxy  Composite". 
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7.2  Carlos  Corleto,  M.S.,  awarded,  June,  1986,  thesis 
title:  "Mode  II  Delamination  Fracture  Toughness  of 
Unidirectional  Gr aphi te/Epoxy  Composites".  Currently 
working  on  Ph.D.  in  area  of  Mode  II  Delamination  with 
expected  graduation  date  of  December,  1988. 

7.3  Mike  Hibbs,  Ph.D.  candidate  with  expected  graduation 
date  of  Fall,  1987,  working  on  strain  field  mapping 
around  crack  tip. 

7.4  Douglas  Goetz,  Ph.D.  candidate  for  May,  1988,  is 
working  on  J-Integral  approach  for  Mode  I  delamination 
of  multiaxial  layups  (supported  half  by  Bradley  and 
half  by  Schapery). 
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APPENDIX  I 

Derivation  of  J  for  the  end-loadeu  split  laminate  tost 
Assumptions : 

(1)  material  homogeneity  with  respect  to  x^  axis. 

(2)  small  rotations  from  fixed  end  of  beam  to  the 
neighborhood  of  the  crack  tip.  Therefore,  the 
following  apply: 

u2,  vertical  displacement  of  the  center  line  of  beam 


U1  =  X1 


du2a 

3u^ 

d  2  u  2 

dx^ 

/ 

CM 

< 

II 

H 

X 

ro 

dx^ 

2  2 

k  -  d  u2/dx^  ,  curvature  ; 


,  slope 


(3)  friction  forces  at  crack  interface  are  negligible  (this 
is  accomplished  in  the  specimen  by  placing  a  teflon  insert 
between  the  crack  surfaces). 

Analysis : 

4  7 

Using  Rice’s  J  integral  as  extended  by  Schapery  , 

J  y[wodx2  -  T.  Oui/3x1)]dL  (Al) 


Jjj  for  the  specimen  configuration  of  Fig.  1.1  will  be  given 
by 

J  j  j  =  +  ^3  +  ^4  +  J  5  +  J  5  +  ^7  ( A  2  ) 


where  J.  correspond  to  the  contribution  of  JI3.  coming  from 
the  counterclockwise  contour  shown  in  Fig.  1.2.  =»  = 

J5 

»  0  since  dx2  =  Ti  =  0  for  these  cases  (see  Fig.  1.2). 
Therefore,  the  only  contribution  to  J  comes  from  J2,  and 
J^.  The  two  former  ones  are  equal  due  to  symmery  and  will  be 
evaluated  close  enough  to  the  crack  tip  so  that  small 
rotation  beam  theory  applies  (assumption  2).  will  be 
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evaluated  far  enough  from  the  crack  tip  and  fixed  end  to 
avoid  nonuniformities  in  the  stress  field.  Let's  consider  J 

l 

first.  Since  TV  =  ai^nj  anc^  noting  that  the  limits  of 
integration  for  path  4  are  -h/2  and  h/2  for  the  lower  and 
upper  limits,  respectively, 


J4  = 


rh/2 

r h/2  3u^ 

9u. 

1  w 
/  0 

dx2  1  (ffllnl  +  T2 

— £  )  dL] 

(A3  ) 

I 

J  axl 

3x. 

-h/2 

-h/2  1 

1 

Using  the  expressions  stated  in  assumption  2,  and  further 
assuming  that 


w  = 
o 


*1 lde 1 1 


neglecting  shear) 


and  since  dL  =  dx2  and  the  unit  outward  normal  n^  =  1,  Eq.  A3 
becomes 


.h/2  £. 


J4  “ 


X2’/ 


-h/2  0 

Now  by  means  of  the  following 


fh/2 

-h/2 


du0  r h/2 


*n  eii  dx2  ' 


2  fn/ 


T2dx2 


-h/2  (A3a) 


e  ^  =  -  k  x  2  ;  d  e  ^  =  -  x 2  dk  ; 


f  T  -7  - 

-h/2  B 


and  performing  a  change  of  variables, 
-k  .h/2  h/2 

J  &  =|  (  — J  *llx2dx2  ^  ^ ^ ~ k  f  (  —  o  1  1 

r\  •/  % 


/  /•"'  ~  2P  du2 

°1 lx2dx2 * dk-k/  <-*llx2dx2^  +  —  -  (A. 3b) 

1  /  B  dx. 

-h/2  -h/2  1 


Further, 


r  h/2 

/  <-*ii 

J  -h/2 


x2dx2  > 


thus , 


1  /* 

—  M 

B  J  0 


k  2P  du- 

dk  -  —  M  +  —  — : 

B  B  dx. 
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T 


T 


Integrating  by  parts  ( u  =  M  du  =  dM  dv  =  dk  v  =  k ) 

1  M  k  2P  du. 


J  =  —  ( M  k 

B 


1 

B 


/ 
*  n 


2pdu2  r 
dxx  J 


M 


dM) 


dM) 


B 


B  dx. 


( A3c  ) 


To  evaluate  or  ^  the  same  procedure  as  that  followed  for 
J4  was  performed.  For  these  cases  one  must  notice  that  dL  » 
-dx^  and  n^  =  -1  and  thus, 


f 


lldeUdx2  + 


/ 


’ll 


3u . 


3x . 


dx . 


h/2  0 

which  can  be  reduced  to 


h/2 


1  ,-M  P  dup 

j  -  j  -  -  /  k  dM - ^ 

z  0  B  J  B  dx . 

o  1 

Addidng  J2,  and 


du^  P 

- -  -  (A4) 

dx^  B 


( A4a  ) 


2  fM  du-,  1  .M  dup 

J  =  -  [  /  k  dM  -  P  — -  1  -  -  l  k  dM  -  2P  — -  ]  (A2a) 

1  B  JO  dx x  L  B  J  0  dx  ^  R 

where  the  subscripts  L  and  R  indicate  evaluation  of  the 
quantities  inside  the  brackets  corresponding  to  the  left 
(paths  2  and  6)  and  right  (path  4)  vertical  segments  of  the 
beam . 

in  order  to  establish  the  path  independence  of  J  ,  it 

will  be  evaluated  at  a  distance  x^  from  loaded  end. 

Left  segment : 

2  2 

since  d  U2/dx^  =  k,  the  slope  of  any  segment  in  the  cracked 
part  of  the  beam  at  a  distance  x^  from  the  loaded  end  is 
given  by 


Jl. 
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uu 

2 

1 

/2Pa 

i  k 

dM]  +  -  i 

rPx 

dx  ^ 

2P 

J2PL 

R  P  o 

Pa 

k  dM] 


Substituting  this  expression  into  Eq.  A4a  and  evaluating  the 
integral  for  x  ^  , 


1  f 

=  -  {  12 

B  J  ( 


•Pa  2Pa 

k  dM  ]  -  (  /  k  dM  ]  } 

0  L  {  0  R 


where  JL  refers  to  the  contribution  to  j  coming  from  the 
left  segment  of  the  beam  (cracked  part).  As  it  can  be  seen, 
this  expression  is  independent  of  path  (i.e.  independent  of 


Right  segment : 

The  slope  of  any  segment  in  the  uncracked  part  of  the  beam  at 
a  distance  x^  from  the  loaded  end  is  given  by 


dx.  2  P  J 


k  dM  ]  R 


Substituting  the  above  expression  into  and  evaluating  the 
integral  for  x^  , 

1  f2Pxl  r  ®  1  /  a 


-  { 
B 


r2Px, 

</  ‘ 
J  9  PT. 


kdM  + 


■-  M 


-  [J  kdM] 

B  JO  R 


where  JR  refers  to  contribution  to  coming  from  right 

segment  of  beam. 

Finally,  by  adding  Eq2.  A3d  and  A5,  we  obtain 


1  r  Pa 

-  {  2  (  kdM] 

B  f  0  L 


-  f 


-2  Pa 

kdM] 

0  R 


which  is  Equation  12  as  shown  in  Theoretical  Model  section, 
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Mike  F.  Hibbs,  '  Ming  Kwan  Tse, 1  and  Walter  L.  Bradley 


Interlaminar  Fracture  Toughness  and 
Real-Time  Fracture  Mechanism  of  Some 
Toughened  Graphite/Epoxy  Composites 

Authorized  Reprint  from  Special  Technical  Publication  937, 1987 

Copyright  American  Society  for  Testing  and  Materials,  1916  Race  Street  Philadelphia,  PA  19103 


REFERENCE:  Hibbs,  M  F .  Tse,  M  K..  and  Bradley.  W  L  ,  “Interlaminar  Fracture 
Toughness  and  Real-Time  Fracture  Mechanism  of  Some  Toughened  Graphite/Epoxy 
Composites,"  Toughened  Composites.  ASTM  STP  937.  Norman  1  Johnsion,  Ed  .  Ameri¬ 
can  Society  for  Testing  and  Materials,  Philadelphia.  1987,  pp  115-130 

ABSTRACT:  Five  graphite/epoxy  composites  containing  toughened  epoxies  prepared  ai 
Dow  Chemical  and  AS-4  graphite  fibers  from  Hercules  have  been  studied  in  Mode  I,  mixed 
mode,  and  Mode  Q  to  determine  their  dclaminanon  fracture  toughnesses  G,  and  the'control- 
ling  micromechanism  of  fracture.  The  G,  values  were  determined  using  spin  laminate 
specimens.  The  dclaminanon  fracture  process  was  observed  in  real  time  in  the  scanning 
electron  microscope  To  increase  de lamination  fracture  toughness,  both  an  improved  inter¬ 
face  as  well  as  higher  resin  toughness  were  found  to  be  required.  Increasing  Mode  II 
loading,  particularly  of  the  more  bottle  systems,  gives  a  significantly  greater  resistance 
to  crack  propagation  as  measured  by  the  total  energy  release  rate  required  to  propagate 
the  crack. 

KEY  WORDS:  composite  material,  delinunauon,  mixed  mode.  Mode  I.  graphiteicpoxy. 
toughened,  in  situ  fracture 


Nomenclature 

B  Specimen  thickness 
C  Compliance  of  test  coupon 
E  Elastic  modulus  in  the  fiber  direction 
G(  Energy  release  rate  for  Mode  1  loading  (opening  mode) 

G a  Energy  release  rate  for  Mode  U  loading  (in  plane  shear) 

G-„  Critical  energy  release  rate  for  stable  crack  growth  for  Mode  !  loading 
Gxsc  Critical  energy  release  rate  for  stable  crack  growth  for  Mode  1J  loading 
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Engineering  Research  Center,  College  Station ,  TX  77843. 

'Senior  research  engineer.  The  Dow  Chemical  Company.  Resms  Research,  B-1215  Bldg.. 
Freeport,  TX  77541 . 
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William  M.  Jordan  '  and  Walter  L.  Bradley 


Micromechanisms  of  Fracture  in 
Toughened  Graphite- Epoxy  Laminates 

Authorized  Reprint  from  Special  Technical  Publication  937,  1987 

Copyright  American  Society  for  Testing  and  Materials,  1916  Race  Street  Philadelphia  PA  19103 


REFERENCE:  Jordan,  W.  M.  and  Bradley,  W  L.,  “Micromechanlsms  of  Fracture  in 
Toughened  Graphite-Epoxy  Laminates,”  Toughened  Composites.  ASTM  STP  937. 
Norman  J.  Johnston.  Ed  ,  American  Society  for  Testing  and  Matenals.  Philadelphia,  1987, 
pp.  95-114 

ABSTRACT :  The  combination  of  low  cross-link  density  and  elastomer  additions  has  been 
seen  to  give  the  most  potent  toughening  for  neat  resins  A  relatively  small  increment  of  the 
additional  neat  resin  fracture  toughening  above  800  J/m;  is  actually  reflected  in  the  delami- 
nauon  fracture  toughness  of  a  composite  Mode  11  delamination  toughness  of  brittle  systems 
may  be  as  much  as  three  times  the  Mode  I  delamination  fracture  toughness,  while  a  ductile 
system  may  have  a  Mode  11  delamination  fracture  toughness  that  is  similar  to  the  Mode  I 
value.  The  energy  absorbed  per  unit  area  of  crack  extension  for  delamination  seems  to  be 
independent  of  ply  orientation  if  proper  accounting  of  the  near  and  far  field  energy  dissi¬ 
pation  is  made. 

KEY  WORDS:  composite  materials,  delamination.  Mode  l.  Mode  It,  mixed  mode,  tough¬ 
ened  epoxy 


Nomenclature 
B  Specimen 

E  Elastic  modulus  in  the  fiber  direction 
Gk  Critical  energy  release  rate  for  stable  crack  growth  for  Mode  I  loading 
Giic  Critical  energy  release  rate  for  stable  crack  growth  for  Mode  II  loading 
Groce  Total  critical  energy  release  rate  for  mixed  mode  loading 
/  Moment  of  inertia  for  cracked  portion  of  split  laminate 
Lc  Crack  length  in  split  laminate 
P,  Asymmetric  load  component  (P.  +  P,)/2  (see  Fig.  1) 

P,  Symmetric  load  component  (P,  -  PL)/2  (see  Fig,  1) 

P '«  Load  applied  to  upper  half  of  split  laminate  (sec  Fig.  1) 


This  work  performed  at  Texxs  A&M  University  was  supported  by  the  Air  Force  Office  of  Scientific 
Research  with  Major  David  Glasgow  as  project  monitor. 

'Research  assistant  and  professor  of  mechanical  engineering,  Mechanical  Engineering  Department. 
Texas  A  AM  University,  College  Station,  TX  77843. 
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Mode  II  delamination  fracture  toughness  that  is  similar  to  the  Mode  I  value.  The 
energy  absorbed  per  unit  area  of  crack  extension  for  deiamination  seems  to  be 
independent  of  ply  orientation  if  proper  accounting  of  the  near  and  far  field 
energy  dissipation  is  made. 
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A  Comparison  of  the  Crack  Tip  Damage  Zone  for  Fracture 
of  Hexcel  F185  Neat  Resin  and  T6T145/F185  Composite 

E.  A.  CHAKACHERY  and  W.  L.  BRADLEY 

Texas  A&M  University 
College  Station,  Texas  77843 

Hexcel  F185  neat  resin  and  T6T145/F185  graphite  fi¬ 
ber-reinforced  composite  were  subjected  to  Mode  1  loading 
in  the  compact  tension  (CT)  geometry  (fibers  parallel  to  the 
crack)  and  the  energy  per  unit  area  of  crack  extension.  J,c. 
determined  to  be  8100  and  1600  J/m2  respectively.  In-situ 
fracture  studies  using  scanning  electron  microscopy  on  a 
CT-type  specimen  of  FI 85  showed  extensive  microcrack¬ 
ing  in  a  damage  zone  ahead  of  the  crack  tip,  which  was 
similar  to  the  microcracking  observed  in  the  whitened  area 
ahead  of  the  crack  tip  in  the  macroscopic  CT  specimens. 

A  simple  calculation  using  a  rule  of  mixtures  approach 
suggests  that  the  diminished  size  of  the  damage  zone  and 
the  presence  of  rigid  fibers  in  the  damage  zone  in  the 
composite  are  not  a  sufficient  explanation  for  the  signifi¬ 
cantly  lower  delamination  toughness  of  the  composite 
compared  to  the  neat  resin.  From  this  it  may  be  inferred 
that  the  strain  to  failure  locally  in  the  damage  zone  ahead 
of  the  crack  in  the  composite  may  also  be  lower  than  that 
which  can  be  tolerated  in  the  neat  resin.  Evidence  for  this 
idea  comes  from  the  observation  that  microcrack  coales¬ 
cence  seems  to  occur  preferentially  at  the  fiber/resin  in¬ 
terface. 


INTRODUCTION 

Delamination  in  graphite  fiber-reinforced 
resin  composites  is  well  known  to  limit  the 
extensive  use  of  these  materials  in  some  struc¬ 
tural  applications.  Much  effort  has  been  di¬ 
rected  in  the  past  decade  towards  improving  the 
delamination  fracture  toughness  of  graphite/' 
epoxy  composites  (1  to  10).  Since  the  delami¬ 
nation  crack  propagates  through  the  interlam¬ 
inar  resin  rich  region,  the  emphasis  has  been 
on  obtaining  epoxy  resins  with  improved  frac¬ 
ture  toughness.  The  addition  of  elastomeric 
modifiers  was  shown  (1  to  6)  to  dramatically 
increase  the  fracture  toughness,  Gfc.  of  the  ma¬ 
trix  resin,  but  resulted  in  a  quite  modest  in¬ 
crease  in  the  delamination  fracture  toughness 
of  the  composite  (7,  8).  This  was  shown  (8  to 
10)  to  be  primarily  due  to  the  reduced  volume 
ot  resin  deformation  in  the  damage  zone  ahead 
of  the  delamination  crack  tip  since  the  fibers 
act  as  rigid  fillers  in  the  ductile  matrix.  Evi¬ 
dence  in  support  of  this  conclusion  was  ob¬ 
tained  (9.  10)  from  the  observation  that  a  de¬ 
crease  In  the  fiber  volume  fraction  resulted  In 
an  Increase  in  the  delamlnation  fracture  tough¬ 


ness.  The  current  investigation  was  undertaken 
to  attempt  to  correlate  the  delamination  frac¬ 
ture  toughness  of  the  composite  with  the  frac¬ 
ture  toughness  of  the  neat  resin.  The  determi¬ 
nation  of  the  fracture  processes  occurring  in 
the  vicinity  of  the  crack  tip  in  both  the  compos¬ 
ite  and  the  neat  resin  was  a  necessary  first  step 
If  these  events  are  nominally  the  same,  then 
there  may  be  a  relationship  between  neat  resin 
toughness  and  delamlnation  toughness. 

To  study  the  relationship  between  neat  resin 
fracture  toughness  and  composite  delamination 
fracture  toughness,  in-situ  observations 
each  fracture  process  have  been  made  m  :t  • 
scanning  electron  microscope  (SEM)  Thes« 
servattons  have  been  correlated  with  « 
scopic  measurements  of  a  composite  made 
this  neat  resin.  A  simple  model  has  tv  . 
posed  to  relate  the  microscopic  ohserv.., 
the  fracture  processes  to  the  observe,- 
scopic  measurements  of  fracture  ;  no. 

The  resin  system  selected  tor  :>  . 
the  F185  formulation  commct.  • 
from  Hexcel  Corp  It  is  a  relative  , 

(9  percent  elongation  In  a 
with  a  reported  fracture  ■ 
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specimen  of  F185  resin  showing  considerable  non-line¬ 
arity. 


gral  analysis  according  to  ASTM  E8 13-81.  A 
well  defined  J-R  curve  was  obtained  from  which 
a  J/c  value  of  8100  J/m2  was  determined.  This 
is  considerably  higher  than  the  value  of  6000 
J/m2  reported  by  the  manufacturer  (11).  This 
difference  may  be  due  to  our  use  of  fatigue 
precracking  rather  than  razor  notching  or  due 
to  our  use  of  a  nonlinear  rather  than  linear 
analysis.  If  we  use  the  linear  elastic  analysis  in 
ASTM  E399-81.  we  obtain  a  K0  value  of  2700 
J/m2  from  the  5  percent  secant  offset  value  of 
PQ.  However.  Pmax  is  much  greater  than  1.1 
Pq.  and  the  specimen  thickness  is  less  than  that 
required  by  ASTM-E399  for  K0  to  be  equal  to 
K,c.  Thus,  the  linear  elastic  analysis  is  invalid 
for  our  specimen  size.  If  we  use  Pmax  and  the 
actual  crack  length  at  Pmax  in  a  LEFM  calcu¬ 
lation.  we  obtain  a  which  has  an  equivalent 

G/c  of  6300  J/m2  which  is  similar  to  the  6000 
J/m2  obtained  by  Bascom  and  coworkers  (2) 
using  a  similar  approach.  We  believe  that  the 
value  of  8100  J/m2  obtained  in  this  work  using 
fatigue  precracked  specimens  and  a  J-integral 
analysis  is  a  more  meaningful  measure  of  the 
fracture  toughness  of  the  FI 85  resin. 

A  J-integral  approach  was  also  used  to  ana¬ 
lyze  the  data  obtained  in  the  CT  specimens  of 
the  composite.  The  J,c  for  transverse  crack 
growth  thus  obtained  was  1600  J/m2  which 
compares  very  well  with  recent  results  by  Bas¬ 
com  (13).  The  macroscopic  test  results  may  be 
summarized  as  follows:  the  neat  resin  fracture 
toughness  was  8100  J/m2.  whereas  the  delam¬ 
ination  and  transverse  cracking  in  the  compos¬ 
ite  gave  toughness  values  of  1900  J/m2  and 
1600  J/m2  respectively.  Thus,  crack  growth  in 
the  composite  parallel  to  the  fiber  direction  is 
much  easier  than  crack  growth  in  the  neat 
resin. 

In-Situ  and  Post  Mortem  Fracture 
Observations  in  SEM 

Fracture  studies  conducted  in  the  SEM  on  the 
F 1 85  neat  resin  indicate  that  the  zone  ahead  of 
the  crack  tip  undergoes  extensive  microcrack¬ 
ing  [Fig.  2).  There  has  been  some  recent  contro¬ 
versy  as  to  whether  these  are  actually  micro¬ 


cracks  in  the  resin  or  microcracks  in  the  150A 
thick  gold-palladium  coating  applied  to  the 
specimens  to  minimize  charging.  To  try  to  verify 
that  these  are  actual  microcracks  in  the  mate¬ 
rial.  we  have  recently  polished  two  sides  of  an 
AS4/3501-6  graphite/epoxy  composite,  coated 
one  side  with  a  1 50A  thick  layer  of  gold-pallad¬ 
ium  in  the  usual  way,  delaminated  the  speci¬ 
men  in  the  SEM,  removed  the  specimen  with 
the  wedge  intact  to  avoid  viscoelastic  recovery 
and  closure  of  the  microcracks,  and  then  coated 
the  other  side  with  gold- palladium  in  the  usual 
way.  and  then  examined  it  in  the  SEM  with  the 
wedge  still  intact.  The  results  of  this  exercise 
are  presented  in  Figs.  3(a)  and  3(b).  The  micro¬ 
crack  morphology  on  the  side  coated  and  then 
deformed  is  essentially  indistinguishable  from 
the  microcrack  morphology  seen  on  the  side 
deformed  and  then  coated.  Thus,  we  believe 
that  the  Au/Pd  coating  is  not  responsible  for 
the  microcracking  seen  on  the  surface  of  our 
composite  specimens. 

The  microcracked  region  ahead  of  the  crack 
tip  of  the  FI 85  neat  resin  (Fig.  2)  extends  60  to 
70  /im  above  and  below  the  crack  tip  and  is 
shaped  somewhat  like  a  kidney  bean.  Crack 
propagation  under  fixed  grip  conditions  oc¬ 
curred  In  a  discontinuous  manner  via  a  time- 
dependent  coalescing  of  the  microcracks  ahead 
of  a  blunted  crack  tip.  With  each  successive 
extension,  crack  advance  of  120  to  150  Mm  was 
observed,  resulting  in  a  very  sharp  crack  tip 
with  relatively  few  microcracks  ahead  of  it.  The 
applied  crack  opening  displacement  was  then 
increased  essentially  instantaneously,  and  with 
this  increase,  the  microcrack  density  ahead  of 
the  crack  tip  would  gradually  increase  with 
time,  accompanied  by  crack  tip  blunting  at  the 
formerly  sharp  crack  tip.  Again,  a  critical  den¬ 
sity  of  microcracks  (or  crack  tip  strain)  is 
reached  and  crack  extension  occurs  by  micro¬ 
crack  coalescence  as  before.  Evidence  of  this 
process  is  seen  in  the  post-mortem  fracto- 
graphic  examination  as  relatively  flat  regions  of 
crack  extension,  separated  by  the  small  lips 
that  probably  arise  from  crack  tip  blunting, 
preceding  the  next  crack  extension  (Fig.  4).  This 
appearance  is  typical  over  the  whole  fracture 
surface,  indicating  that  the  phenomena  ob¬ 
served  at  the  surface  is  typical  of  the  bulk  frac¬ 
ture  behavior. 

If  crack  extension  does  indeed  occur  by  mi¬ 
crocrack  coalescence,  then  the  relatively  flat 
regions  between  lips  seen  in  Fig.  4(a)  and  (b) 
should  be  found  on  examination  at  higher  mag¬ 
nification  to  be  composed  of  many  small  facets 
with  pledges  between  them.  This  is  exactly  what 
is  observed  at  lO.OOOx  magnification,  as  shown 
in  Fig.  4(c).  It  is  worth  noting  that  even  at  this 
very  high  magnification,  the  individual  facets 
do  not  seem  to  have  fractured  by  a  brittle  cleav¬ 
age.  The  surface  is  relatively  textured,  with  cav- 
iation  presumably  at  the  rubber  particle  addi¬ 
tions.  which  suggests  that  the  microcracking  is 
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deformation  have  the  same  beneficial  effect  of 
redistributing  the  load  away  from  the  crack  tip 
and  lowering  the  local  stresses  at  the  crack  tip. 
Thus,  from  a  physical  point  of  view,  either  pro¬ 
cess  enhances  toughness. 

In-situ  delamination  studies  on  T6T145/' 
FI 85  composite  show  that  the  resin  micro- 
cracks  quite  extensively  in  the  vicinity  of  the 
crack  tip.  The  microcracking  extends  from  the 
Interlaminar  region  and  into  the  ply  about  50 
Mm  above  and  below  the  crack  (Fig.  5).  Here,  as 
in  the  resin,  the  crack  propagates  by  microcrack 
coalescence  and  In  the  same  discontinuous 
manner.  At  each  stage,  however,  the  crack  tip 
blunting  is  apparently  less  than  in  the  FI 85 
resin  and  the  discontinuous  extension  varies 
from  40  to  100  Mm.  The  coalescence  of  micro- 
cracks  occurs  predominantly  at  the  fiber  resin 
interface.  Thus  it  appears  that  the  presence  of 
the  fibers  prevents  the  development  of  the  full 
resin  toughness  from  being  realized  due  to  pre¬ 
mature  microcrack  coalescence  at  the  fiber 
resin  Interface. 

The  presence  of  the  rigid  fibers  also  restricts 
the  height  above  and  below  the  delamination 
plane  over  which  deformation  and  microcrack¬ 
ing  occurs  (100  fiin  in  the  composite,  compared 
to  130  ^m  in  the  resin).  However,  the  region  of 
microcracking  ahead  of  the  crack  tip  is  actually 
increased  from  25  Mm  in  the  resin  to  about  50 
m m  in  the  composite. 

The  fracture  surface  of  the  delaminated  spec¬ 
imens  showed  evidence  of  cavitation,  voiding, 
and  microcracking  (Fig.  6).  The  resin  region  of 
the  composite  delamination  fracture  surface 
shows  more  coarse  cavitation  than  in  the  neat 
FI 85  resin  fracture  surface,  where  voids  were 
relatively  fine  and  more  homogeneously  dis¬ 
persed.  Furthermore,  the  voids  seem  to  be  more 
dense  in  the  resin  adjacent  to  the  fibers  (Fig. 


Fig.  4.  Fracture  surface  of  FI  85  resin.  Crack  growth  Is 
from  left  to  right,  (a)  330x.  Away  from  the  free  surface. 
Arrows  indicate  a  lip  formed  by  crack  tip  blunting  after 
a  growth  sequence,  (b)  330x.  Just  adjacent  to  the  free 
surface,  (c)  lO.OOOx.  Detail  of  a  relatively  flat  region  In 
(al.  Cavitation  due  to  phase  separated  CTBN.  Facets  sep¬ 
arated  by  ledges  probably  correspond  to  microcrack  coa¬ 
lescence. 


tip  can  be  made  equally  well  In  either  case.  The 
possibility  also  exists  that  microshear  bands 
formed  perpendicular  to  the  surface  can  give 
rise  to  sharp  contours,  which,  when  viewed  on 
edge,  appear  as  microcracks  on  the  surface.  It 
should  also  be  noted  that  microcracking  and 


Fig.  5.  lOOOx.  In  situ  delamlnatlon  of  T6T145/F185  com¬ 
posite  showing  microcracking  In  the  damage  zone  ahead 
of  the  crack  tip.  Note  microcracking  is  more  dense  adja¬ 
cent  to  the  fibers.  Preferential  microcrack  coalescence 
near  the  resin  fiber  Interface  where  the  mlcroc  racks  arc 
Inclined  to  the  primary  crack. 
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m2  in  the  composite  delamination  fracture 
toughness.  The  first  two  of  these  three  factors 
can  be  quantified  based  on  actual  observations. 

If  one  assumes  that  microcracking  begins  for 
strains  above  a  threshold  strain,  then  the  extent 
of  the  microcracked  zone  can  be  used  to  quan¬ 
tify  the  magnitude  of  the  h  in  Eq  2. 

A  cross-section  of  the  composite  prepared  me- 
tallographically  to  reveal  the  microstructure 
(Fig.  8 )  may  be  used  to  determine  the  volume 
fraction  of  fibers  in  the  hypothetical  tensile  . 
specimen  ahead  of  the  crack  tip.  Since  the  vol¬ 
ume  fraction  is  quite  nonuniform.  the  micro- 
structure  was  divided  into  three  regions:  the 
resin  rich  region  between  plies  with  a  fiber 
volume  fraction  of  19  percent;  the  ply  region 
with  a  volume  fraction  of  fibers  of  76  percent 
and  a  transition  zone  with  a  volume  fraction  of 
approximately  33  percent.  The  relative  heights 
of  these  three  regions  are  shown  in  the  sche¬ 
matic  in  Fig.  9.  along  with  the  height  of  the 


Ftg.  8  600x.  Cross  section  of  T6T145/F 185  showing  in¬ 
terlaminar  resin  rich  region. 


L 


resin  2h  130  us 

- -  o 

-  composite  2h  100  un 

Fig.  9.  Schematic  showing  relative  height  of  damage 
zone  in  FI85  resin  and  T6T145/F185  composite. 


microcracked  zone  in  the  FI 85  resin.  Using  a 
simple  rule  of  mixtures  approach,  and  taking 
account  of  the  smaller  microcracked  (and  de 
formed)  zone,  one  may  estimate  the  delamina¬ 
tion  fracture  toughness  to  be  4000  ,J/m2  as 
shown  below: 

Glc  =  ^  [hA(l  -  VAj)  +  hB  (1  -  VBj) 

^  +^1-^)1 
This  calculation  implicitly  assumes  that  the  lo¬ 
cal  strain  to  fracture  in  the  FI 85  resin  and  the 
composite  are  the  same,  and  that  the  stress 
distribution  is  also  similar,  at  least  on  average. 

Since  the  measured  value  of  delamination 
fracture  toughness  in  this  system  is  1900  J/m2, 
the  estimated  value  of  4000  J/m2  is  seen  to  be 
quite  excessive.  It  also  strongly  suggests  what 
the  tn-situ  fractography  has  already  indicated, 
namely,  that  preferential  microcrack  nuclea- 
tion  at  the  fiber/resin  interface  leads  to  pre¬ 
mature  failure,  preventing  the  realization  of 
4000  J/m2  that  might  otherwise  be  possible. 
This  suggests  that  more  attention  to  the  inter- 
phase  region  is  required  if  a  greater  fraction  of 
the  intrinsic  toughness  of  the  neat  resin  is  to 
be  manifested  in  delamination  fracture  tough¬ 
ness  of  the  fiber-reinforced  composite  material. 

SUMMARY 

1.  The  delamination  fracture  toughness  of 
the  T6T145/F185  composite  is  1900  J/m2. 
whereas  the  fracture  toughness  of  the  neat 
FI 85  resin  is  8100  J/m2. 

2.  Microcracking  in  the  resin  appears  to  be  a 
significant  deformation  mode  for  load  redistri¬ 
bution  at  the  crack  tip. 

3.  Crack  propagation  occurs  in  a  discontin¬ 
uous  manner  through  microcrack  coalescence 
ahead  of  the  crack  tip. 

4.  In  the  composite,  the  microcracking  is 
coarser  and  coalescence  seems  to  occur  prefer¬ 
entially  at  the  resin-fiber  interface. 

5.  The  height  of  the  deformation  zone,  above 
and  below  the  plane  of  crack  propagation,  is 
greater  in  the  resin  than  in  the  composite. 

6.  Comparison  of  the  relative  heights  of  the 
damage  zone  yielded  4000  J/m2  as  a  calculated 
delamination  G,c  for  the  composite,  which  is  an 
over-estimation  by  a  factor  of  two. 

7.  It  may  be  inferred  that  the  strain  to  failure 
in  the  composite  is  lowered,  probably  due  to 
premature  microcrack  coalescence  at  the  resin 
fiber  interface. 
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INTRODUCTION 

It  has  been  recognized  for  some  time  that  increasing  the  resin  toughness  in 
composite  materials  does  not  give  a  proportionate  increase  in  the  deLami  nation 
fracture  toughness  of  these  materials. *  » ^  it  has  been  suggested  that  the  reason 
for  this  smaller  increment  of  increase  in  delaraination  toughness  compared  to  the 
increment  of  increase  in  toughness  in  the  neat  resin  is  due  to  reduction  in 
crack  tip  strain  possible  in  the  composite  due  to  interference  of  the  fibers. 

The  addition  ot  fibers  to  a  ductile  resin  system  could  potentially  reduce 
the  crack  tip  strain  for  at  least  three  reasons. ^  First,  failure  in  the 
composite  might  occur  at  the  fiber/raatrix  interface  prior  to  the  extraction  of 
the  full  potential  ductility  out  of  the  resin.  Second,  the  presence  of  fibers 
could  give  local  constraint  in  the  vicinity  of  the  crack  tip,  lowering  the  total 
strain  to  failure  observed.  Third,  the  fibers  can  act  as  rigid  filler,  which 
could  reduce  the  plastic  zone  size  in  resins  which  have  a  neat  resin  plastic 
zone  size  which  is  greater  than  the  size  of  the  resin  rich  region  between  plies 
in  the  composite. 

C  h  a  k  a  c h  e  r  y  and  Bradley^  and  Jordan  and  Bradley^  have  observed  qualitatively 
the  above  mentioned  redvietion  in  plastic  zone  size,  comparing  crack  growth  in  a 
ductile  resin  to  «ie lamination  in  a  graphite/epoxy  composite  made  from  the  same 
resin.  Chakachery  and  Bradley'’  attempted  to  quantify  with  a  simple  phenomeno¬ 

logical  model  the  reduction  in  fracture  toughness  in  going  from  a  tough  resin  to 
d  e lamination  In  a  composite  made  from  this  resin.  They  were  unsuccessful  ly  in 
tully  accounting  for  the  observed  reduction  in  fracture  toughness,  most  likely 
because  their  model  implicitly  assumed  that  the  local  strain  field  In  each  case 
was  similar,  with  the  difference  in  toughness  resulting  primarily  from  the 
presence  of  rigid  filler;  namely,  carbon  fibers.  Obviously,  a  more  realistic 
mode l  wo old  need  ro  take  into  account  the  difference  in  the  strain  field  around 


the 


a.ido 


:  k  tip  tor  neat  resin  compared  to  de lamination  of 


composite 


this  resin. 

The  purpose  of  this  study  has  been  to  develop  a  technique  which  would 
permit  the  direct  measurement  of  the  strain  field  around  a  crack  tip  with 
sufficient  resolution  to  see  how  distinctive  the  fields  are  for  cracks  in  neat 
resin  versus  del  ami  nation  cracks  in  composites  made  from  the  same  resin.  The 
approach  used  has  been  to  burn  a  very  fine  dot  map  onto  the  polished  surface 
(perpendicular  to  the  plane  of  the  crack)  of  a  specimen  using  the  electron  beam 
in  a  scanning  electron  microscope.  The  exact  coordinates  for  each  of  these 
points  are  determined  before  and  after  loading  the  specimen,  allowing  the  deter¬ 
mination  of  the  displacement  of  each  point.  This  displacement  field  may  then  be 
differentiated  to  give  the  strain  field  around  the  crack  tip. 


EXPERIMENTAL  PROCEDURE 

With  considerable  trial  and  error,  a  technique  has  been  developed  to  burn  a 
very  systematic  dot  pattern  onto  the  polished  surface  of  the  specimens.  The 
beam  of  the  electron  microscope  had  to  be  adjusted  to  give  a  dots  that  were 
sufficiently  large  to  remain  easily  recognizable  after  considerable  deformation 
but  small  enough  to  give  adequate  resolution.  The  very  regular  dot  pattern  seen 
in  Fig.  1 A  was  burned  onto  the  surface  using  a  JEQL-J5  scanning  electron  micros¬ 
cope  with  the  assistance  of  an  image  processing  system  manufacture  red  by  K  on¬ 
er  o  n .  The  spacing  seen  in  Figure  1  A  is  approximately  10  microns.  A  finer  dot 
spacing  is  possible,  but  not  very  useful  unless  the  size  of  the  dots  can  also  be 

reduced,  while  retaining  their  visibility  after  large  deformation. 

The  resin  system  selected  for  use  In  this  study  was  H  excel  Fid1).  Tills 
resin  was  selected  because  it  is  known  to  be  relatively  ductile  (10*  elongation 
in  ri  tensile  test),  and  therefore  would  be  expected  to  give  significant  strain 
in  the  vicinity  of  the  crack  tip  prior  to  failure.  Furthermore,  this  system  has 
a  large  reduction  in  toughness  in  0  o i n  g  from  the  neat  resin  to  del  .mind  lion  of 
the  composite  (d00()J/m~  to  2  000J/m^)t>.  Thus,  it  was  considered  ideal  lor  this 
study. 

H  excel  F  1 6  3  resin  was  cast  Into  rectangular  plates.  It  was  then  machined 
into  compact  tension  specimens  which  were  subsequently  fatigue  precracked.  Then 
the  specimens  were  polished  on  a  mlcroprocesser  to  a  0.03  micron  finish. 
Finally,  the  dot  map  was  burned  on  the  surface  using  the  electron  beam  in  the 
SEM,  as  shown  in  Figure  1  A.  The  crack  tip  is  seen  at  the  bottom  of  Figure  IA, 


T 


V 


just  below  the  dot  map.  After  loading  which  resulted  in  a  small  amount  of  crack 
growth  which  moved  the  crack  tip  just  into  the  dot  map,  the  deformation  field  is 
clearly  seen  through  the  distortion  of  the  dot  map  (see  Figure  IB).  The  dis¬ 
placement  tie  Id  and  the  associated  strain  field  were  determined  using  the  image 
processing  s  v  s  t  o  in  in  conjunction  with  a  computer  program  written  to  numerically 
ditferentinl  the  displacement  field  data. 

RESULTS  AND  DISCUSSION 

The  reduced  results  are  presented  in  Figures  2-4.  The  normal  strains  as 
well  as  the  principal  shear  strain  are  seen  as  a  function  of  their  respective*  x  - 
y  coordinate  locations  around  the  crack  tip,  which  has  a  coordinate  or  x » b  9 . 1  , 
v=  57.48.  The  compact  tension  specimen  was  loaded  In  the  x-direct  ion  with  crick 
growth  in  the  v-direetion. 

It  should  be  emphasised  that  the  three  dimensional  p  l  o  C  t  i  n  g  routine  used  t  .> 
present  the  results  graphically  does  local  smoothing  prior  to  plotting.  The 
polynominal  curve  fit  for  the  smoothing  operation  interpolates  as  well,  giving  a 
much  finer  grid  in  the  graphical  results  than  was  actually  utilized  (cuiur.ht 
actual  grid  size,  Figure  1A,  to  grid  size  seen  in  the  three  dimensional  plotting 
(Figures  2-4). 

The  elastic  strain  at  yielding  in  F  L  8  5  as  measured  in  a  tensile  test  i* 
about  1.2*.  Thus,  the  region  of  nonlinear  deformation  around  the  crack  tip  is 
seen  to  bo  quite  extensive.  If  it  is  remembered  that  carbon  fibers  typical  have 
diameters  or  about  7  microns,  then  the  nonlinear  deformation  zone  would  extend 
to  at  least  five  fiber  diameters  above  and  below  the  plane  of  del  a  mi  nation  if 
the  composite  has  a  similar  strain  field  around  the  tip  of  a  growing  delamina¬ 
tion  crack  to  that  observed  in  Figs.  2-4  for  the  FIBS  neat  resin.  Qualitative 
indications  of  a  deformation  zone  of  this  size  in  a  delaminating  T  b  T  1  4  ■>  /  F  1  8  5 
composite,  but  without  direct  measurement  of  the  strain  field,  have  been  noted 
by  Bradley  and  Cohen  . 

A  second  important  observation  is  the  magnitude  of  the  local  strain  at  the 
crack  tip.  The  elongation  measured  in  a  tensile  test  is  10%,  with  no  apparent 
necking.  Thus,  one  might  have  assumed  for  micromechanics  modelling  that  the 
local  strain  to  failure  would  be  10%,  or  possibly  less  due  to  the  effects  of 
constraint.  However,  it  Is  clear  that  the  tensile  elongation  greatly  under¬ 
estimates  the  true  local  strain  to  failure  at  the  tip  of  a  growing  crack  in 
FIS').  The  measured  value  of  48%  should  be  considered  a  lower  bound  since  It  is 


measured  over  a  It)  micron  gage  at  the  crack  tip  in  a  very  steep  strain  gradient 
and  would  increase  to  an  even  higher  value  as  crack  advance  gave  ductile  frac¬ 
ture  to  this  material.  The  critical  value  for  local  strain  to  failure  for 
delamination  crack  growth  is  probably  somewhat  lower  than  the  value  for  crack 
growth  in  the  neat  resin  for  reasons  previously  discussed. 

One  concern  that  should  be  mentioned  is  the  question  of  whether  the  surface 
strain  measurements  are  in  any  way  representative  of  the  bulk  strain  distribu¬ 
tion.  In  materials  where  cavitation  or  microcracking  does  not  occur,  one  would 
expect  a  significant  difference  between  the  strain  distribution  for  the  plane- 
stress  state  of  stress  that  exists  at  the  surface  and  the  strain  distribution 
that  wouLd  exist  for  the  plane-strain  state  of  stress  in  the  bulk.  However, 
Kidd  with  rubber  particle  toughening  does  cavitate,  relaxing  the  through  thick¬ 
ness  stress.  Thus,  for  this  system  we  believe  the  surface  strain  distribution 
is  representative  of  the  bulk  s  t  r  a  i  n  discribution. 

Comparable  results  for  crack  growth  in  a  delaminating  composite  were  not 
available  at  cl;  e  time  of  the  preparation  of  this  manuscript  but  should  be 
available  tor  the  presentation  at  the  meeting 
SUMMARY 

A  technique  for  the  direct  measurement  of  the  strain  field  around  a  growing 
crack  has  been  developed.  Results  on  K  1  8  5  resin  indicate  a  large  non-linear 
Jetor  na:  ion  cone  .irounJ  the  crack  tip  and  a  surprisingly  large  local,  crack,  tip 
strain  to  failure. 
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ABSTRACT 

The  Stress  field  ahead  of  the  crack  tip  of  a  split  laminate  specimen 
.oa  dec  oncer'  moce  I  and  Mode  II  de  lam:  rat  ion  has  been  determined  by  means 
of  a  finite  element  analysis.  The  results  have  then  been  compared  to  the 
observed  damage  zone  developed  ahead  of  the  crack  tip  for  the  same  loading 
conditions.  A  direct  correspondence  has  been  found  between  the  stress 
fields,  damage  zones  observed,  and  resistance  to  delamination  of 
composites  for  this  modes  of  delamination  fracture. 

INTRODUCTION 

Fiber-reinforced  composite  materials  usually  offer  poor  resistance  to 
delamination.  As  a  result,  a  complete  understanding  of  the  delamination 
process  is  needed  to  properly  design  composite  structures  and  develop 
materials  with  improved  fracture  toughness  characterist ics .  An  important 
area  of  investigation  is  the  determination  of  the  stress  field  ahead  of 
the  crack  tip  for  the  various  modes  of  delamination  fracture  (i.e.  mode  I, 
mode  II,  and  mode  III),  by  means  of  finite  element  analyses.  Furthermore, 
these  stress  fields  can  be  related  to  the  corresponding  observed  damage 
zones  that  develoo  ahe±d  of  the  growing  crack  as  determined  from  in-s:tu 
observations  in  the  scanning  electron  microscope.  This  combination  of 
analyses  with  direct  c iservation  of  delamination  has  been  very  useful  in 


cevelopmg  ar.  overall  uncer  stare. ''.g  of  " r  -  r*e.at  .or  .£•* 

rr.acrosoop:  o  res  istarce  :o  crack  propagation  arc  roc-c  of  ce  1  a.T. ;  rat  .  ~r. 
fraco-re . 

lr.  tr.s  investigation ,  tne  stress  field  --.esc  :f  tr.e  crac*  tip  0 
split  laminate  specimen  loadec  under  rncce  1  arc  nose  1 1  del  am  mat .  ;r. 
conditions  das  been  determined  by  means  of  a  finite  e.e.T.ert  analysis.  Tne 
results  were  then  compared  to  tne  observed  carnage  zone  deve.opeo  anead  of 
tr.e  cracK  tit  for  the  same  loading  conditions. 

MATERIALS  AJID  METHODS 

The  determination  of  the  stress  field  aneac  of  tne  crack  tip  of  the 
split  laminate  specimen  loaded  under  mode  !  anc  moce  II  was  accomplished 
oy  first  generating  a  two  dimensional  mesh  '52m  rmr  long  anc  2.54  mm  thick 
consisting  of  79"  nodes  and  29*  elements  as  shown  in  Figure  'a.  Triangular 
six-noded  elements  were  used  around  the  crack  tip  with  mid-side  nodes 
displaced  to  the  quarter  point  [1],  and  a  substantial  refinement  of  the 
mesh  around  the  crack  tip  was  made  to  overcome  the  difficulty  imposed  by 
the  stress  singularity  present  at  the  crack  tip.  This  has  been  a  linear 
analysis  made  with  a  finite  element  algorithm  developed  by  Henriksen  [2]. 
The  algorithm  is  based  on  a  nonlinear  code  updated  with  Lagrangian 
formulation  using  six  and/or  eight  noded  isoparametric  elements  with  two 
degrees  of  freedom  per  node.  Elastic  constants  for  a  unidirectional 
laminate  of  AS4/3502  graphite/epoxy  composite  to  be  used  in  the  analysis 
were  obtained  from  Hercules  [3]. 

Mode  I  loading  was  simulated  by  applying  a  symmetric  load  at  the 

cracked  end  of  the  mesh  as  shown  in  Figure  lb.  Mode  II  loading  was 

introduced  by  asymmetrically  loading  the  cracked  end  (see  Figure  1c).  The 
load  level  applied  corresponds  to  approximately  the  load  at  onset  of  crack 
growth  for  these  two  modes  of  failure,  as  determined  from  experimental 
measurements  in  combination  with  beam  linear  beam  theory  (4).  Stress 

contour  plots  were  obtained  from  the  output  data  and  the  stress 
distribution  was  plotted  as  a  function  of  distance  ahead  of  the  crack  tip. 

Two  sizes  of  split  laminate  specimens  were  used  for  direct 

observation  of  deiamination  in  the  scanning  electron  microscope,  one  38.1 
mr  long,  5.08  mm  wide  and  1.14  mm  thick,  and  the  ocher  30.48  mm  long,  5.08 
mm  w:d-,  and  1.14  mm  thi"k.  Hercules  unidirectional  AS4/3501-6 
graphite -epoxy  composite  material  was  used  for-  mode  I  and  mode  II 
observations  of  the  damage  zones  respectively.  The  real  time  observations 


of  ce  lamina  Cion  were  made  ir.  a  UEOL-35  scanning  electron  microscope.  Mode 
1  de  lam  mat  ion  was  achieved  b_,  pushing  a  wedge  into  the  precracked  portion 
of  the  longer  split  laminate  specimens,  using  the  tensile  stage  of  the 
scanning  electron  microscope.  The  wedge  was  sufficiently  blunt  to  ensure 
that  it  remained  well  away  from  the  crack  tip,  giving  esentialiy  pure  mode 
I  loading.  Mode  II  delamination  was  achieved  by  means  of  a  specially 
designed  three  point  bend  fixture  installed  to  the  tensile  stage  of  the 
scanning  electron  microscope.  All  surfaces  observed  were  coated  with  a  150 
A  thick  gold/'paiiadium  film  to  minimize  charging  effects  associated  with 
the  nonconduct  ive  nature  of  the  eooxy  matrix  of  the  composite.  The 
experimental  results  were  recorded  on  standard  tn-x  film. 

RESULTS  AVD  DISCUSSION 

Figures  2  anc  3  show  the  finite  element  results  for  both  mode  I  anc 
mode  II  loading.  Figure  2a  is  a  Syy  (normal)  stress  contour  plot  from  the 
vicinity  of  the  crack  tip  for  moce  I  loading.  The  normal  stress  Syy  which 
for  this  loading  and  geometry  is  the  principal  normal  stress  S£ ,  rapidly 
decreases  ahead  of  the  crack  tip.  Figure  2b  corresponds  to  the  Sxy  (shear) 
stress  contour  plot  for  a  mode  II  condition.  Note  how  the  shear  stress 
drops  off  much  more  slowly  with  distance  ahead  of  the  crack  tip  than  does 
the  normal  stress  for  mode  I  loading  (Figure  2a).  Furthermore,  the  shape 
of  the  stress  field  is  more  narrow  and  elongated.  For  this  loading 
condition,  Sxy  =  maximum  shear  stress  =  S • .  Figure  3  shows  the  stress 
distribution  as  a  function  of  distance  ahead  of  the  crack  tip  for  mode  I 
and  mode  II.  In  the  case  of  the  normal  stress  for  mode  I,  the  stress  drops 
off  rapidly  until  it  is  compressive  at  a  distance  of  0.76  tnn  ahead  of  the 
crack  tip.  Finally,  it  gradually  approaches  a  zero  stress  level  at 
approximately  8  mm  from  the  crack  tip  which  is  maintained  all  the  way  to 
the  end  of  the  beam.  It  should  be  noted  that  the  compressive  stresses 
observed  are  not  expected  to  significantly  influence  the  fracture  mode  at 
the  crack  tip  because  it  develops  far  enough  away  from  the  crack  tip.  In 
the  case  of  the  shear  stress  distribution  ahead  of  the  crack  tip  for  mode 
II,  it  monotonically  decreases  to  a  constant  value.  As  it  can  be  seen,  the 
shear  stress  ahead  of  the  crack  tip  for  mode  II  loading  decays  much  slower 
than  the  normal  stress  for  mode  I  loading.  These  results  indicate  that  the 
stress  concetration  at  the  crack  tip  is  distr.buted  over  a  larger  distance 
for  mode  II  loading  than  for  mode  I  loaning.  For  mode  II,  the  extent 
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Figure  2.  Stress  contour  plots. 


of  the  stressed  material  above  and  below  the  center  plane  being  much 
smaller  than  for  mode  I,  as  seen  in  the  stress  contour  plots  in  Figure  2. 

Figure  4  corresponds  to  the  observed  damage  zones  developed  ahead  of 
the  crack  tip  for  both  mode  I  and  mode  II  delamination  of  AS4/3501-6.  The 
mode  I  damage  zone  is  characterized  by  extensive  microcracking  of  the 
resin  extending  at  least  90  microns  ahead  of  the  macroscopic  crack  tip  and 
approximately  7.5  microns  above  and  below  the  primary  plane  of  crack 
advance.  In  the  case  of  mode  II,  the  damage  zone  extends  much  further 
ahead  of  the  crack  tip  (at  least  180  microns)  than  the  mode  I  damage  zone, 
and  is  approximately  5  microns  in  extent  above  and  below  the  primary  crack 
plane. 


w 


T 


Figure  3.  Stress  distribution  ahead  of  crack  tip  for  moce  I  and  mode  II. 

A  careful*  comparison  of  the  stress  fields  of  Figures  2,3  and  the 
damage  zones  shown  in  Figures  4  for  mode  I  and  mode  II  delamination, 
reveals  a  direct  correspondence  between  them.  The  more  gradual  decrease  In 
the  magnitude  of  the  shear  stress  with  distance  from  the  crack  tip  for 
mode  II  (Figure  3),  along  with  the  elongated  shape  of  the  shear  stress 
contour  plot  seen  in  Figure  2b  indicates  that  the  stress  concentration  at 
the  crack  tip  is  distributed  over  a  larger  range,  causing  the  size  of  the 
damage  zone  ahead  of  the  crack  tip  to  be  elongated  (Figure  4b).  In 
contrast,  the  rapid  decay  of  the  magnitude  of  the  normal  stress  for  mode  I 
with  distance  from  the  crack  tip  (Figure  3)  and  the  more  circular  shape  of 
the  normal  stress  contour  plot  of  Figure  2a  indicate  that  the  stress 
concentration  at  the  crack  tip  is  much  more  localized,  causing  the  extent 
of  the  damage  zone  ahead  of  the  crack  tip  to  be  smaller  than  for  mode  II 
delamination.  Furthermore,  the  extent  of  the  stressed  material  above  and 
below  the  plane  of  delamination  is  seen  in  the  stress  contour  plots  to  be 
much  smaller  for  mode  II  than  for  mode  I.  Thus,  the  damage  zone  (Figure  4  ) 
changes  from  wide  to  narrow  as  one  goes  from  mode  I  to  mode  II.  These 
findings  indicate  that  the  greater  resistance  to  delamination  for  mode  II 
loading  (CTjc  =  570  J/m2)  than  for  mode  I  loading  (Gjc  =  190  J/m^)  in 
AS4/3502,  a  similar  material  as  AS4/3501-6,  may  be  explained  in  part  as 
resulting  from  the  different  decay  rate  of  the  stress  field  and  the 
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Figure  5.  Von-Mises  stress  for  mode  I  and  mode  II  loading. 


obviously  from  the  edge  of  the  specimen  where  a  plane  stress  condition 
occurs,  "herefore,  in  order  to  have  a  first  order  approximation  as  to  the 
effect  tie  difference  in  the  state  of  stress  at  the  surface  of  the 


specimen  where  our  direct  observations  are  made,  and  at  the  center  were 
plane  strain  conditions  prevail,  the  Von-M:ses  stress  for  mode  I  and  mode 
II  was  plotted  as  function  of  distance  ahead  of  the  crack  tip.  The  results 
in  Figure  4a  indicate  that  for  mode  I,  the  difference  in  the  Von-mises 

stress  (or  root  mean  square  shear  stress)  in  plane  stress  versus  plane 
strain  is  approximately  10%  .  In  the  case  of  mode  II,  the  difference  is 
negligible  as  evidenced  in  Figure  4b.  Therefore,  it  appears  that  for  these 
orthotropic  materials,  the  difference  between  plane  stress  and  plane 

strain  is  not  too  significant  and  our  surface  observations  of  damage  zone 
size  and  detail  should  not  be  very  different  than  subsurface  behavior. 
Previous  studies  in  which  the  post  mortem  appearance  of  the  fracture 

surface  near  a  free  edge  to  that  at  the  center  of  a  specimen  are 

consistent  with  this  interpretation  [6],  A  final  observation  from  Figure  4 
is  that  the  stress  field  ahead  of  a  crack  tip  decays  much  more  slowly  for 
and  orthotropic  material  with  the  crack  parallel  to  the  fibers  than  for  an 
isotropic  material . 

CONCLUSION 

Ir.  summary,  it  has  beer,  found  that  the  different  size  and  shape  of 
the  damage  zone  observed  when  comparing  mode  I  to  mode  II  delamination  is 
consistent  with  the  difference  in  the  size  and  shape  of  the  stress  field 
predicted  by  the  finite  element  analysis.  Furthermore,  these  differences 
are  also  consistent  with  the  difference  in  resistance  to  delamination  of 
brittle  graphite/epoxy  composites  for  the  two  types  of  delaraination 
studied.  Also,  there  appears  to  be  a  very  small  difference  between  plane 
stress  and  plane  strain  states  of  stress  for  orthotropic  materials. 
Therefore,  the  damage  zones  observed  during  in-situ  delaraination  in  the 
scanning  electron  microscope  should  not  vary  significantly  from  the  actual 
damage  zone  that  develops  at  the  center  of  the  specimens  (plane  strain 
damage  zone) . 
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ABSTRACT :  Multi-specimen  J  integral  tests  which  are  defined  in 
the  Standard  Test  for  Jlc,  A  Measure  of  Fracture  Toughness  ( ASTM 
E  813-81)  have  been  conducted  on  side-grooved  compact  tension 
specimens  of  medium  and  high  density  polyethylene.  Specimens  were 
precracked  with  a  razor  blade,  as  has  been  common  practice  to 
date,  or  were  fatigue  pre-cracked.  Tests  were  performed  over  a 
50X  range  of  crosshead  rates  in  the  quasi-static  regime.  Results 
indicated  that  the  razor  notching  gives  a  J-R  curve  which  is 
significantly  higher  chan  that  obtained  using  fatigue  precracked 
specimens.  A  significant  variation  in  the  J-R  curve  was  noted  for 
changes  in  the  displacement  rate  due  to  the  viscoelastic 
character  of  the  polyethylene.  The  orientation  of  the  specimens, 
which  were  cut  from  actual  pipe,  was  found  to  have  a  significant 
effect  on  the  fracture  toughness.  The  fracture  toughness  results 
were  correlated  with  f ractog raphic  results  and  the  behavior  of 
side-grooved  vs.  non-side-grooved  specimens  was  compared. 

KEY  WORDS :  Polyethylene,  Fracture  (materials).  Compact  Tension 


Test,  Fatigue  (materials) 


INTRODUCTION 


Thermoplastics  are  seeing  increased  use  as  structural 
materials.  With  this  increased  structural  use  a  need  has  arisen 
for  a  method  of  determining  useful  fracture  mechanics  parameters 
for  failure  prediction.  Thermoplastics  such  as  polyethylene  (PE) 
which  have  a  high  degree  of  ductility  and  exhibit  strong 
viscoelastic  behavior  present  a  problem  for  linear  elastic 
fracture  mechanics  (LEFM)  testing  procedures.  Previous 
investigators  [1,2]  have  used  the  concepts  of  LEFM  with  tough  PE 
in  an  effort  to  determine  a  valid  plane  strain  fracture  toughness 
(Klc).  They  encountered  difficulty  in  doing  so  at  temperatures 
near  normal  operating  conditions  for  PE  piping.  However,  it  has 
been  shown  [3]  that  an  elast ic-plast ic  J  integral  approach  as 
defined  in  the  Standard  Test  for  J-,  ,  A  Measure  of  Fracture 
Toughness  ( ASTM  E  813-81)  will  yield  plane  strain  toughness 
values  for  tough  PE  which  meet  the  requirements  of  the  standard 
for  valid  J  values  when  tested  in  the  normal  operating 
temperature  regime. 

The  question  then  arises  as  to  the  validity  of  a  non  linear 
elastic  analysis  such  as  the  J  integral  when  applied  to  a 
viscoelastic  material.  There  is  no  currently  approved  standard 
for  determining  the  fracture  toughness  of  viscoelastic  materials. 
However,  R .  A .  Schapery  ]  4  ]  has  shown  that  tie  J  integral  may  be 
applied  to  fracture  in  viscoelastic  media.  Schapery  assumes  in 
his  analysis  that  the  work  of  fracture  in  the  process  zone  is 
independent  of  the  crack  growth  rate  and  that  the  significant 
viscoelastic  effect  (rate  dependence)  is  through  the  time 


dependent  behavior  in  the  stress  field  which  surrounds  the 

process  zone. 

A  number  of  questions  remain  unresolved  in  the  application 
of  fracture  mechanics  testing  of  tough  thermoplastics  such  as  PE. 
This  study  examined  the  effects  of  several  experimental 
parameters  on  the  apparent  J.  tougnness  of  PE  tested  at  room 
temperature.  Test  rate,  precracx  acuity,  and  specimen  thickness 
were  examined.  Previous  investigators  [1-3;  have  employed 
primarily  a  three  point  tend  configuration  in  their  tests.  In 
this  study  a  compact  tension  ■  CT '  con.f  igurat  ion  as  defined  in 
A  STM  E  813-81  was  .sec  for  al.  tests.  Due  to  the  highly  drawn 
nature  of  pipe  materials,  a  directional  dependence  in  tougnness 
of  actual  pipe  material  was  suspected.  Experiments  were  performed 
to  confirm  this.  In  addition  to  alignment  of  molecules,  extrusion 
of  pipe  usual.y  introduces  residual  stresses  which  may  influence 
the  effective  toughness  of  pipe  in  service.  Previous 
investigators  have  worked  primarily  from  sheet  stock  or 
flattened,  annealed  pipe  in  preparing  specimens.  Since  actual 
field  performance  is  of  interest  with  pipe  materials,  oot.n 
annealed  and  unannealed  specimens  were  tested  m  this  study. 

PROCEDURE 

Equipment-Eractu  re  testing  was  performed  on  a  servo- 
hydraulic  closed  loop  test  stand  with  load-lire  di  spl  aceme.o  t 
measured  by  a  linearly  varying  differential  transformer  mounted 
on  the  specimen  holding  pms.  Data  aquisit.cn  and  reduction  were 
accomplished  with  a  micro-computer 


connected 


to  an  analog  t  o 


digital  signal  processer.  Specimens  were  observed  during  resting 
through  an  optical  stereo  microscope  and  post  mortem  fractography 
was  performed  in  a  scanning  elecrron  microscope  (SEM). 

Specimen  Preparation-2.54  cm  CT  specimens  with  uncracked 
ligament  lengths,  b,  of  1.5  cm  were  prepared  from  a  high  density 
pipe  grade  polyethylene  PE  3306  described  in  the  Standard 
Specification  for  Thermoplastic  Gas  Pressure  Pipe,  Tubing,  and 
Fittings  (AST.M  E  2513  —  75b)  .  Son-annealed  specimens  were  produced 
in  two  thicknesses,  0.635  cm  and  0.556  cm,  which  were  reduced  20 
percent  by  side-grooving  to  0.503  cm  and  0.445  cm  respectively. 
Annealed  specimens  were  cut  from  stock  which  had  been  heated  to 
130  °C,  flattened,  and  then  annealed  for  two  hours  at  that 
temperature.  Specimens  cut  from  this  stock  were  0.572  cm  thick 
after  being  reduced  20  percent  by  sice-grooving.  Two  add.tional 
specimens  were  cut  to  a  thickness  of  0.508  cm.  One  was  reduced  in 
tniCKness  20  percent  with  sidecrocves.  Wit n  the  exception  of  one 
set,  specimens  were  cut  wit  n  t.oe  starting  notch  in  the 
longitudinal  pipe  direction;  thus,  they  are  termed  transverse 
specimens  in  the  standard  ASTV.  nomenclature.  The  remaining  set  of 
specimens  was  cut  to  test  transverse  pipe  properties  and  are 
termed  longitudinal  specimens  in  the  AST.M  nomenclature .  The  angle 
of  the  sidegrooves  was  45  degrees  and  where  fatigue  precracking 
was  used  specimens  were  precracked  before  side  grooves  were 
machined.  Razor  notched  specimens  were  side  grooved  before 
notching . 


Fat ique- Specimens  were  fatigue  precracked  in  accordance  with 
the  procedure  outlined  in  ASTM  E  813-81  with  the  exceptions  noted 
below.  It  was  found  that  fatigue  at  40  percent  of  the  limit  load 
specified  by  the  standard  produced  no  noticeable  sharpening  of 
the  crack  tip  and  no  crack  advance.  Therefore  the  requirement 
that  the  final  stage  of  crack  growth  be  at  this  level  was 
ignored.  Specimens  were  fatigued  at  a  rate  of  ten  Hz  and  a  load 
of  70  percent  of  ?r  for  30  000  cycles  to  produce  the  precrack. 
Sidegrooves  were  cut  after  precracking  because  it  was  found  that 
the  fatigue  crack  front  tended  to  run  ahead  at  the  surface  of  the 
specimen  when  they  were  machined  prior  to  fatiguing. 

Razor  N’o t  ch  i  r.c-Ra  zo  r  notches  were  made  using  the  test  stand 
and  specially  designed  fixtures.  A  new  razor  blade  was  used  for 
each  specimen  to  obtain  one  sharpest  possible  crack  tip. 

Tests-Tests  were  run.  at  an  ambient  laboratory  temperature  of 
2*  °  C .  The  standard  crosshead  rate  for  tests  was  0.0254  cm/min. 
One  set  of  tests  was  run  at  0.00254  and  another  at  0.127  cm/min 
to  examine  rate  effects.  Immediately  after  loading,  specimens 
were  immersed  in  liquid  nitrogen  and  then  broken  in  the  test 
stand.  The  liquid  nitrogen  fracture  produced  clear  boundaries 
between  the  actual  crack  growth,  the  Dugdale  type  process  zone, 
and  the  undamaged  material.  Approximate  measurements  of  the  crack 
advance  were  made  during  the  tests  using  an  optical  microscope  to 
permit  a  good  distribution  of  crack  lengths  which  would  result  in 
a  more  reliable  J-R  curve.  Final  crack  length  measurements  were 
made  from  micrographs  taxer,  m  t.ne  SEM.  To  permit  viewing  in  the 
SEM,  one  half  of  each  specimen  was  sputter  coated  with  150A  of  a 
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gold-palladium  alloy.  The  crack  advance  visible  on  the  resulting 
micrographs  was  measured  at  eleven  points  evenly  spaced  across 
the  specimen  width.  The  mean  of  the  nine  inner  points  and  the 
average  advance  at  the  two  surfaces  was  taken  as  the  crack 
advance.  The  micro-computer  used  to  record  load  displacement  data 
was  also  used  to  calculate  the  areas  under  the  load  displacement 
curves.  J^c  was  found  from  the  resistance  curve  plotted  according 
to  the  procedure  in  ASTM  E  813-81. 

RESULTS 

J  Integral  Testing-Resistance  curves  li<e  the  one  shown  in 
Figure  1  were  generated  for  each  test.  The  results  are  summarized 
in  Table  1.  The  non-s i deg rooved  thin  specimen  (0.508  cm)  showed  a 
high  degree  of  crack  tip  tunneling  which  made  it  unusable  in  the 
J  determination.  However,  the  identical  specimen  sidegrooved  20 
percent  to  a  thickness  of  0.4064  cm  yielded  a  straight  crack 
front.  The  sidegrooved  thin  specimen  yielded  a  J  vulue  that  fell 
just  below  the  J-R  curve  for  the  thicker  specimens  used  in  the 

i 

remainder  of  the  tests. 

ASTM  E813-81  establishes  minimum  dimensions  for  specimens 
which  are  necessary  in  testing  of  metals,  for  the  J  analysis  to 
be  valid.  The  requirements  are: 

a.  b ,  3  >  25  J-^/S 

b.  dJ  <  S 
da  y 

These  have  been  calculated  for  each  of  the  tests  and  are  listed 
in  Table  1.  Though  the  yield  strength  of  PE  can  vary  greatly  with 


strain  rate,  all  tests  easily  satisfied  the  second  requirement. 
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However,  according  to  the  standard,  the  unannealed  razor  notched 
CT  specimens  were  slightly  smaller  than  the  specified  minimum. 
This  was  also  the  case  in  the  high  crosshead  rate  test  using 
annealed  material.  These  values  are  reported  for  comparison  since 
the  standard  does  not  strictly  apply  to  this  material.  Tests 
satisfying  these  requirements  were  not  possible  since  thicker 
stock  was  noc  available. 

Fatigue  precracking  resulted  in  a  drop  in  J^c  of  33  percent 
when  compared  to  razor  notched  specimens  in  both  annealed  and 
unannealed  samples.  Annealing  decreased  the  measured  J^c  by  47 
percent  in  both  fatigue  precracked  and  razor  notched  tests. The 
longitudinal  specimens  yielded  a  J,c  which  was  11  percent  greater 
than  that  determined  for  the  transverse  specimen  tests. 

Blunting  Line-In  several  of  the  multi-specimen  tests,  one 
test  was  terminated  before  the  onset  of  crack  advance  across  the 
entire  crack  front.  The  resulting  points  should  have  fallen  on 
the  blunting  line.  As  can  be  seen  in  Figure  2,  the  high  rate  test 
required  a  steeper  blunting  line.  When  a  yield  stress  measured  at 
a  strain  rate  which  corresponded  to  the  rate  of  the  fracture  test 
was  used,  the  appropriate  blunting  line  resulted. 

F r ac t oq r aph y -Compa r i son  of  the  surfaces  of  the  low  rate  and 
medium  rate  fracture  tests  shows  a  marked  similarity  in 
morphology  but  a  significant  difference  in  the  scale  of  features. 
(Figures  3  and  4).  The  difference  in  size  of  fibrils  and  voids 
corresponds  to  the  trend  identified  by  Chan  and  Williams  [5]  and 
Lee  and  Epstein  [6’.  The  lower  rate  has  smaller  voids  and  shorter 
fibrils.  The  average  diameters  of  the  voids  were  measured  from 


SEM  micrographs  by  Lowry  [7]  and  are  listed  in  Table  2. 

Of  special  interest  was  the  observation  that  the  Dugdale 
type  process  zone  varied  in  length  with  test  rate.  The  average 
lengths  of  process  zones  in  three  annealed  razor  notched  specimen 
tests  are  listed  in  Table  2.  This  variation  was  not  what  one 
would  expect  if  the  yield  stress  rate  sensitivity  was  primarily 
responsible  for  increases  in  toughness  while  the  work  done  in  the 
process  zone  remained  constant  with  changing  rate.  In  fact,  since 
the  yield  strength  increases  with  strain  rate,  the  process  zone 
should  have  gotten  smaller  as  strain  rate  increased.  This  would 
have  been  due  to  the  increased  constraint  induced  by  the  higher 
yield  strength  of  the  surrounding  material.  This  is  clearly 
illustrated  by  the  Dugdale  equation 

r/a  =  sec  ( <l<x/2Sy )  -  1 

Since  Sy  is  increasing  in  the  surrounding  material  with  the 
increasing  rate,  r  should  be  decreasing.  Pest  mortem  examination 
indicated  that  this  is  not  the  case  in  this  material. 

DISCUSSION 

Crack  advance  was  observed  to  follow  the  process  described 
by  S.K.  Bhattacharya  [8].  A  Dugdale  type  region  formed  ahead  of 
the  crack  tip.  Within  this  process  zone  voiding  begins.  These 
voids  grow  and  coalesce  creating  a  fibrillated  structure.  Further 
crack  tip  opening  causes  fibril  failure  and  the  crack  advances. 

Annealinq-The  47  percent  drop  in  J^c  which  occured  with 
annealing  is  extremely  significant.  This  may  be  due  to  an 
accelerated  physical  aging  caused  by  the  anrealing  process.  It  is 
also  possible  that  annealing  relieved  residual  stresses  which 
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inhibit  crack  advance  in  the  unannealed  material.  J.G.  Williams 
[9]  has  measured  significant  residual  stresses  in  extrudea  pipe. 
In  fact,  the  annealed  specimens  yielded  a  symmetric  crack  front 
while  unannealed  material  exhibited  uneven  crack  advance 
suggestive  of  the  presence  of  residual  stresses.  These  residual 
stresses  could  have  been  responsible  for  the  higher  resistance  to 
crack  growth  in  the  original  pipe  material.  It  is  also  possible 
that  annealing  modified  the  crystalline  structure  of  the 
material.  Density  measurements  indicated  a  change  in 
crystallinity  of  less  than  one  percent.  However,  spherulite  si2e 
may  have  been  changed  by  the  heat  treatment,  thus  altering  the 
toughness.  The  annealing  technique  is  commonly  used  when  working 
with  pipe  material  but  these  results  indicate  that  material  so 
treated  may  not  reflect  true  pipe  behavior. 

Rate-The  effect  of  rate  on  the  magnitude  of  Jlc  was  much 
greater  than  expected.  The  300  percent  increase  in  magnitude  of 
Jlc  observed  with  the  increase  in  test  rate  from  the  intermediate 
rate  to  the  high  rate  was  surprising.  From  linear  viscoelasticity 
it  can  be  shown  that  crack  growth  rate  is  related  to  stress 
intensity  by: 

a  =  K2«1  +  Vn) 

if  fracture  energy  is  not  a  function  of  rate.  For  PE  n  can  be 
taken  to  be  approximately  0.2  and  the  equation  becomes: 

k;c  -  a‘/n 

By  using  the  relationship  between  K  and  J  for  linear  elastic 
materials  we  can  say  that  some  equivalent  "J^c"  is  thus 
proportional  to  a1^.  Thus  a  fifty  fold  increase  in  test  rate 


(crack  growth  rate)  will  yield  a  90  percent  increase  in  J-,c. 
Clearly  this  is  not  adequate  to  account  for  the  increase  in 
measured  J  with  the  increasing  test  rate.  Thus  Schaper^'s  [4] 
assumption  that  the  fracture  energy  is  independent  of  1  crack 
growth  rate  does  not  hold  in  this  system.  Stated  another  way,  the 
far  field  viscoelastic  effects  are  not  sufficient  to  account  for 
the  300  percent  increase  in  measured  J1c.  Some  two  thirds  of  the 
increase  remain  unaccounted  for. 

The  changes  in  measured  J  _  which  correspond  to  changes  in 
the  length  and  surface  morphology  of  the  process  zone  suggest 
that  the  fracture  process,  and  therefore  the  energy  dissipated  in 
the  process  zone,  is,  in  fact,  rate  dependent.  The  J  integral  can 
be  related  to  other  material  and  fracture  parameters.  If  we  take 
the  standard  relationship  between  J  and  CTOD: 

J  =  m'Sy’dCT0D 

and  assume  that  m  is  rate  independent  then  J  can  be  increased  by 
increasing  the  yield  stress,  Sy,  or  the  CTOD. 

As  has  already  been  stated,  yield  stress  increases  with 
increasing  strain  rate.  The  strain  rate  dependence  of  the  stress- 
strain  behavior  of  PE  is  demonstrated  in  Figure  5.  The  increase 
in  yield  strength  and  rate  of  strain  hardening  as  the  strain  rate 
is  increased  are  clearly  apparent.  However  the  increase  in  yield 
stress  over  the  range  of  strain  rates  used  i'r.  these  tests  is  less 
than  25  percent.  This  is  not  nearly  enough  to  account  for  the  200 
percent  increase  in  J  which  remains  after  far  field 
viscoelasticity  is  taken  into  account. 

From  the  above  equation  it  is  clear  that  an  increase  in  CTOD 


can  account  for  the  increase  in  J.  The  voids  in  the  process  zone 
are  ellipsoids,  elongated  in  the  direction  of  the  load  applied  to 
the  specimen.  However,  if  it  is  conservatively  assumed  that  the 
voids  in  the  process  zone  are  spherical,  then  the  diameter  of  the 
voids  on  the  fracture  surface  can  be  used  as  a  lower  bound  for 
the  CTOD .  Comparison  of  the  average  diameters  of  the  voids 
produced  by  the  various  test  rates  indicates  a  900  percent 
increase  in  void  size  over  the  range  of  rates  tested.  This  is 
much  more  than  necessary  to  explain  the  increase  in  J  but  is 
evidence  that  the  majority  of  the  increase  in  the  measured  J  with 
increasing  strain  rate  comes  from  an  increasing  CTOD. 

The  increase  in  CTOD  has  two  possible  complementary 
explanations.  The  first  arises  from  the  thermally  activated 
nature  of  the  fracture  process.  A  high  rate  of  crack  growth 
necessitates  the  breaking  of  more  bonds  and  void  formation  from 
fewer  sites  than  a  lower  rate  which  allows  time  for  thermally 
activated  disentanglement  of  polymer  chains.  Thus  high  rates 
produce  larger  and  fewer  voids  and  fibrils  with  a  correspondingly 
larger  CTOD  which  results  when  the  larger  masses  of  material  are 
drawn  to  failure.  The  second  possible  explanation  for  the 
increase  in  CTOD  is  the  possibility  that  any  decrease  in  final 
elongation  of  the  fibrils  due  to  the  n.igher  rate  may  be 
overwhelmed  by  an  increase  in  the  amount  of  material  drawn  into 
the  process  zone.  The  higher  yield  strength  of  the  process  zone 
material  under  faster  load  rates  could  cause  additional  low 
strength  material  to  be  drawn  in  from  the  surrounding  region. 
This  process  is  illustrated  in  Figure  6.  Material  in  the  fibrils 
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has  been  highly  elongated  and  its  strength  would  correspond  to 
point  B  in  Figure  5.  The  material  on  the  edge  of  the  process  zone 
is  less  severely  deformed  and  would  perhaps  correspond  to  the 
material  condition  at  A  in  Figure  5.  At  the  higher  strain  rates 
the  difference  in  strengths  is  much  more  pronounced,  as  can  be 
seen  in  Figure  5,  so  more  material  will  be  drawn  into  the  fibril. 

Pr ecr ack ing-The  33  percent  drop  in  toughness  between  razor 
notched  and  fatigue  precracked  specimens  can  probably  be 
explained  in  terms  of  the  process  zone.  Fatigue  precracking 
introduces  a  process  zone  ahead  of  the  actual  crack  tip.  Razor 
notching  does  not  do  so  and  thus  razor  notched  specimens  require 
additional  work  input  to  cause  voiding  and  fibrillation  prior  to 
crack  advance. 

The  suitability  of  razor  notching  vs.  fatigue  precracking 
for  these  materials  can  be  determined  only  on  the  basis  of  the 
actual  nature  of  expected  service  damage.  The  observed  decrease 
in  process  zone  length  with  decreasing  rate  would  suggest  that  in 
the  case  of  long  term  static  field  failures  a  small  process  zone 
is  to  be  expected.  Thus  razor  notching  which  produces  virtually 
no  Dugdale  type  zone  would  be  most  appropriate.  However,  if  the 
service  load  will  be  in  fatigue  or  a  mixture  of  fatigue  and 
static  loading,  then  fatigue  precracking  is  appropriate. 

Thickness-  The  results  of  the  thickness  comparisons  would 
seem  to  indicate  that  very  thin  specimens  can  still  yield  valid 
J's  if  they  are  sidegrooved.  The  amount  of  tunneling  observed  on 
non-s idegrooved  CT  specimens  indicates  that  sidegrooving  is 
essential  for  J  integral  testing  of  PE  when 


specimen  thickness  is  approximately  0.5  cm. 

Orientation-The  20  percent  increase  in  Jlc  which  was 
observed  when  specimen  orientation  was  changed  from  the 
longitudinal  to  the  transverse  pipe  direction  clearly 
demonstrated  the  existence  of  residual  alignment  effects.  Since 
extrusion  will  tend  to  align  polymer  chains  parallel  to  the  axis 
of  the  pipe,  crack  growth  perpendicular  to  this  axis  should 
encounter  greater  resistance.  This  was  the  case. 

Blunting  Line-The  pre-initiation  blunting  point  measured  for 
the  high  rate  test  clearly  indicates  that  some  account  must  be 
taken  of  the  strain  rate  dependence  of  the  yield  stress  of  the 
material  when  the  blunting  line  is  calculated.  If  the  appropriate 
yield  stress  is  used  the  blunting  line  will  more  accurately 
reflect  the  actual  behavior  of  the  specimen. 
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Conclusions 

Annealing  of  pipe  materials  prior  to  testing  caused  a  47 
percent  decrease  in  the  measured  fracture  toughness. 

Fatigue  precracking  yields  a  lower  J-,  toughness  than  does 
razor  notching.  This  drop  is  probably  due  to  the  presence  of  a 
preformed  process  zone  in  the  fatigued  specimen. 

Specimen  orientation  in  unannealed  material  influences  the 
measured  J-,c.  Longitudinal  specimens  yield  a  higher  apparent 
toughness  than  do  transverse  specimens. 

Fractography  indicates  that  CTOD  increases  with  increasing 
rate.  This  increase  is  reflected  m  an  increase  in  the  measured 
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TABLE 


SUMMARY  OF  Jlc  TEST  RESULTS 


Jlc 

kN/m 

dJ/da 

kN/m^ 

25 ( Jlc)/Sy 

cm 

Fatigue  Precrack 

unanneaied 

3.08 

74 

0 .41 

annealed 

1.65 

35 

0.216 

Razor  Notched 

unannealed 

4.67 

44 

0 .61 

annealed 

2.43 

35 

0.33 

annealed  high  rate 

6 .35 

43 

0 .84 

annealed  low  rate 

1.94 

12 

0.25 

longitudinal 

5.25 

0.69 

TABLE  :: 

PROCESS  ZONE  LENGTHS 

RATE _ PROCESS  ZONE  LENGTH _ VOID  DIAMETER 

low  (0.00254  cm/mi n)  0.266mm  13  Mm 

0 . 5  4  mm  1 9  Mm 


medium  (0.0254) 
high  (0.127) 


1.23  mm 


131  Mm 


T 
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ABSTRACT: 

The  mode  II  delamination  fracture  toughness  of  a  ductile  and  a  brittle 
unidirectional  graphite/epoxy  composite  has  been  studied  using  the  end-notched 
flexure  (ENF)  test  and  the  end-loaded  split  laminate  (ELS)  test.  The  stress 
field  in  the  vicinity  of  the  crack  tip  of  a  split  laminate  beam  under  mode  II 
(ELS  test)  and  mode  I  conditions  has  been  determined  by  means  of  a  finite 
element  analyses.  Also,  the  micromechanics  of  mode  II  fracture  have  been 
studied  during  in-situ  and  post-mortem  observations  of  the  fracture  process. 
Both  the  ENF  and  ELS  test  give  similar  values  for  GIIc.  However,  since  elastic 
material  behavior  is  assumed  in  the  analysis,  the  Gjjc  results  for  a  ductile 
composite  are  somewhat  uncertain  because  a  small  permanent  deformation  was 
observed.  The  ELS  test  has  been  found  to  provide  a  pure  shear  stress  state  in 
the  vicinity  of  the  crack  tip.  The  formation  of  hackles  provides  a  more 
tortuous  path  for  the  crack  'leading  to  an  increased  resistance  to  delamination 
under  mode  II  conditions  compared  to  mode  I  for  brittle  composites.  However, 
extensive  resin  deformation  and  yielding  p’ay  a  more  significant  role  in  the 
fracture  resistance  of  tough  epoxy  composites  for  mode  II  loading. 


INTRODUCTION 


One  of  the  most  important  factors  Hindering  the  use  of  fiber  reinforced 
composite  materials  in  structural  applications  is  their  inherently  poor  damage 
tolerance  for  deamination.  The  resistance  of  composites  to  Gelamination  can 
be  well  character' zee  by  tne  de'amination  fracture  toughness,  measured  as 
energy  dissipated  oer  unit  area  of  crack  growth.  Where  behavior  is  essentially 
linearly  e’astic,  tne  fracture  tcugnness  can  be  measured  as  an  energy  release 
rate,  G,  of  the  materia1.  Delaminaticn  in  composites  can  occur  due  to  tensile 
stresses  (mode  I),  in-plane  shear  stresses  (mode  II),  and  tearing  stresses 
(mode  III).  As  a  result,  a  compete  understanding  of  these  delamination 
processes  is  neeceo  to  properly  design  composite  structures  and  develop 
material  with  improved  fracture  toughness  characteristics.  Several  areas  must 
then  be  accressec  i~  c'de1'  to  throucr’y  understand  delamination.  First,  a 
critical  assessment  of  current  methods  usee  to  measure  resistance  to 
del  ami  nation  of  composites  's  needed  to  verify  the  validity  of  the  analysis 
for  the  various  test  configurations.  Furthermore ,  tne  geometry  independence  of 
the  measured  energy  re'ease  rates  must  be  established  by  comparing  measured 
va’ues  f or  different  test  geemetr-es.  r''n’te  e’ement  analyses  may  be  used  to 
determine  the  state  of  stress  at  tHe  crack  tip  where  beam  theory  does  not  give 
details  about  the  stress  field.  Second,  a  study  of  the  physical  mechanisms 
involved  during  the  fracture  process  must  be  undertaken  by  means  of  in-situ 
and  post-mortem  fracture  observations.  Investigations  of  this  nature  should 
help  in  the  understanding  of  the  re1 at ionshi p  between  the  macroscopic 
resistance  to  crack  propagation  and  the  micromechanics  of  fracture. 

Recent  investigations  indicate  that  increasing  mode  II  loading, 
particularly  for  brittle  systems,  results  in  a  significantly  greater 
resistance  to  crack  propagation  as  measured  by  the  total  energy  release  rate 
required  to  propagate  the  crack  [!’.  This  indicates  the  need  to  throughly 
understand  mode  II  delamination  in  order  to  properly  characterize  this  mode  of 


fa i lure. 


Several  test  conf iaurat'crs  have  been  proposed  for  measuring  tre 
resistance  to  delamination  unde*-  mode  II  loading  including  the  end-notched 
flexure  (ENF)  test  and  the  end-loaded  split  laminate  (ELS)  test.  The  enc- 
notched  flexure  test  configuration  can  be  seen  in  Figure  la.  This  test 
consists  of  a  three  point  bend  loading  configuration  by  which  a  split  laminate 
beam  specimen  is  de1  am1' rated.  The  (ELS)  test  configuration  snown  in  Figure  Id 
also  uses  a  spin  laminate  beam  specimen  loaded  at  the  clacked  end  and  fixed 
at  the  opposite  end.  Both  the  ENF  ana  ELS  test  methods  evaluate  GT^  frcm 
equations  derived  using  linear  beam  theory  and  assume  linear  elastic  mate^al 
behavior  [2,3!.  A  current  review  of  the  literature  indicates  t*at  bctn  test 
configurations  are  currently  being  used  to  measure  G>*c  .12,3,4!. 

During  a  recent  finite  e1emer,t  analysis  [5!,  the  ENF  test  was  found  to  be 
a  pure  mode  II  fracture  test  within  the  constraints  of  small  deflection 
theory.  The  study  also  revea’ed  that  the  ’’"ter  laminar  normal  stress  is 
identically  zero  along  the  beam  center  ’ire  a"d  the  inter1 aminar  shear  stress 
exhibits  the  expected  singu’arity.  No  ana'ysis  of  this  nature  has  yet  been 
done  on  the  ELS  test  configurator. 

As  pointed  out  earlier,  under  increasing  mode  II  loading,  particularly 
for  more  brittle  systems,  a  significantly  greater  resistance  to  crack 
propagation  has  been  observed.  Post-mortem  observation  studies  in  the  scanning 
electron  microscopy  indicate  that  this  increase  t  fracture  toughness  is  the 
result  of  the  development  of  hack'es  in  the  fractured  surfaces  (!'.  Therefore, 
the  hackle  formation  process  as  we11  as  damage  zone  processes  ahead  of  the 
crack  tip  need  further  study  by  means  of  additional  in-s’tu  and  post-mortem 
delamination  observations. 

The  objectives  of  this  investigation  have  been  (1)  to  measure  the  mode  II 
delamination  fracture  toughness  of  a  ductile  and  a  brittle  graphite/epoxy 
composite  using  the  (ENF)  test  and  the  (ELS)  test  and  compare  the  results,  (2) 


to  determine  tne  stress  fie'd  in  the  vanity  ot  the  crack  t;p  of  the  (E-Si 
specimen  under  noce  II  and  mcce  I  loading  conditions  by  means  ot  a  '^’te 
element  analysis,  and  (3)  to  study  the  micromechanics  of  fracture  us;ng  rea'- 
time  and  post-mortem  observations  oe  the  fracture  process  and  come' ate  t^e^ 
to  the  macroscooica' ’>  obtained  energy  reiea5e  r-ates  fo r  a  bctt'e  a^d  d-ct;'e 
compos i te. 

ANALYSIS 

nC- No* CU9C  ■  '  ^vj^p  /  £  S'  *  ' 

The  EN-  test  consists  c*  'ozc^c  a  sp';t  ’ami  mate  beam  specimen  using  a 
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2  b  (  2LJ+  3a  ) 

where  a  -is  the  i^tia1  crack  length,  b  1-s  the  width,  C  is  the  measured 
compliance,  L  is  ha’f-span  length,  and  P  is  the  applied  load  at  the  center 
pin.  The  strain  energy  re’ease  G y  t  can  a'so  be  obtained  by  means  of  an 
experimental  method  for  determining  the  re’ationship  between  compliance  and 
crack  length.  In  this  case,  G-j  is  given  by  (6! 

3  P  S  a  d 

Gtt  =  - - - O-  (2) 

“  2  t> 


where  m  is  obtained  •■"rcm 


C  =  ma0  +  const. 


(3) 


This  expression  is  the  result  of  a  least  squares  fit  of  experimental  data  from 
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data  ’s  octa^ec  frpm  'cac-ce'' 'ect  :on  c-^ves  at  d-'<'ere,~t  c-a  c«  'e^ct"s  zf  ore 

of  tre  tested  cc-ocr's  z*  *’dtr  b„. 

*»/ 

Er'd-Loadec  Soi:t  -a-j-’ate  [ _i_$2  "est 

The  energy  re'ease  ’•ate  *or  o-^e  "cce  ”  t~!s  case  ;s  determined  by 
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13!.  From 

this 

ana'>'s;s,  G't  car  be  exp-essea  by 


P  is  the  load  aoo’'ed  at  tre  c-acxed  e^c  zf  the  bean  (see  Hgure  4) .  The  axia' 
modulus  Ej’i  car  a’so  be  shown  to  be 


where  i  is  the  teta'  bean  t'?  deflection. 

Area  Metrpc 

In  a'l  the  methods  of  data  reduction  previously  described,  the  critical 
energy  re7ease  rate  was  calculated  by  knowledge  of  the  crit’cal  loac  Pc  at 
onset  of  crack  growth,  and  tre  crack  length  'a1  pr-or  to  crack  extension. 
Therefore,  G  j  r  c  7s  a  measure  of  the  energy  required  to  init;ate  crack 
propagat’d.  However,  as  the  crack  propagates,  G T r c  does  not  necessarily 
rema;n  constant,  in  *h;cn  case  a  propagation  and  arrest  va’ue  of  Gjjc  can 
ar;se.  Ar  alternate  approach  for  evaluating  the  fracture  toughness  is  the  area 
method.  In  th;s  case,  GTr  ’s  interpreted  as  the  energy  reqUired  to  create  a 
new  cracked  surface  area.  Grr  can  then  be  civen  by 
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(6) 


u  ' s  the  area  between  tne  1cad-def 'ect’on  curve  for  ’oading  and  unloading  for 
a  very  sma’’  change  of  crack  growth  (see  cigure  2).  An  important 

acvantace  cf  t^s  ~ef ’cc  ;s  tnat  on'y  e-ast’C  material  behavior  is  required  to 
predict  3t  t  .  There*cre,  for  geomef^ ca' ’y  non-linear  ano/or  non- linear  elastic 
material  responses,  this  netncc  g;ves  an  average  energy  release  rate  for  the 
observed  crack  extension,  -o'*  u^stab’e  crack  growth,  it  g’ves  an  average  value 
for  Gpc  «h^ch  tyo'ca','y  *a’  's  between  the  initiation  and  arrest  values 
measured  for  l;rear  be*'av'or'.  "^e  ’cad-displacement  record  should  return  to 
the  origin  to  guarantee  tnat  ~c  sigr;Pcant  far  f i e 1 d  damage  is  included  in  U 
in  Equation  (6).  T"e  irc'-sPr  c*  energy  dissipation  in  far  field  damage  in 
the  fracture  energy  ter-  „  ;n  Ec.  6  -cu'd  give  an  erroneously  high  estimate  of 
G m c •  f or  mafer;a's  wet  Z- -  ;s  ; '-ceoencert  of  crack  growth  rate,  and  crack 
growth  distance  (i.e.,  systems  Pth  minimal  fiber  bridging  and/or  platic 
wake),  the  average  ano  initiav'cr  G--  shculd  be  identical.  Furthermore,  where 
s 1  ow  stab’e  crack  growth  occurs,  the  average  G j j c  calculated  from  the  area 
methed  shcu’d  eq^a’  the  3-rc  for  ^PaPon  calculated  from  linear  beam 
theory.  Another  advantage  oc  this  metHod  is  that  Gj*  is  obtained  without 
knowledge  of  the  materia’ 's  elastic  properties  and  it  is  always  given  by 
Equafon  (6),  nega*-d ’ess  of  the  test  configuration  used. 

EXPERIMENTAL  PROCEDURE 

So’ it  laminate  beam  specimens  approximately  1:52.4  mm.  long,  25.4  mm. 
w;de,  and  4. S3  mm.  thick  for  the  end-notched  (ENF)  test  and  292.1  mm. long, 
25. 4  mm.  wide,  and  3.05  mm.  th;ck  for  the  end-’oaded  so’ it  laminate  (ELS)  test 
were  cut  from  24  ply  precracked  laminated  panels  of  Hercules  AS4/3502  and 
uexce!  T6C145/F1S5  graphite/epoxy  composites.  The  starting  crack  was 
introduced  on  both  oaneis  by  means  of  a  tePon  insert  0.127  mm.  thick. 
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Aluminum  and  brass  tabs  were  bcrced  to  tre  cracKec  end  of  the  ELS  test 
specimens  with  a  mixture  of  EA9C1  sfuct.ra-  ac^es^ve  and  31  curing  agent.  An 
UN -10- A  gray  brittle  coating  was  a'so  applied  to  the  lateral  sides  of  the 
specimen  to  fac’"itate  mon’tor’ng  crack  1ccat;ons. 

End-Notched  Flexure  "est'^c 

r;ve  spec'mens  Cc r  eac1"  g,-aph'te/epc.xy  system  were  tested  using  a  tnree 
po’nt  bend  fixture  *;th  a  tota’  1engtn  c-‘  101.6  mm.  attached  to  a  c’csed  loco 
servc-hycrau 1 ’ c  test  machine  .  5ri:r  to  testing,  the  specimens  were  pre- 
cracked  in  order  to  prcv;ce  a  sna^p  c^ack  t'p  frpm  *n’cn  to  initiate  tre  mode 
II  fracture.  _oac  co;nt  c’sp'acerne°ts  were  measured  by  tre  r‘am  d'sp'acement 
and  ^ads  were  me"- tc-ec  by  a  2225  *«.  'cac  ze". .  2ea’  time  ana’cg  pipes  o*'  tne 
1cad-def  1  ect :or  curve  we^e  "ace  on  an  x-y  ,~eczr'zer.  Tne  tests  were  ccrcucteo 
by  pcsit;cn'nc  tre  stec^e^s  ;n  tne  C'ree-oe  '"t  bene  fixture  sue-  tr.at  tre 
location  c*  t"e  c-ack  t'o  ; s  app'-c.x<Tate’y  T'cway  between  tre  center  a^c  ere 
o*  the  outs’*  icac*"gs  P'"S.  ."cer  z^z'  ace^e,-z  centre’  cone* tions  and  a  ram 
rate  of  2.5A  mm/rrn,  the  seec^e^s  *e-e  1eadec  unti'  de’aT;rat'on  c^ack  growth 
occurred.  ^ken,  tne  spec;me,'s  we^e  u^'eaded. 

Since  machine  st*"e<9  d ; sp’  aceme"t  was  used  to  measure  specimen 

deflection,  a  mac";r*e  cc*-p';aree  exper-'~ert  was  performed  to  make  the 

necessary  deflection  corrections  to  the  data  obtained  during  the  mode  II 

delamination  tests. 

Gjrc  calculations  were  mace  using  Equation  (1)  and  measured  compliances 
after  appropriate  riach-'ne  ccmp1ia”ce  corrections  had  been  made.  It  should  be 
noted  that  corrections  associated  with  shear  compliance  were  not  taken  into 
consideration  because  no  re’iao’e  values  of  the  shear  modulus  were 


available.  However, 

these  are  usually  sma11 

and 

this  omission  was 

not 

thought  to  have  a 

significant  effect  on 

tne 

measured  Gjjc  values. 

The 

measured  compliances  were  obta;ned  by  curve  fitting  tne  best  straight  line  to 
the  data  during  the  loading  part  of  the  test.  Equations  (2)  and  (3)  were  also 


used  to  calculate  Gi^.  'ro  use  these  two  ecuatso"S,  a  cc^pl  -'a^ce  ca'  -prat’ 
was  cerfonred  w’th  one  of  the  tested  specimens.  ri,-a'  ,  3t-_  va ' -es 

ca’cu’ated  by  the  area  method  were  dete^mi^eo  by  rears  c’  Ec^at'cr  '  t) . 
energy  absorbed  ;n  the  creation  of  a  new  cracxec  s-r^ace  *as  approx ^rateo 
numerica^y  by  integrating  the  area  under  the  loac/de* lect ;bn  curves  obtained. 
£nd- Loaded  Sol't  .arr'^ate  rest':nQ 

Two  soec'>e,'s  *or  e&c"  gnat"’’  te/eocxy  system  cnosen  f;r  tvs  study  were 
tested  ’ n  a  mode  I/mode  II  ^ixtu^e  attached  to  a  closed  'ccd  servc-nycrau 1 ;c 
machine.  A  rede  I  erec^ac*  was  ;ntrcducec  to  the  spec^e^s  to  sharpen  tne 
i n i t i a 1  crack  cheated  by  tne  tefcn  insert  a"d  to  prov'ce  an  i  ni  ti al  crac< 
"ength  to  beam  ’e^ct”  (a/-)  '•at'O  c’  approx ' rate 'y  0.55,  wre^-a  stab'e  crac< 
growth  i s  expectec  7  :  . 

Mcde  II  ce'  am;"at  ;cn  testing  *as  dc^e  unce’*  c;  sp"  acerent  contra' 
conditions  at  a  ■'ate  of  10.16  mm/min  and  6.35  mm/min.  ;_oads  were  recorded  wi tn 
a  445  N  load  cs’ 1 .  O'so' ace~er'ts  ..ere  ~'cr' ’ tc^ec  u$;rg  the  n=m  c ; so' acement  of 
the  test  raen’^e.  ’-e  tests  were  conducted  by  leading  tne  specimens  u"til  the 
crack  was  avowed  to  grow  approx -'mate’y  13.16  mm  to  38.1  mm.  At  that  po;nt  the 
test  was  stopped  and  tne  new  crack  location  marked  a^d  recorded.  This 
procedure  was  repeated  seve’-a1  f'nes  with  the  same  specimen  until  the  crack 
t’P  was  approx imate’y  12.7  to  25.4  mm  from  the  uncracked  end. 

The  critical  energy  reiease  rate  G  j  t  c  was  calculated  by  means  of 
Eauaticns  (4)  and  (6).  The  axial  modulus  was  ca1cu',ated  from  the  Icad- 
def’ection  test  data  anc  Equation  (5).  The  energy  absorbed  in  the  creation  of 
a  new  cracked  surface  was  approximated  numerically,  by  integrating  the  area 
unden  each  loop  ;n  the  ’oad-dei7 lect;on  curves. 

Finite  E  ’’e^gnt  Ana1ys  ’  s 

The  stress  field  in  the  vicinity  of  the  crack  tip  of  a  split  laminate 
beam  specimen  unden  mode  II  and  mode  I  was  determined  using  a  finite  element 
analysis.  This  was  a  linear  analyst  made  with  a  finite  element  algorithm 
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deve’oced  by  Herr’xsen  [3;.  The  a'gor^t^n  's  based  on  a  nonlinear  coce  uocatec 
«'  t*  Lac^ana'  a”  fcrmu’  at  ion  using  s'x  and/or  e’ght  nodec  ’’ soparaxetr  ic 
e’emerts  w’-t"  t*c  degrees  of  ‘reecom  per  °ode. 

To  overcome  the  difficulty  ’mposec  by  tne  singu'arity  at  t"e  crac<  v:p,  a 
sobstant’a'  ’•e* 'rement  of  the  mesh  i-n  the  vicinity  of  the  crack  tip,  ana  the 
use  of  tr-'angu’ a**  s;x  noded  elements  arounc  tne  crack  t’o  w;t"  m’d-s^de  ncces 
a ; sp 1 aced  to  t^e  c-a^te^  po'^t  n as  oee~  co^e.  ’"ese  crac<  t;p  e'ements  contain 
a  s^gu’ar'ty  c*  tne  ‘orm  l/.*-  ,  prov'c'ng  a  stress  fie’a  wh-'cn  agrees  w;th 
the  theoret’ca'  stress  s;rgu’ar’'ty  c*  ’'near  fracture  mechanics.  These 
spec'a'’  ncdes  a ’sc  co^ta'n  rigic  body  action  and  constant  strain  modes, 
thus  sat’S'y1'r'G  a''  tne  "ecessa-y  ccrc't'cr,s  ‘or  monotonic  convergence  ;  9 ; . 

The  generated  ~esh  cprs;sts  cf  791  rpces  and  294  e’ements  ana  can  be  seen 
'n  Figure  3.  vcce  1  'oacng  *as  s^j'atec  by  app’y-'ng  a  symmetric  'oac  at  tne 
crac<ed  end  o*  tne  -esn.  Voce  II  ’oac^g  was  ’ntrccuced  by  asymmetrical  ly 
1oadinc  the  enacted  e~d  (E_S  test,'. 

The  materia’  propert'es  usee  i-i  t^s  ‘;r;te  e'ement  analjS’S  *.ere  from 
typica1  properties  c*  un'c'rect'ona1  AS-/35C2  graphite/epoxy  composite.  Linear 
beam  theory  as  we’’  as  ’'"ear  elastic  material  behavior  were  ass-meo  'n  the 
ana'ys's.  Stress  contour  p'ets  were  obtained  from  the  output  data  and  tne 
stress  distribution  was  o’otted  as  a  ‘unction  of  distance  aheac  of  the 
crack  tip. 

Scanning  E’ectrpr  Microscopy 

Sp’it  laminate  specimens  approximate 'y  30.5  mm.  long,  5.1  mm.  w’oe,  and 
1.02  mm.  thick  of  Hercu'es  AS4/35C2,  a"d  30.5  mm.  long,  5.1  mm.  wide  and  1.27 
mm.  th^ck  of  Hexcel  T6C145/F 1S5  graphite/epoxy  composites  were  cut  from  3  ply 
and  6  ply  laminated  panels  of  these  two  materials.  Specimen  dimensions  were 
determined  by  the  space  allowed  by  the  tensile  stage  of  the  scanning  electron 
microscope,  and  stab'e  crack  growth  cr^ter^a  [61.  All  specimens  were  polished 
with  0.3  alumina  powder  and  cleaned  u’trasonical ly  in  a  scanning  electron 
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microscope  "rec'1  113  so’ut'on  nam’ess  to  the  i"es’n. 

‘■'ode  II  observations  of  the  de1  aminat’cn  fracture  process  cere  aeterm':red 
us’pc  rea'!-time  cose^vat 'C^s  4r  the  scar,°;r'c  e’ectrcn  microscope.  Mode  II 
ce 1  ami nati or.  was  achieved  by  means  of  a  esoec'a’,y  designed  three  point  oerc 
fixture  instal 1ec  to  the  tens’ le  stage  o*  t"e  scanning  electron  microscope,  ~o 
confirm  that  these  real -t’me  observations  «ere  representative  of  the  bu'< 

u6  ^  i  n<J  ^  '  QH  S’3*  dV  '  O  t*  CO  $  *  —  **  6  ^  Stud  4  6  S  rt'0^0  ^0  »*  +  Q  r'r71Gd 

A' I  surfaces  were  coated  with  a  150  angstroms  thick  gold/pal lacium  f  1 In 
to  minimize  charging  ejects  associated  with  the  noncor.ducti ve  nature  of  tne 
epoxy  resins  of  tne  compos''tes .  T^e  exoe^;mer,ta’  results  were  recorded  on  both 
v i ceo  tape  and  an  standard  t^-t  f;,m. 

RESULTS  AND  DISCUSSION 
End  Nctched  ~ 1  exur g  ( E nc )  'rest  ^esu 1 ts 

Tab’e  1  s^Ows  tne  ENF  test  resu'ts  for  both  grapni te/ecoxy  systems 
invest;gated.  wet”cc  I  corresponds  tc  3r-c  values  obtained  using  beam  theory 
equations  (Equat^C”  (1))  anc  measured  cq^o ' ' ances.  Method  II  results  are  based 
on  the  compliance  calibration  (Ecuatic^s  (2)  and  (3)).  Method  3  refers  to  cata 
reduct;cn  based  on  the  area  method  (Equation  (6))  and  the  energy  dissipated 
under  the  load-de* ’ect ion  curve.  QeTanrnation  crack  growth  occurred  at  typical 
load  values  of  556  to  612.5  n  fcr  AS-/2532  graphite/epoxy  and  was  unstable  fen 
this  material.  rigure  4  shews  a  typica1  load-deflection  curve  of  AS4/3502 
graphite  epoxy  sowing  the  sudden  drop  of  load  as  unstable  delamiration 
occurred.  This  ind;cates  that  the  Gjrc  values  measured  using  the  area  method 
wil1  net  only  ;nc'ude  the  energy  ’'ecu’'r‘ed  to,  create  a  crack  surface  but 
will  a’so  inc'ude  energy  dissipated  ;n  dynamic  crack  growth  processes.  The 
very  sma-1  nor’ J "earity  just  prior  tc  crack  extension  <s  prebab’y  associated 
w;th  the  development  oe  a  crack  tip  damage  zone. 

"he  average  c^fcd’  energy  re’ease  wate  for  AS4 /3 502  using  methods 

I,  II,  and  III  oe  data  reduction  were  560  J/m*- ,  648  J/m",  and  613  J/nO, 


respective’^,  "fese  values  are  a’ '  w''tn’n  15  percent  c'f^ererce  of  eacr  otter, 
wtvch  can  oe  cons ice’-ed  withi"  tne  expe^imenta'  e^^or  associated  w’tn  tne 
ff event  methods  cf  data  reduct  <c'.  >  Cact,  ‘•‘ur-i  and  O'Brien  nave  i ndicatec 
that  Grrc  measurements  may  vary  «Jt"  method  cf  data  reduction  up  to  9«  *itn 
no  consistent  pattern  ;n  tne  var-acc"  [  1C j . 

Figure  5  stows  a  typ'ca’  oad-cef iect’on  curve  obtainec  'or  ncce  II 

de ’ ami  nat ;  or  testae  of  T6Cl-5/rl£5  t,,aj"ite/e30xj .  CracK  gvcwtn  was  stao'e  as 
seen  ’ n  tne  upper  pa^t  zf  tne  c-rve.  -clever,  tne  po’nt  at  wnich  era c<  growtn 
occuvrec  ;  s  "ot  we’’  aer*'nec .  ^"e'S'ce,  tne  Tax’ mum  1  oad  reacnec  curing  eacn 
test  was  used  to  estimate  t*e  c'set  :f  c' act  cv'Cwtr'.  ”  n  e  s  e  va'ues  we^s  oetween 
1390  ano  1631  S'.  Tne  curve  a's:  '"e'eates  an  ’"e’ast-c  serav^cr  of  tne 
nater  ’  a  ’  t"e  -r'’oad  c-r  *e  - 1  o  -ct  '■“ty"  to  tue  cr;t'  n  1 .  ’’’hpre'pre,  a'l 

tne  assumptions  o'  tne  test  -ave  tee"  «;o'ateo.  It  snc-'d  be  noted  that  tne 
’oad-cef ’ect icn  curves  at  c:'ce'rer'Z  c^acv  'edge's  usee  fcr  tne  copp’iance 
cal1' oration  a  ’so  snowed  "on- ’  i”ea’"  oe-av  ;or.  7-e  Gr  »c  values  predicted  by 

methods  I,  II,  a^c  III  a^e  a’’  a"t'o-c.,s.  -c-ever,  the  energy  ■•e'ease  rate 

obta;ned  by  tne  a^ea  met "od  gives  a  better  approx imat' on  of  tne  critical 

energy  re'ease  rate  G  r  t  c  because  at  ''east  tne  geometrically  non-Iirear 
behavior  of  tne  specimen  ’'s  taxer  ;^to  cons; derat i on  by  this  metneo.  The 

ca1cu1ated  average  G  •  m  cf  Idle  j/~-  sno^’c  oe  considered  an  upper  bound 
estimate  of  the  mode  II  critical  energy  release  rate  of  T6Cld5/F135  because 
the  inelastic  behdv;or  is  no t  acccu”tec  for  in  this  ana'ysis.  A  J-;ntegral 

approach  needs  to  be  deve’eped  to  obta;n  a  more  meaningful  measure  of  the  moce 
II  d6 1 3^  i  nd  t ;  on  -C'jchr,0ss  c  *  t^s  system. 

End-Loaded  Split  ^aminate  (E~S)  Test  Fesu'ts 

The  ELS  test  results  are  seen  in  'abie  2.  Method  I  corresponds  to  Gmc 

calculated  us’'ng  Eauation  (d)  a^d  w;th  l,,  ca’culated  from  Equation  (5)  ard 

test  data,  ^etnod  2  corresponds  to  the  a^ea  ■"etHcd  (Equation  (6)).  An  average 
Em  of  125.5  GPa  for  ASd/35C2  and  109. 9  Gpa  for  T 6C 1 4 5/ F 1 S5  were  calculated. 


cigure  6  shows  a  typica'  'cac-cef  1ect'or>  curve  cbta'"ed  a* ter  testae  ere  c 
the  AS4/3502  soecrers.  As  see"  'n  tn;s  f igure,  unscab'e  crac<  growth  ;s 
observed  cor  the  * ; r s t  crack  extension  oesp;te  t-'e  (a/L)  rat(c  z*  at  least 
0.55  i  nt  reduced  to  a''  E.S  seec'^e^s  cc’  ~g  the  node  I  precrack  cc^e  pr5or  to 
nece  II  test'rg,  wn’c"  according  to  the  stability  cons ’cerat 'cns  cf  tn’s 
leading  arc  geometry  sneu'e  have  been  stab’e  )3).  This  is  probably  t"e  result 
of  the  acd't;cr<i'  e-'e-gy  storage  ’n  tne  ♦ 'xtures.  s'cn1  inear’ ty,  escec'a1 ’y  at 


1  so  1  scene  nt  s  oroyser  t~ar 


ray  a’sc  be  associated  wit.n  f 


lim’tat;crs  cf  ''"ear  bean  tnec-y  usee  ’"  tne  stab^ity  analysis.  An  average 


GtV-  c f  59: 


“i  /m*-  ,  ""dbT 


:d  1}  and  593  Z/'T'<~  (method  II)  were  founc  for  c. s  —  /2 


nese  res:  ts  a^e  *- 


d’  '*ererce  cf  each  ot^ec  -cwever. 


co"s  ;de"ao'e  scatter  ,s5S  observe:  t”e  G-*.  data  from  eacn  spec 'ten  (Tao'e 

2'1.  ~n  ■  s  scatter  -s  be'  ;evec  to  be  t*e  resu't  o*  tne  n-'gh  degree  c  error 
•"vc’ved  -easur-rg  crack  ’ccar!cr's  w-iicn  is  dcne  by  monitoring  tne  crack  as 
:t  propagates.  >  co"trast,  tn;s  scatter  was  very  sna’’  in  the  ENF  resjlts.  A 


n  ^  4  ■*  _  r*  a  *  p  ^  -  n  r->  (-  :  y*  \  a  *■  ' 


:s  seen  in  Ficure  7.  Stable 


crac<  g-owt-  ;s  observe:  *or  a’’  ’cad/c'cad  'occs.  >i$  is  evicencec  by  tne 
s-oct"  decease  'cad’"c  "ate  *c "owe:  by  a  decrease  ’n  ’oad  until  the  crack 
was  arrestee.  It  sncu'd  be  -otec  t-at  ;t  ’s  not  known  when  crack  growth 


Started  iu  '  ’  r  C  0-"'S 


is,  w"*ch  was  a  cb"rcn  urcerta’nty  for  this  material 


for  boC  types  o'  node  II  tests.  At  ’east  some  cf  the  observed  nonlinear 
be"av'or  ’ s  due  to  ;re’ast;c  te^av'c  s;"ce  the  tbe  curves  fail  to  returned  to 
the  or-g-r  as  t”e  speerne"  ; s  c'oaced.  Nor'  'near  se"av'or  associated  with  the 
' taf  b"S  o*  ’ ;  "ear  tnecry  are  a’so  observed  for  def1ections  orgater 


'he  average  "rcce  II  energy  release  r.a:e  va’ues  for  76C195/F135  were  found 

n  o 

to  be  22.-  5  2/^  (^et"od  I)  anc  2260  J/t  (•"et^oc  2).  However,  a  systematic  or 
mono  ton  - c  -'"crease  G--c  for  •"etnod  I  can  be  observed  throughout  the  test. 


's  ’  s  thp  ro  :  o*  0"e 


the  'oac/def lect ion  curve  seen  in 


Figure  7.  Although  tn;s  ron' ; near’ ty  ’ess  severe  than  the  nonl inearity 
observed  fo r  th;s  •Tiate''’a'  tested  us;rg  the  EN^  test,  •  t  is  sufficient  to 
invalidate  tne  data  analyses  done  usmc  ’inear  beam  theory  (method  I).  The 
scatter  of  the  ’•esu'ts  fo^  •"ethcd  '.1  (a^ea  method)  ;s  more  random  and  is 
associated  with  t"e  excen'menta’*  -'cu’-ty  in  obtaining  a  precise  measure  of 
the  change  in  crac<  length.  However,  these  va’ues  are  an  ove^est^ation  of 
Gt  t-  because  the  curves  did  "ct  return  to  the  crig;n  as  tne  specimens  were 
un’oadec . 

A  comparison  of  the  resu’ts  predated  using  both  test  ccnf igurations 
indicate  that  <rcr  «S^/35C2  (a  brittle  composite)  s’mi’ar  results  were  obtained 
for'  the  EN-  a np  E-S  tests.  In  contrast,  r'or  1601-15/^135  both  the  ENF  and  ELS 
test  do  ”ot  give  va’;c  nesu'ts  s’nce  tne  assumptions  assoc;atec  with  each  test 
method  wens  v'c'atec.  ~nen=""pre  t  ;  r  oncer  to  properly  cnaracten''ze  tne 
fracture  toughness  c’  tK’s  mater* al  and  toucner  grapni te/epoxy  composites  in 
ge"eni’,  an  analyses  wh-'ch  a”ows  Ccn  mater* a’  ana  geometric  noni;near’ty  sucn 
as  a  l-'ntegna'  ;s  needed.  It  is  worth  noting  that  s;mi'ar  resu’ts  were 
obtained  cor  t-^e  3 n -  a"c  E_S  «rer  cata  was  ana’ysec  us''ng  tne  area  method. 

However,  bet n  res.'ts  sneu'd  be  considered  upper  bc-nos  of  the  mode  II 

crit’ca1  energy  ro’ease  rate  zf  the  mater;a’. 

Fjn-'te  E’ement  *ra\,  s*’s 

Figure  S  shows  the  ''"-'te  e’emert  resu’ts  of  a  split  laminate  beam 
specimen  tested  uncer  cure  ■node  I  conditions  (dcub’e  cantilever  beam  test). 
rigure  Sa  is  a  stress  contour  of  tne  norxa’  (Syy)  stness.  A  perfectly 
symmetric  stress  condition  ;s  shewn  thncugn  the  beam  sickness  indicating  the 
exoected  mode  I  lcac;ng.  Tne  st-ess  ;s  seen  to  nap’a’y  'release  as  the  crack 

tip  is  approached.  Figure  3b  shows  the  Syy  stress  distribution  ahead  of  the 

crack  tip  at  the  midpiane  of  the  specimen.  The  Syy  stress  is  a  principal 
normal  stress  a’eng  this  p’ane  and  it  rapidly  drops  off  with  distance  from  the 
crack  tip  until  ’t  is  compressive  at  a  distance  of  0.76  mm  ahead  of  the  crack 


13 

t"D.  r-ina”>,  ’t  grac-a”y  approaches  a  zero  stress  'eve1  at  approx -irately  3 
mm  fro«i  pre  crdc<  t;p,  *nicn  -s  ma-ota^med  a’’  tne  way  to  toe  end  of  the  team. 
Simitar  perav-pr  *as  observed  when  tne  stress  ^-e’c  ahead  of  the  cracx  tip  was 
calcu’atec  by  "’ea-s  o;  a  --'g’-  o^ce"-  o’ate  tneo^y  |1I|.  An  ;dent’ca'  loading 
case  out  assuming  *sctropic  mate^a’  be"avior,  a'so  -'nd;cated  that  compressive 
stresses  deve’opec  *r»  the  specimen  anead  of  the  crac<  tip.  "his  suggests  that 
these  co!,’press’ve  stresses  ane  due  to  geo-etry  ana  'oading  rather  than 
materia’  arvsotnppy.  -'clever,  tn-s  compne^s^ve  stress  state  is  not  expected  to 
sigrv^-'cantly  e'-’ect  de’  am-ravon  because  -t  develops  far  enough  away  from  tne 
crack  tip.  The  das-ed  ' - ne  a’sc  seen  tn  --cure  So  -no-cates  the  magnitude  of 
t-e  snean  stness  Sx_,,  wr';c''  -s  seen  to  oe  zero  everywhere  atcng  the 
ce  ’  a"1- nat' pn  p'd-e  •c'  tne  mcce  1  ’cac^g.  ’"erefore,  this  conf-rms  that  a 
pu**e  ~-c  ce  1  ’cac'ng  mac-cscco- :a "  j  toes  -nceec  c;ve  a  pure  mode  I 
de’a~;nacc-  st-ess  mic,'osccp;ca,,y  ar:.-c  t-e  crack  tip. 

r;gce  9  s ”Ow s  tne  results  ce  a  so’;t  ’am- "ate  beam  spec-men  subjected  to 
an  asymmetric  ’cac;r-g  case  cf  pee  -cce  II  cc-d it’on.  Similar  to  the  normal 
stress  ccer  ~cce  I  co~c:  "i-  z-s ,  t-e  s-ear  st-ess  ;s  per^ect'y  symmetric  across 
tne  beam  thic<ress  suggest^  a  pee  moce  II  loading  condition  on  the 
speemen,  as  seen  in  ricure  9a.  uowever,  tKe  snape  of  the  stress  f i e 1  d  appears 
more  narrower  and  elongated  tKan  for  "-ode  I  loading.  Figure  9b  snows  the 
stress  f-!e!d  as  a  Arctic"  of  distance  ahead  of  the  crack:  t-'p.  In  th;s  case, 
the  shear  stress  can  be  seen  to  mcnotonica’ ’y  decrease  to  a  constant  value 
which  extends  all  the  way  to  the  end  of  the  beam.  The  dashed  line  ino;cates 
the  normal  stress  distribution  a’ong  the  delaminatic-n  plane  to  be  zero;  i.e., 
macroscopic  mode  II  ’oading  does  indeed  appear  to  give  a  microscop-'c  mode  II 
stress  along  the  delamination  plane.  As  can  be  seen,  the  shear  stress  ahead  of 
the  crack  tip  for  mode  II  loading  decays  much  more  slower  than  the  normal 
stress  for  mode  I  loading,  furthermore,  for  mode  II  loading  the  extent  of  the 
stress  field  above  and  below  the  delamination  plane  is  much  smaller  for  mode 


14 


II  than  for  mode  I,  as  seer  the  stress  contours  plots.  This  results  are 
consistent  with  the  SEy  observations  of  damage  zone  size  and  shape  to  be 
reported  in  the  next  section. 

Scanning  Electron  Microscopy  Observations  cf  Mode  II  de'anrnatlon 

Figure  IQ  corresponds  to  rea1  time  observations  under  mode  II 
delamination  of  AS4/3502  graphite/eooxy.  Figure  10a  shows  part  from  the  damage 
zone  ahead  of  the  crack  tip  which  developed  as  the  specimen  was  being  loaded 
(!eft  side  of  micrograph).  This  damage  zone  is  characterized  by  extensive 
microcracxing  and  it  was  found  to  oe  at  least  180  -m  long  and  extending  about 
5um  above  and  5  _m  below  the  o'are  of  deiamination,  and  resembling  the 
shape  of  a  Dugdaie  type  damage  zone,  long  and  narrow.  This  eiongated  damage 
zone  is  at  least  twice  as  long  as  the  damage  zone  reported  ahead  of  the  cracx 
tip  for  mode  I  delaminate  c*  the  sim^ar  material  AS4/3501-6  [12].  This  is 
probably  because  for  mode  II,  the  resin  rich  region  between  plies  behaves  like 
a  soft  material  between  rigid  platters  with  all  the  strain  being  localized  in 
th;s  region,  and  also  because  of  the  s'ower  rate  of  decay  of  the  stress  field 
under  mode  II  compared  to  mode  I  >"evea''ed  by  the  finite  element  analysis.  The 
hackee  formation  process  often  associated  with  mode  II  deiamination  of 
composites  made  with  relatively  britt'e  resins  can  also  be  seen  in  Figure  10. 
First,  microcracks  start  to  form  at  approximately  45  deg.  to  the  plane  of  the 
plies  (Figure  10a),  perpendicu lar  to  the  principal  normal  stress  direction. 
Then,  macroscopic  crack  extension  occurs  as  the  sigmoidal  shaped  microcracks 
beg;n  to  coalesce,  forming  the  hackles  (Figures  10a  and  10b).  During  this 
process,  some  resin  deformation  must  have  accompanied  the  coalescence  of 
mi crccrack s  into  hack’es  as  the  gap  between  adjacent  hackles  results  from 
rotat;on  at  the  base  when  microcrack  coalescence  takes  place.  Note  how  the 
fina’  hackle  or’entaicn  is  steeper  than  the  45  deg.  angle  of  orientation  of 
the  i^tM1  microcracks  ^rom  which  the  hackles  form.  No  clear  indication  of 
c^ack  growth  direction  is  given  by  the  hackle  orientation  on  the  fractured 
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surfaces.  Arrows  at  the  upper  left  hand  corner  of  all  micrographs  indicate 
crack  growth  direction,  figure  11  further  illustrates  the  difficulty  in  trying 
to  infer  crack  growth  direction  trom  hack’e  direction  alone  [13!. 

In-situ  deiamination  observations  of  T6C145/F185  graphite/epoxy  are  shown 
in  Figure  12,  where  the  damage  zo^e  deve^o^ed  anead  of  the  crack  tip  as  the 
material  was  being  loaded  is  seen.  This  damage  zone  was  approximately  1000 
•_m  long  and  extended  at  least  25  „m  above  and  25  _m  below  the  plane  c* 
delamination  near  the  macroscopic  crack  tip,  progressively  decreasing  to  about 
5 urn  at  the  end  of  the  damage  zone.  The  damage  zone  is  characterized  by 
extensive  microcracking.  It  is  a'lso  approximate^  15x  larger  by  area  than  tne 
damage  zone  for  AS4/3502,  indicating  that  much  more  energy  is  dissipated  in 
the  delamination  process  *or  this  tougher  materia’  (T6C145/F135) .  Note  also 
that  the  mi crccrac<s  do  "ot  cca’esce  to  form  hackles  as  the  macroscopic  crack 
is  created.  Instead,  the  resin  extensively  deforms  as  the  localized  strain 
increases  until  a  complete  seoaration  of  the  surfaces  occurs  due  to  local 
yielding  (Figure  12b).  Note  the  large  strain  undergone  by  the  resin  in  areas 
well  above  and  below  the  resin  rich  area  where  the  crack  formed.  The  extensive 
resin  deformation  and  large  damage  zone  observed  explains  the  non-linear 
behavior  observed  in  the  load-deflection  curves  for  this  material  (Figures  5a 
and  7).  Furthermore,  the  much  greater  resistance  to  delamination  of 

T6C145/F 185  compared  to  AS4/3502  measured  for  mode  II  loading  (2260  J/m^  vs. 

2 

540  J/m  )  is  consistent  with  these  observations.  The  long  narrow  shape  of  the 
damage  zone  for  both  systems  is  consistent  with  what  one  would  expect  based  on 
the  finite  element  analysis. 

Post-mortem  fractography  of  the  fractured  surface  under  mode  II 
delamination  of  AS4/3502  can  be  seen  in  Figure  13.  The  most  significant 
artifact  observed  in  these  micrographs  are  the  hackles  formed  on  the 
fractured  surface.  Note  the  feather  like  appearance  of  the  hackles  which  have 
deveiooed  at  a  very  steep  angle  (large’-  than  the  45  deg.  angle  at  which  the 


mi crocracks  were  observed  to  develop  during  the  in-situ  observations).  See 
Figure  13a.  Therefore,  the  hackles  undergo  some  rotation  before  full 
separation  of  the  fractured  surfaces  occurs.  This  is  in  good  agreement  with 
the  final  hackle  angle  orientation  from  the  in-situ  observations  (Figure  10) 
Figures  13b  and  13c  show  the  detailed  ^'veT'  pattern  markings  that  develop  on 
the  surface.  No  clear  indication  of  crack  growth  direction  is  seen  from  the 
river  pattern  markings  on  the  hackees  or  the  direction  in  which  the  hackles 
point  (Figure  10a) .  In  fact,  the  direction  the  hackees  and  the  river  pattern 
markings  point  wi  l1  be  determined  by  whether  microcrack  coalescence  occurs  on 
the  upper  or  lower  boundary  of  the  resin  rich  region  between  pl4es,  and  this 
is  a  random  process. 

Post-mortem  observat4cns  the  fractured  surface  under  mode  II 

delamination  of  T6C  1^5 /p  1S5  composite  ca"  oe  seen  4n  r;gur‘e  14.  Note  the 

extensive  rancom  defo^at io-'  jnce^go^e  by  t~e  res‘n  ideating  that  the  resin 

fared  due  to  yie’dirg  instead  c*  t-e  b^-tt’e  :_,:e  fa4'uro  facing  to  hackle 

formation  of  the  mere  o^tt'e  grapn;te  etc*.,  s.ste-'s  s-ch  as  AS4/3502. 
Correlating  these  f;gures  ^t"  =  ‘c-'e  12  t-e  •  -  -  .  - 1  -  observations  of  this 

material,  it  can  be  ,sr--':ec  t-a:  t-e  -■ c  ;  (''n-s’tu)  do  not 

general’y  coalesce  to  fo^m  racx’es.  3dt-er,  a  -m.:"  •ncr‘e  duct4  'e  ^ra1  fracture 
is  noted. 

As  ’ t  has  been  sHown  a’^eady,  the  cr;t;ca'  e"er:7  re'ease  rate  for 

T6C145/F185  is  appro*  :mate’y  five  t4re.  ’a--cer  ’md"  fo r  AS4/3502.  Therefore, 
it  can  be  inferred  that  a’trcug*  the  formation  of  nack'es  provides  a  more 
tortuous  path  for  tne  crack  ’eac4^  to  ar  ^crease  ;n  fractured  toughness  in 
brittle  composites,  extens’ve  resin  deformat;cn  arc  yie'ding  p’ay  a  more 
significant  ro’e  in  tke  fract-red  toughness  resistance  of  tougher 
graphite/epoxy  systems.  In  fact,  the  reported  fracture  toughness  of 
T6C145/F185  under  mode  I  conditions  where  little  to  no  hackle  formation  is 
observed  is  similar  to  the  fracture  toughness  of  the  material  under  mode  II 
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conditions 


CONCLUSIONS 


(1)  Both  the  E N F  and  ELS  test  give  similar  values  for  the  critical  mode  II 
energy  release  rate  G*jc.  However,  since  elastic  material  behavior  is  assumed 
in  the  analysis,  the  G*-c  results  for  the  ductile  ccmoosite  are  somewnat 
uncertain  because  a  sma1!  permanent  deformation  was  observed. 

(2)  Different  size  and  snaoe  of  damage  zone  observed  when  comparing  moce  I  to 
mode  II  ;s  cons;stent  w^tn  the  f^e^ce  ;n  the  size  and  shape  of  the  stress 
field  predicted  by  finite  element  analyst. 

(3)  The  damage  zone  deve1oped  a^eac  or'  the  crack  tip  at  the  onset  of  cracx 
growth  unde1"  mode  II  de’ arrrnat  ion  concisions  is  much  larger  by  area  for  the 
compos;te  made  using  a  more  cuctiie  epc.xy  (r135)  than  for  the  composite  made 
with  a  brittle  epoxy  { 3 5 C 2 ) . 

(4)  The  formation  of  hacxles  p^cvices  a  more  tortuous  path  for  the  crack 
iead;r,c  to  an  increased  resistance  to  de'amination  under  mode  II  conditions 
compared  to  mode  I  for  britt’e  composites.  However,  t"e  extensive  resin 
deformat’lon  and  yielding  plays  a  mere  significant  roie  in  the  fracture 
toughness  vesistance  o*  tougKer  composites  g;ving  essentia’iy  similar  fracture 
toughness  va’ues  under  mode  I  and  mode  II  conditions. 


T 
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Tab’e  1. 

-  Enc-notc.^ed 

-’exure  ( E NF )  Test 

Resul ts . 

r' 

<2  r  t 

c  (3/^) 

Mater -'a  I 

Soec inen 

I 

T  T 

T  T  r 

AS-V3502 

1 

525 

643 

£13 

2 

503 

630 

C  y  s 

■j 

560 

632 

613 

4 

5-3 

643 

630 

Av.  543 

648 

613 

★ 

T6C145/F 135 

1 

- 

- 

224C 

2 

- 

- 

2695 

3 

- 

- 

2485 

4 

- 

- 

1995 

Av. 

- 

2393 

*  No  ca'c-'at'cn  of  3^^.  *s  ?2sio'e  due  to  significant  ron linearity  in 
’cac/de*’ '  =  :t fon  curves .  T*  one  uses  linear  portion  of  curves  to  measure 
core  ’  i  arce  arc  P  'or  crac:<  extension  (ignoring  nonlinear  behavior  rear 
?!7iax},  Gj'C  lower  ocunc  va'ues  or  1225  nay  be  calculated. 
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Table  2.-  End-lcaoed  So' 

;t  Laminate 

(ELS)  Resu 

Its. 

Grv  ( 

J/pl*"  \ 

Materia'  Specimen 

r 

II 

AS4/35C2  1 

620 

cc  5 

C' 

513 

7? 

c^S 

t  /  7 

0 

59  5 

525 

595 

525 

5CS 

503 

Av. 

595 

595 

215  3 


2030 

2/pon 

2275 

2415 

0  C 

2363 

p '  ■:  c 

i  C7  : 

2235 

2113 

!  ££  0 

2643 

2223 

2242 

L  -v  -  0 

2033 

259C 

2 1  70 

2713 

2765 

2SGC 

2323 

2275 

2260 

*  Nonlinearity  is  too  great  to  sat-'s-'y  linear  oeam  tneory  assumptions  implicit 
ip  tfts 
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Figu-e  5 

Figure  6 

Figure  7 

Figure  S 
Figure  9 
Figure  1 
Figure  i: 
Figure  1 
Figure  1 
figure  1 


FIGURE  CAPTIONS 


-  Mode  II  test  configurations. 

-  Alternate  I "terpretat ;cn  c’  G. 

-  Finite  Element  Mesh. 

-  Tipica1  I oad-de*1 ect’ on  curve  of  AS4/3502  tested  using  the  ENF 
test. 


T^ca'  ,cad-d9f"’ec:icn  curve  of  T6C145/F135  tested  using  the  ENF 
test  conf ’gurat ;on. 

Tip;cal  1cac-deflection  curve  of  AS4/35Q2  tested  usin.a  the  ELS 


test  con^igurat ;on. 

-  ^ip-'ca1 

’cac-def'ec 

t  ’  C-r  curve  cf 

T6CI45/F1S5  tested 

test  co 

nf ’gurat ;cn. 

_  F  i n  ’  t  e 

6<1-6^6rt  '’IQdg 

I  results. 

-  r ;  n  ’  t  e 

0"10rn0n^  tic  C° 

II  resu'ts. 

voce  I 

I  de 

1  drvi ndt ^  Cn  C f 

AS  4/ 3502 . 

.  -  Moce  I 

I  in-s't-  ce 

"  arr4  ndt :  on  of 

AS 4/ 3501-6  !  I31. 

vcde  I 

I  in-s't-  de 

1  an' nat i c"  oc 

TSC1 A5/F 1S5. 

wode  II  pcst-mprter"  ‘ractcgrapry  o£  AS -4/3502 . 

I.-  Mcce  II  post-norfen  fractograp^y  o*  ”6C!45/F  135. 


figure  10 


Figure  11 


Figure  12 
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Figure  14 


